WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 6 
G06K 15/332 



Al 



(11) International Publication Number: 
(43) International Publication Date: 



WO 99/34322 

8 July 1999 (08.07.99) 



(21) International Application Number: PCT/US98/27624 

(22) International Filing Date: 23 December 1998 (23.12.98) 



(30) Priority Data: 

60/068,834 



24 December 1 997 (24. 1 2.97) US 



(71)(72) Applicant and Inventor: MILLS, Randell, L. [US/US]; 
1780 Greebriar Court, Yardley, PA 19067 (US). 

(74) Agent: MONACO, Daniel, A.; Seidel, Gonda, Lavorgna 
& Monaco, P.C., Suite 1800, Two Perm Center Plaza, 
Philadelphia, PA 19102 (US). 



(81) Designated States: AL, AM, AT, AU, AZ, BA, BB, BG, BR, 
BY, CA, CH, CN, CU, CZ, DE, DK, EE, ES, FI, GB, GD, 
GE, GH, GM, HR, HU, ID, IL, IN, IS, JP, RE, KG, KP, 
KR, KZ, LC, LK, LR, LS, LT, LU, LV, MD, MG, MK, 
MN, MW, MX, NO, NZ, PL, PT, RO, RU, SD, SE, SG, 
SI, SK, SL, TJ, TM, TR, TT, UA, UG, UZ, VN, YU, ZW, 
ARIPO patent (GH, GM, KE, LS, MW, SD, SZ, UG, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, CH, CY, DE, DK, ES, FI, FR, 
GB, GR, IE, IT, LU, MC, NL, PT, SE), OAPI patent (BF, 
BJ, CF, CG, CI, CM, GA, GN, GW, ML, MR, NE, SN, 
TD, TG). 



Published 

With international search report 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 



(54) Title: A METHOD AND SYSTEM FOR PATTERN RECOGNITION AND PROCESSING 



JO 



,12 



Data 



Input Layer 



.14 



Association 
Layer 



4-> 



Memory 



(57) Abstract 

The present invention 
provides a method and system 
(1) for pattern recognition and 
processing. Information represen- 
tative of physical characteristics 
or representations of physical 
characteristics is transformed into 
a Fourier series in Fourier space 
(12) within an input context of 
the physical characteristics that 
is encoded in time as delays 
corresponding to modulation of 
the Fourier series at corresponding 
frequencies. Associations (14) 
are formed between Fourier series 
by filtering the Fourier series 
and by using a spectral similarity 
between the filtered Fourier series 
to determine the association based 
on Poissonian probability. The 
associated Fourier series are added 
to form strings (16) of Fourier 
series. Each string is ordered by 
filtering it with multiple selected 
filters to form multiple time order 
formatted subset Fourier series, and 
by establishing the order through 
associations with one or more 
initially ordered strings to form an 
ordered string. Associations are 

formed between the ordered strings to form complex ordered strings (18) that relate similar items of interest. 
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A METHOD AND SYSTEM FOR PATTERN RECOGNITION AND 

PROCESSING 

Cross Reference to Related Applications 
5 ~ 

This application claims the benefit of United States provisional 

application Serial No. 60/068,834, filed December 24, 1997. 

Backgro und Of the Invention 

1 0 

Attempts have been made to create pattern recognition systems 
using programming and hardware. The state of the art includes neural 
nets. Neural nets typically comprise three layers-an input layer, a 
hidden layer, and an output layer. The hidden layer comprises a series of 

1 5 nodes which serve to perform a weighted sum of the input to form the 

output. Output for a given input is compared to the desired output, and a 
back projection of the errors is carried out on the hidden layer by 
changing the weighting factors at each node, and the process is reiterated 
until a tolerable result is obtained. The strategy of neural nets is 

2 0 analogous to the sum of least squares algorithms. These algorithms are 

adaptive to provide reasonable output to variations in input, but they can 
not create totally unanticipated useful output or discover associations 
between multiple inputs and outputs. Their usefulness to create novel 
conceptual content is limited; thus, advances in pattern recognition 

2 5 systems using neural nets is limited. 

SUMMARY OF THE INVENTION 

The present invention is directed to a method and system for 

3 0 pattern recognition and processing involving processing information in 

Fourier space. 

The system of the present invention includes an Input Layer for 
receiving data representative of physical characteristics or 
representations of physical characteristics capable of transforming the 
3 5 data into a Fourier series in Fourier space. The data is received within an 
input context representative of the physical characteristics that is 
encoded in time as delays corresponding to modulation of the Fourier 
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series at corresponding frequencies. The system includes a memory that 
maintains a set of initial ordered Fourier series. The system also includes 
an Association Layer that receives a plurality of the Fourier series in 
Fourier space including at least one ordered Fourier series from the 
5 memory and forms a string comprising a sum of the Fourier series and 
stores the string in memory. The system also includes a String Ordering 
Layer that receives the string from memory and orders the Fourier series 
contained in the string to form an ordered string and stores the ordered 
string in memory. The system also includes a Predominant Configuration 
1 0 Layer that receives multiple ordered strings from the memory, forms 
complex ordered strings comprising associations between the ordered 
strings, and stores the complex ordered strings to the memory. The 
components of the system are active based on probability using 
weighting factors based on an activation rates. 

1 5 One aspect of the present invention is directed to inputting 

information as data to the system within an input context and associating 
the data. This aspect of the invention includes encoding the data as 
parameters of at least two Fourier components in Fourier space, adding 
the Fourier components to form at least two Fourier series in Fourier 

2 0 space, the Fourier series representing the information, sampling at least 

one of the Fourier series in Fourier space with a filter to form a sampled 
Fourier series, and modulating the sampled Fourier series in Fourier 
space with the filter to form a modulated Fourier series. This aspect of 
the invention also includes determining a spectral similarity between the 

2 5 modulated Fourier series and another Fourier series, determining a 

probability expectation value based on the spectral similarity, and 
generating a probability operand having a value selected from a set of 
zero and one, based on the probability expectation value. These steps are 
repeated until the probability operand has a value of one. Once the 

3 0 probability operand has a value of one, the modulated Fourier series and 

the other Fourier series are added to form a string of Fourier series in 
Fourier space, and the string of Fourier series is stored in the memory. 

Another aspect of the present invention is directed to ordering a 
string representing the information. This aspect of the invention utilizes 
3 5 a High Level Memory section of the memory that maintains an initial set 
of ordered Fourier series. This aspect of the invention includes obtaining 
a string from the memory and selecting at least two filters from a 
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selected set of filters stored in the memory. This aspect also includes 
sampling the string with the filters such that each of the filters produce a 
sampled Fourier series. Each Fourier series comprises a subset of the 
string. This aspect also includes modulating each of the sampled Fourier 
5 series in Fourier space with the corresponding selected filter such that 
each of the filters produce an order formatted Fourier series. 
Furthermore, this aspect includes adding the order formatted Fourier 
series produced by each filter to form a summed Fourier series in Fourier 
space, obtaining an ordered Fourier series from the memory, determining 
1 0 a spectral similarity between the summed Fourier series and the ordered 
Fourier series, determining a probability expectation value based on the 
spectral similarity, and generating a probability operand having a value 
selected from a set of zero and one, based on the probability expectation 
value. These steps are repeated until the probability operand has a value 

1 5 of one. Once the probability operand has a value of one, this aspect 

includes storing the summed Fourier series to an intermediate memory 
section. Thereafter, this aspect includes removing the selected filters 
from the selected set of filters to form an updated set of filters, removing 
the subsets from the string to obtain an updated string, and selecting an 

2 0 updated filter from the updated set of filters. This aspect further 

includes sampling the updated string with the updated filter to produce a 
sampled Fourier series comprising a subset of the string, modulating the 
sampled Fourier series in Fourier space with the corresponding selected 
updated filter to produce an updated order formatted Fourier series, 

2 5 recalling the summed Fourier series from the intermediate memory 

section, adding the updated order formatted Fourier series to the 
summed Fourier series to form an updated summed Fourier series in 
Fourier space, and obtaining an updated ordered Fourier series from the 
memory. This aspect further includes determining a spectral similarity 

3 0 between the updated summed Fourier series and the updated ordered 

Fourier series, determining a probability expectation value based on the 
spectral similarity, and generating a probability operand having a value 
selected from a set of zero and one, based on the probability expectation 
value. These steps are repeated until the probability operand has a value 
3 5 of one or all of the updated filters have been selected from the updated 
set of filters. If all of the updated filters have been selected before the 
probability operand has a value of one, then clearing the intermediate 
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memory section and repeating the steps starting with selecting at least 
two filters from a selected set of filters. Once the probability operand has 
a value of one, the updated summed Fourier series is stored to the 
intermediate memory section and steps beginning with removing the 
5 selected filters from the selected set of filters to form an updated set of 
filters are repeated until one of the following set of conditions is satisfied: 
the updated set of filters is empty or the remaining subsets of the string 
is nil. If the remaining subset of strings is nil, then the Fourier series in 
the intermediate memory section is stored in the High Level Memory 
1 0 section of the memory. 

Another aspect of the present invention is directed to forming 
complex ordered strings by forming associations between a plurality of 
ordered strings. This aspect of the invention includes recording ordered 
strings to the High Level Memory section, forming associations of the 

1 5 ordered strings to form complex ordered strings, and recording the 

complex ordered strings to the High Level Memory section. A further 
aspect of the invention is directed to forming a predominant 
configuration based on probability. This aspect of the invention includes 
generating an activation probability parameter, storing the activation 

2 0 probability parameter in the memory, generating an activation 

probability operand having a value selected from a set of zero and one, 
based on the activation probability parameter, activating any one or 
more components of the present invention such as matrices representing 
functions, data parameters, Fourier components, Fourier series, strings, 

2 5 ordered strings, components of the Input Layer, components of the 

Association Layer, components of the String Ordering Layer, and 
components of the Predominant Configuration Layer, the activation of 
each component being based on the corresponding activation probability 
parameter, and weighting each activation probability parameter based on 

3 0 an activation rate of each component. 

Brief Description of the Drawings 

FIGURE 1 is a high level block diagram illustrating an embodiment of 
3 5 the present invention; 

FIGURE 2 is a detailed block diagram illustrating an Input Layer, an 
Association Layer, and a memory layer of the embodiment of FIGURE 1; 
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FIGURE 3 is a flow diagram of an exemplary transducer data structure 
of a time delay interval subdivision hierarchy wherein the data from a 
transducer having n levels of subcomponents numbering integer m per 
level is assigned a master time interval with n + l sub time intervals in a 
5 hierarchical manner wherein the data stream from the final n th level 
transducer element is recorded as a function of time in the n + l th time 
coded sub memory buffer in accordance with the present invention; 

FIGURE 4 is a detailed block diagram illustrating an String Ordering 
Layer and the memory layer of the embodiment of FIGURE 1; 
1 0 FIGURE 5 is a detailed block diagram illustrating a Predominant 

Configuration Layer and the memory of the embodiment of FIGURE 1 in 
relation to the Input Layer, the Association Layer, and the String 
Ordering Layer of the present invention; 

FIGURE 6 is a schematic drawing of the "P or M element response" 

1 5 comprised of a series of seven "impulse responses" in accordance with the 

invention; 

FIGURE 7 is a schematic drawing of the Fourier Transform H[k p ,k z ] of 
the system function h(p,z) corresponding to the "impulse response" in 
accordance with the invention; 

2 0 FIGURE 8 is a schematic of h(t) given by Eq. (37.51) where a = l and 

TV = 100 in accordance with the invention; 

FIGURE 9 is a schematic of H(f) given by Eq. (37.50) where a=l and 
Af = 100 in accordance with the invention; 

FIGURE 10 is a schematic of h(t) given by Eq. (37.51) where a = 10 and 

2 5 TV = 100 in accordance with the invention; 

FIGURE 11 is a schematic of H(f) given by Eq. (37.50) where a = 10 

and Af = 100 in accordance with the invention; 

FIGURE 12 is a schematic of h{t) given by Eq. (37.51) where a = l and 
AT = 500 in accordance with the invention; 

3 0 FIGURE 13 is a schematic of H(f) given by Eq. (37.50) where cc = l and 

iV = 500 in accordance with the invention; 

FIGURE 14 is a schematic of h(t) given by Eq. (37.51) where a = 10 and 
TV = 500 in accordance with the invention; 

FIGURE 15 is a schematic of H(f) given by Eq. (37.50) where <* = 10 
3 5 and TV = 500 in accordance with the invention; 

FIGURES 16A, 16B, and 16C illustrate plots of the probability P A ((j>) 

(Eq. (37.106a)) of association of the corresponding Fourier series based on 
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a first active association ensemble coupling with a second association 
ensemble as a function of frequency difference angle, <p s , coupling cross 
section amplitude, & 2 , and phase shift, 5 S =0 wherein the parameter 
& 2 =0.01, 0*25, and 1.00, respectively, in accordance with the invention; 

5 FIGURES 17A, 17B, 17C, and 17D illustrate plots of the probability 

P A ((p) (Eq. (37.106a)) of association of the corresponding Fourier series 

based on a first active association ensemble coupling with a second 
association ensemble as a function of frequency difference angle, <f> s , and 
phase shift, S s , for the coupling cross section amplitude, /? s 2 =0.25, wherein 
10 the parameter <5 5 =0, 0.25 /r, 0.50 ;r, and respectively, in accordance with 
the present invention; 

FIGURE 18 is a flow diagram of an exemplary hierarchical relationship 
between the characteristics and the processing and storage elements in 
accordance with the present invention; 

1 5 FIGURE 19 is a flow diagram of an exemplary hierarchical relationship 

of the signals in Fourier space comprising FCs, SFCs, groups of SFCs, and a 
string in accordance with the present invention; 

FIGURE 20 is an exemplary layer structure in accordance with the 
present invention, and 

2 0 FIGURE 21 is a flow diagram of an exemplary layer structure and 

exemplary signal format in accordance with the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

2 5 The present invention is directed to systems and methods for 

performing pattern recognition and association based upon receiving, 
storing, and processing information. The information is based upon 
physical characteristics or representations of physical characteristics and 
a relationship of the physical characteristics, hereinafter referred to as 

3 0 physical context, of an item of interest. The physical characteristics and 

physical context serve as a basis for stimulating a transducer. The 
transducer converts an input signal representative of the physical 
characteristics and the physical context into the information for 
processing. The information is data and an input context. The data is 
3 5 representative of the physical characteristics or the representations of 

physical characteristics and the input context corresponds to the physical 
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context based upon the identity of a specific transducer and its particular 
transducer elements. The physical context maps on a one to one basis to 
the input context. The information defines a Fourier series in Fourier 
space that represents the item of interest. In other words, a Fourier 
5 series in Fourier space represents the information parameterized 
according to the data and the input context. In addition, the input 
context maps on a one to one basis to an Input Layer section of a 
memory. Thus, there is a one to one map of physical context to input 
context to Input Layer section of a memory. The representation of 
1 0 information as a Fourier series in Fourier space allows for the mapping. 

As illustrated in Figure 1, at a high level, the system 10 includes 
several function specific layers. These include an Input Layer 12, an 
Association Layer 14, an String Ordering Layer 16 and a Predominant 
Configuration Layer 18. The Input Layer 12 receives the data within the 

1 5 input context and transforms the data into the Fourier series in Fourier 

space representative of the information. The system 10 also includes a 
memory 20 for storing information. The Input Layer 12 also encodes the 
input context as delays in time corresponding to a modulation factor of 
the Fourier series at corresponding frequencies. The Association Layer 

2 0 14 receives a plurality of Fourier series in Fourier space, including at 

least one ordered Fourier series from the memory 20, forms a string 
comprising a sum of the Fourier series and stores the string to the 
memory 20. The String Ordering Layer 16 receives the string from the 
memory 20, orders the Fourier series contained in the string to form an 

2 5 ordered string and stores the ordered string in the memory 20. The 

Predominant Configuration Layer 18 receives multiple ordered strings 
from the memory 20, forms associations between the ordered strings to 
form a complex ordered string, also referred to as a predominant 
configuration string, and stores the predominant configuration string to 

3 0 the memory 20. The memory 20 may be partitioned in several distinct 

sections for storing different types of information or distinctly classified 
types of information. Specifically, the memory may include a High Level 
Memory section, an intermediate level memory section, etc. as will be 
described in more detail below. 
3 5 The following example illustrates how the present invention 

processes the physical characteristics of an item of interest, specifically a 
triangle. In flat geometry, the physical characteristics of a triangle are 
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three connected lines at angles aggregating to 180°. The physical 
characteristics provide spatial variations of light scattering. In one 
embodiment, a light responsive transducer (not shown) of the system 10 
transduces the light scattering into the data. An exemplary transducer is 
5 a charge coupled device ("CCD") array. One data element at a point in 

time may be a voltage of a particular CCD element of the CCD array. Each 
CCD element of the CCD array has a spatial identity. The physical context 
for the triangle is the relationship of the lines at the corresponding angles 
providing a spatial variation of light scattering. The input context is the 

1 0 identity of each CCD element that responds according to the physical 

context. For example, a CCD element (100,13) of a 512 by 512 CCD array 
will uniquely respond to light scattered by the lines and angular relations 
of the triangle relative to the other CCD elements of the CCD array. The 
response is stored in a specific memory register of an Input Layer section 
15 of the memory 20. The specific memory register is reflective of the input 
context. In the present invention, a Fourier series in Fourier space 
represents the information of the triangle parameterized according to the 
voltage and the CCD element identity. 

Referring to FIGURE 2, in the first step, the Input Layer 12 receives 

2 0 the data from the transducer (not shown). A Fourier transform processor 

22 encodes each data element as parameters of a Fourier component in 
Fourier space and stores the data parameter values to the Input Layer 
section 24 of the memory 20. Each Fourier component of the Fourier 
series may comprise a quantized amplitude, frequency, and phase angle. 
2 5 For example the Fourier series in Fourier space may be: 

4k 

T 2 



1*-- 1 + 

K p 



f 2k 
PoJ 



sin 



2k 

k z — n 

z o« J 



having a quantized amplitude, frequency, and phase angle, wherein a Qm is 
a constant, k p and k z are the frequency variables, n 9 m, and M are 
integers, and A^, N m ^ 9 p 0w , and z Qtn are the data parameters. 
30 In a first embodiment, the data parameters N and N of the 

Fourier series component are proportional to the rate of change of the 
physical characteristic. Each of the data parameters p 0m and z 0m of each 

Fourier component is inversely proportional to the amplitude of the 
physical characteristic. In the triangle example, the amplitude of the 
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voltage at a given CCD element relative to the neighboring CCD element 
defines the rate of change of the voltage which is converted into the data 
parameters and N^. The inverse of the amplitude of the voltage of 

each CCD element is converted into the data parameters p 0 and z 0h . As 
5 illustrated in FIGURE 3 and described above, for each CCD element, the 
Fourier series, parameterized accordingly, are stored to a specific sub 
register 27 of a specific register 26 of the Input Layer section 24 of the 
memory 20. Since the structure of a Fourier series is known in the art, 
only the parameters need to be stored in a digital embodiment. 

1 0 The number and types of transducers that may supply information 

to the system is only limited by available technology, hardware and 
economics, as is the number m of corresponding registers 26 for each 
transducer. Each register 26 may have any number d of subregisters 27, 
where the number d of subregisters of one register 26 is not necessarily 
15 the same as other registers .26. For example, "Level 1 " register "1 " may 
have thirty "Level 2 " subregisters 27 and "Level 1 " register "2 " may 
have one-hundred subregisters 27. Furthermore, each "Level 2 " register 
may have any number e of subregisters, where the number e of 
subregisters of one register 27 is not necessarily the same as other 

2 0 registers 27. For example, "Level 2 " register "1 " may have fifty "Level 

n" subregisters 29 and "Level 2 " register "2 " may have seventy "Level 
«" subregisters 29. Still further, each "Level n" register 29 may have any 
number / of time buffer elements 31, where the number / of time 
buffer elements 31 is not necessarily the same as other time buffer 

2 5 elements 31. 

In a second embodiment, each of the data parameters N m ^ and N nu 

of the Fourier series component is proportional to the amplitude of the 
physical characteristic. Each of the data parameters p 0 ^ and ^ of each 

Fourier component is inversely proportional to the rate of change of the 

3 0 physical characteristic. As in the first embodiment, for each CCD element, 

these parameters are stored in a specific sub register of the Input Layer 
section of the memory. 

In a third embodiment, each of the data parameters N and N of 

the Fourier series component is proportional to the duration of the signal 
3 5 response of each transducer. Each of the data parameters p K and z Qm of 

each Fourier component is inversely proportional to the physical 
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1 0 

characteristic. As in the first embodiment, for each CCD element, these 
parameters are stored in a specific sub register of the Input Layer section 
of the memory. 

As an alternative example, the Fourier series in Fourier space may 

be: 



kl 



k -n 



sin 



2k 
■ n 

J 



having a quantized frequency, and phase angle, wherein a 0m is a constant, 
k p and k z are the frequency variables, «, m, and M are integers, and N , 
N mzo > Po m > and z 0m are the data parameters. As described with respect to 

1 0 the previous example, for each CCD element, these parameters are stored 

in a specific sub register of the Input Layer section of the memory. 

The physical context is conserved by mapping with a one to one 
basis between the physical context and the input context based on the 
identity of each transducer. The input context is conserved by mapping 
15 on a one to one basis to the Input Layer section 24 of memory 20. In an 
embodiment, the input context is encoded in time as a characteristic 
modulation frequency band in Fourier space of the Fourier series. The 
characteristic modulation frequency band in Fourier space represents the 
input context according to the identity of a specific transducer of the 

2 0 relationship of two transducer elements. The modulation within each 

frequency band may encode not only input context but context in a 
general sense. The general context may encode temporal order, cause 
and effect relationships, size order, intensity order, before-after order, 
top-bottom order, left-right order, etc. all of which are relative to the 

2 5 transducer. 

Still referring to FIGURE 3, the transducer has n levels of 
subcomponents. Each transducer is assigned a portion 26 of the Input 
Layer section 24 of the memory 20. The memory 20 is arranged in a 
hierarchical manner. Specifically, the memory is divided and assigned to 

3 0 correspond to a master time interval with « + l sub time intervals. The 

hierarchy parallels the n levels of the transducer subcomponents. The n 
th level transducer sub component provides a data stream to the system 
10. The data stream is recorded as a function of time in the n + \ sub time 
interval. The time intervals represent time delays which correspond to 
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1 1 

the characteristic modulation frequency band in Fourier space which in 
turn represents the input context according to the specific transducer or 
transducer subcomponent. 

An exemplary complex transducer which may be represented by a 
5 data structure comprising a hierarchical set of time delay intervals is a 
CCD array of a video camera comprising a multitude of charge coupled 
devices (CCDs). Each CCD comprises a transducer element and is 
responsive to light intensity of a given wavelength band at a given 
spatial location in a grid. Another example of a transducer is an audio 

1 0 recorder comprising transducer elements each responsive to sound 

intensity of a given frequency band at a given spatial location or 
orientation. A signal within the band 300-400 MHz may encode and 
identify the signal as a video signal; whereas, a signal within the band 
500-600 MHz may encode and identify the signal as an audio signal. 
15 Furthermore, a video signal within the band 315-325 MHz may encode 
and identify the signal as a video signal as a function of time of CCD 
element (100,13) of a 512 by 512 array of CCDs. 

In one embodiment, the characteristic modulation having a 
frequency within the band in Fourier space is represented by e ~ nm . The 

2 0 modulation corresponds to the time delay S(t-t 0 ) wherein / is the 

frequency variable, t is the time variable, and t Q is the time delay. The 
characteristic modulation is encoded as a delay in time by storing the 
Fourier series in a specific portion of the Input Layer section of the 
memory wherein the specific portion has n + 1 sub time intervals. Each 
2 5 sub time interval corresponds to a frequency band. 

In an alternative embodiment, the characteristic modulation, 
having a frequency within the band is represented by e ~ Jk ^ p ^ +p) -\ Thus, 
the Fourier series in Fourier space may be: 



X 1 ^«oN mm N m y k >^^ s in 

m=l n=-oo 1 | K z \ 2 2 j 

kl 



sin 



p 2 2 



3 0 wherein p, m = v, m t, m is the modulation factor which corresponds to the 
physical time delay t p =v Jb t Jb is the modulation factor which 

corresponds to the specific transducer time delay , v t and v A are 
constants such as the signal propagation velocities, a 0 is a constant, k 
and k z are the frequency variables, «, m, and M are integers, and N m , 
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N ™« ' Po„> and z 0m are data parameters. The data parameters are selected 
in the same manner as described above. 

Transducer strings may be created by obtaining a Fourier series 
from at least two selected transducers and adding the Fourier series. 
5 Transducers that are active simultaneously may be selected. The 
transducer string may be stored in a distinct memory location of the 
memory. The characteristic modulation, having a frequency within the 
band in Fourier space can be represented by e~ n ^ a which corresponds to 
the time delay S(t-t 0 ) wherein / is the frequency variable, t is the time 
1 0 variable, and t Q is the time delay. 

Recalling any part of the transducer string from the distinct 
memory location may thereby cause additional Fourier series of the 
transducer string to be recalled. In other words the Fourier series are 
linked. Fourier series, in addition to those of transducer strings may be 

1 5 linked. In order to achieve linking of the Fourier series, the system 

generates a probability expectation value that recalling any part of one of 
the Fourier series from the memory causes at least another Fourier series 
to be recalled from the memory. The system stores the probability 
expectation value to memory. The system generates a probability 

2 0 operand having a value selected from a set of zero and one, based on the 

probability expectation value. The system recalls at least another Fourier 
series from the memory if the operand is one. The probability 
expectation value may increase with a rate of recalling any part of any of 
the Fourier series. 

25 The system may be initialized by learning. The relationship 

between the data and the data parameters such as p n and N of each 

component of the Fourier series is learned by the system by applying 
standard physical signals. In the case of the triangle example, the 
standard physical signals are the scattered light from the physical 

*0 characteristics of the triangle. The physical signals are applied to each 
transducer together with other information that is associated with the 
standard. A data base is established. This information that is associated 
with the standard is recalled and comprises input into the Association 
Layer and the String Ordering Layer. 

5 5 The data parameters and the input context are established and 

stored in the Input Layer section 24 of the memory 20. 
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Referring again to Figure 2, several the parameterized Fourier 
components are input to the Association Layer to form associations of the 
Fourier series. The Fourier components may be stored in a Fourier 
component section 30 of a temporary memory section 28. The Fourier 
5 components are added to form multiple Fourier series which in turn may 
be stored in a Fourier series section 32 of the temporary memory section 
28. At least one of the Fourier series stored in the Fourier series section 
32 is input to a filter 34 wherein the filter 34 samples and modulates the 
Fourier series. The filtered Fourier series is input to a spectral similarity 
1 0 analyzer 36. The spectral similarity analyzer 36 determines the spectral 
similarity between the filtered Fourier series and another Fourier series 
stored in the Fourier series section 32 of the temporary memory section 
28. A spectral similarity value is output from the spectral similarity 
analyzer 36 and input to a probability expectation analyzer 38. The 

1 5 probability expectation analyzer 38 determines a probability expectation 

value based on the spectral similarity value. The probability expectation 
value output from the probability expectation analyzer 38 is input to a 
probability operand generator 40. The probability operand generator 40 
generates a probability operand value of one or zero based upon the 

2 0 probability expectation value. The probability operand value is output to 

a processor 42. If the probability operand value is zero, the processor 42 
sends another Fourier series from the Fourier series section 32 of the 
temporary memory section 28 to the filter 34 and begins the process 
again. If the probability operand value is one, the filtered Fourier series 
2 5 and the other Fourier series are added to form a string and the string is 
stored in a string memory section 44. 

The filter 34 can be a time delayed Gaussian filter in the time 
domain. The filter may be characterized in time by: 



3 0 wherein J — is a delay parameter, a is a half-width parameter, and / is 




2 



a 



V2tt 



e 




the time parameter. The Gaussian filter may comprise a plurality of 
cascaded stages each stage having a decaying exponential system 
function between stages. The filter, in frequency space, can be 
characterized by: 
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wherein and a are a corresponding delay parameter and a half- 

width parameter in time, respectively, and / is the frequency parameter. 
The probability distribution may be Poissonian. Thus, the probability 
expectation value can be based upon Poissonian probability, 
probability expectation value may be characterized by 



The 



(P - Pi J^; 2 (^^)]cos(5 J+ 2sin0 f ) 



wherein P is the maximum probability of at least one other Fourier 
series being associated with a first Fourier series, p Ts is a probability of at 

least one other Fourier series being associated with a first Fourier series 
in the absence of coupling of the first Fourier series with the at least one 
other Fourier series, p) is a number that represents the amplitude of 
spectral similarity between at least two filtered or unfiltered Fourier 
series, <p s represents the frequency difference angle between at least two 
filtered or unfiltered Fourier series, and 8 S , is a phase factor, p) may be 
characterized by 



# = (8*) a 



1 



a\a] 



2k v at + a: 



A/, M s 

X Oo^*, X Oo^exp- 



m s = l 



at 



at 



a, a„ 2v m 2v_. 



and - correspond to delay parameters of a first and s-th time 



delayed Gaussian filter, respectively, a x and a s corresponding half-width 

parameters of a first and s-th time delayed Gaussian filter, respectively, 
M, and M s are integers, and a 0 ^ are constants, v mx and v nij are 

constants such as the signal propagation velocities, and N m , N m , p ft , and 

P 0h)s are data parameters. The data parameters are selected in the same 

manner as described above. <p s may be characterized by 
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K 



a, ~i 2v 



a, 



and 



correspond to delay parameters of a first and s-th time 



delayed Gaussian filter, respectively, a, and a s corresponding half-width 

parameters of a first and s-th time delayed Gaussian filter, respectively, 
5 M x and M s are integers, a Q ^ and are constants, v mi and v mj are 

constants such as the signal propagation velocities, and N ni 9 N , p n , and 

Po ms are data parameters. The data parameters are selected in the same 

manner as described above. 

An exemplary string with a characteristic modulation having a 

10 frequency within the band represented by e ~ jk ^ pjb - +p '^ [ s: 



S M s ™ 

XXX 



sin 



k p -n 



Ik 



Po 



N. 



„A>,„ 



sin 



k z —n 



2k 



N. 



= v ft / ft is the modulation factor which 



and Vr. are 

J u s.m 



wherein p t = v, t t is the modulation factor which corresponds to the 

's,m *s,m 's r m ± 

physical time delay t t 

corresponds to the specific transducer time delay f A , v, 
constants such as the signal propagation velocities, a 0 ^ is a constant, k p 
and k z are the frequency variables, n y m, s, Af„ and 5 are integers, and 
N s. mpo > N s, mzo > Po, M > and Zo Mjm are data parameters. The data parameters are 
selected in the same manner as described above. 

An exemplary string with each Fourier series multiplied by the 
Fourier transform of the delayed Gaussian filter represented by 



e v e , „ e .v € !«.«, that established the association to 

form the string is: 
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XX X -^rr** N. 

s=l m~l n — 



1+ l 



e -*>("A- + "-) sin 



k p -n 



2k 



Po 



N. 



„A>,„ 



J,Hi J 



(f 



sin 



2k 



si. are delay 



wherein v jpo and v^ 0 are constants such as the signal propagation 
velocities in the p and z directions, respectively, \ NspQ and 
parameters and a spo and a sz0 are half-width parameters of a 
5 corresponding Gaussian filter in the p and z directions, respectively, 
Pi sm = v t s Ji 5yfn is the modulation factor which corresponds to the physical 
time delay * w Pjb sm ==v jb Jjb is the modulation factor which corresponds 



to the specific transducer time delay t 



and v 



are constants such 



as the signal propagation velocities, a 0 is a constant, k n and k are the 
frequency variables, «, m, j, M„ and 5 are integers, and N tm , M , p n , 
and z 0 m are data parameters. The data parameters are selected in the 
same manner as described above. 

Therein, the Association Layer forms associations between Fourier 
series and sums the associated Fourier series to form a string. The string 
is then stored in the string memory section. 

The next aspect of the present invention is the ordering of the 
strings stored in the string memory section 44. The ordering may be 
according to any one of the following: temporal order, cause and effect 
relationships, size order, intensity order, before-after order, top-bottom 
2 0 order, or left-right order. Referring to FIGURE 4, the method for ordering 
the strings stored in the string memory section 44 entails the following: 

a. ) obtaining a string from the string memory section 44 and storing 
the string to a temporary string memory section 46; 

b. ) selecting at least two filters 48, 50 from a selected set of filters 

52; 

c. ) sampling the string with the filters 48, 50, each of the filters 
forming a sampled Fourier series, each Fourier series comprising a subset 
of the string; 
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d. ) modulating each of the sampled Fourier series in Fourier space 
with the corresponding selected filter 48, 50, each forming an order 
formatted Fourier series; 

e. ) adding the order formatted Fourier series to form a summed 
5 Fourier series in Fourier space; 

f. ) obtaining an ordered Fourier series from the High Level Memory 
section 54; 

g. ) determining a spectral similarity with a spectral similarity 
analyzer 56 between the summed Fourier series and the ordered Fourier 

1 0 series; 

h. ) determining a probability expectation value, with a probability 
expectation value analyzer 58 based on the spectral similarity; 

L) generating a probability operand, with a probability operand 
generator 60 having a value selected from a set of zero and one, based on 

1 5 the probability expectation value; 

j.) repeating steps b-i until the probability operand has a value of 
one as determined by the processor 42; 

k.) storing the summed Fourier series to an intermediate memory 
section 62; 

2 0 1.) removing the selected filters from the selected set of filters 52 to 

form an updated set of filters 52; 

m.) removing the subsets from the string to obtain an updated 
string; 

n.) selecting an updated filter 64 from the updated set of filters; 

2 5 o.) sampling the updated string with the updated filter to form a 

sampled Fourier series comprising a subset of the string; 

p.) modulating the sampled Fourier series in Fourier space with the 
corresponding selected updated filter to form an updated order 
formatted Fourier series; 

3 0 q.) recalling the summed Fourier series from the intermediate 

memory section 62; 

r.) adding the updated order formatted Fourier series to the 
summed Fourier series from the intermediate memory section to form an 
updated summed Fourier series in Fourier space; 
3 5 s.) obtaining another ordered Fourier series from the High Level 

Memory section 54; 
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t.) determining a spectral similarity between the updated summed 
Fourier series and the another ordered Fourier series; 

u.) determining a probability expectation value based on the 
spectral similarity; 

v.) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value; 

w.) repeating steps n-v until the probability operand has a value of 
one or all of the updated filters have been selected from the updated set 
of filters as determined by processor 42; 

x.) if all of the updated filters have been selected before the 
probability operand has a value of one, then clearing the intermediate 
memory section and returning to step b; 

y.) if the probability operand has a value of one, then clearing the 
intermediate memory section and storing the updated summed Fourier 
series to the intermediate memory section; 

z.) repeating steps 1-y until the one of the following set of 
conditions is satisfied: the updated set of filters is empty, or the 
remaining subsets of the string of step m.) is nil as determined by the 
processor 42; 

aa.) storing the Fourier series of intermediate memory section to 
the High Level Memory section 54. 

Each filter of the set of filters can be a time delayed Gaussian filter 
having a half-width parameter a which determines the amount of the 
string that is sampled. Each filter of the set of filters can be a time 

delayed Gaussian filter having a delay parameter J— which corresponds 



to a time point. Each Fourier series of the ordered string can be 
multiplied by the Fourier transform of the delayed Gaussian filter 



established the correct order. The ordered string can be represented by: 




represented by 
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N._ N. 



-j(-i)'-'J5<-"»-i(-t;)'-e'-»''» 



1+- 
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1 5 



20 



k p -n 



2k 



Po 



t.m J 



(f 



sin 



2k 



wherein 



V s,m t 0 s ^ m ) 



v s P0 and V ^ 0 are constants such as the signal propagation 



velocities in the p and z directions, respectively, ^ei and JSZ are delay 
parameters and a spo and a sz0 are half-width parameters of a 
corresponding Gaussian filter in the p and z directions, respectively, 
p <s.m = \ m r <s, m is the mod ulation factor which corresponds to the physical 
time delay r^, = v^j^ is the modulation factor which corresponds 
to the specific transducer time delay , i 



and 



are constants such 



as the signal propagation velocities, a 0gm is a constant, k p and k z are the 
frequency variables, n, m, ^, M,, and 5 are integers, and N , N , p n 
and z o„„ are data parameters. The data parameters are selected in the 
same manner as described above. 

The probability expectation value may be based upon Poissonian 
probability. The probability expectation value is represented by 



n 



Pr. + ( P ~ P,shp[-K 2 { 1 C ° S2 ^ )]cos(^ + 2sin ft) 



wherein P is the maximum probability that at least one other Fourier 
series is active given that a first Fourier series is active, p r is a 

probability of a Fourier series becoming active in the absence of coupling 
from at least one other active Fourier series, f} 2 s is a number that 
represents the amplitude of spectral similarity between at least two 
filtered or unfiltered Fourier series, <p s represents the frequency 

difference angle between at least two filtered or unfiltered Fourier series, 
and 8 S , is a phase factor. fi] may be represented by 
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wherein p = v t and p, = v, t, are the modulation factors which 
corresponds to the physical time delays t, m and t tm , respectively, 
Pjb nh =v /b m} t /b mi and P/b m ^ v fb m t fb m are the modulation factors which 
corresponds to the specific transducer time delay t, and , 

■f t? m l J"m s 

respectively, v, , v, , v. , and v. are constants such as the signal 



propagation velocities, 



and 



correspond to delay parameters of a 



first and s-th time delayed Gaussian filter, respectively, a, and 

corresponding half-width parameters of a first and s-th time delayed 
Gaussian filter, respectively, M x and M s are integers, a^, a 0m are 

constants, v Mi and v mj are constants such as the signal propagation 
velocities, and A^, N m ^ , p 0 ^ , and p 0 are data parameters. The data 
parameters are selected in the same manner as described above. <p s may 
be represented by 



K 



a 



2v„ 



"A. 



v. 



'mi ) 



P'.„ 



2v„ 



"A. 



a. 



2v *. V A„ 



-+- 



•mi y 



wherein p, = v, r, and p f = v, r are the modulation factors which 
corresponds to the physical time delays t, m and t t , respectively, 

^ = XVm, and ~ v f b m t fi>m are ^ e m °d ulat i° n factors which 
corresponds to the specific transducer time delay r A and r A , 
respectively, v, m , v, m , , and are constants such as the signal 



propagation velocities, 



and 



correspond to delay parameters of a 



first and s-th time delayed Gaussian filter, respectively, a x and a s 
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corresponding half-width parameters of a first and s-th time delayed 
Gaussian filter, respectively, Af,, and M s are integers, a 0 and a 0 are 

constants, v ffl[ and v mj are constants such as the signal propagation 
velocities, and N mi , N m f , p 0 bi , and p 0 ^ are data parameters. The data 

5 parameters are selected in the same manner as described above. 

The String Ordering Layer produces an ordered string of Fourier 
series, wherein the ordered string is stored in the High Level Memory 
section. 

The next aspect of the present invention is the formation of a 
1 0 predominant configuration by forming complex ordered strings through 
the association of ordered strings. Referring to FIGURE 5, the method for 
forming the complex ordered strings from strings stored in the string 
memory section entails the following. The Predominant Configuration 
Layer 18 receives ordered strings from the High Level Memory section 

1 5 54 and forms more complex ordered strings by forming associations 

between the ordered strings. The complex ordered strings are stored in 
the complex ordered string section 72 of the memory 20. 

The Predominant Configuration Layer 18 also activates components 
within the Input Layer 12, the Association Layer 14, and the String 

2 0 Ordering Layer 16. The layers of the present invention may be treated 

and implemented as abstract data types in the art of computer science 
relating to object-oriented programming. The components of the layers 
therefore refer to all classes, instances, methods, attributes, behaviors, 
and messages of the layer abstractions as defined above. A class is the 

2 5 implementation of an abstract data type (ADT). It defines attributes and 

methods implementing the data structure and operations of the ADT, 
respectively. Instances of classes are called objects. Consequently, classes 
define properties and behavior of sets of objects. An object can be 
uniquely identified by its name and it defines a state which is 

3 0 represented by the values of its attributes at a particular time. The 

behavior of an object is defined by the set of methods which can be 
applied to it. A method is associated with a class. An object invokes a 
method as a reaction to receipt of a message. 

Thus, the components of a layer comprise all entities in anyway 
3 5 related to or associated with the layer such as inputs, outputs, operands, 
matrices representing functions, systems, processes, methods, and 
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probability parameters. In a digital embodiment, activation results in the 
recall of the component from memory and may further result in 
processing steps such as matrix multiplication of matrices representing 
functions. Activation involves generating an activation probability 
5 parameter. The activation probability parameter is a parameter 

responsible for activating any component of the system and is dependent 
on a prior activation history of each component in the system. 

The Predominant Configuration Layer 18 includes an activation 
probability parameter generator 66. The activation probability 
1 0 parameter generator 66 receives a listing of prior activation frequencies 
of all of the available components of the present invention such as 
matrices representing functions, data parameters, Fourier components, 
Fourier series, strings, ordered strings, components of the Input Layer, 
components of the Association Layer, components of the String Ordering 

1 5 Layer, and components of the Predominant Configuration Layer from an 

activation frequency memory section 68. The activation probability 
parameter generator 66 also receives a listing of all active components 
from the processor 42. Alternatively, the activation probability 
parameter generator 66 may receive a listing of all active components 

2 0 directly from the active components. The activation probability 

parameter is stored in memory 20. The activation probability parameter 
is input to an activation probability operand generator 70. The activation 
probability operand generator 70 generates a probability operand value 
of one or zero based upon the activation probability parameter. The 

2 5 probability operand value is output to the processor 42. Any one or 

more of the components are activated when the probability operand 
corresponding to each component has a value of one as determined by 
the processor 42. Thus, the activation of each component is based on the 
corresponding activation probability parameter. Each activation 

3 0 probability parameter is weighted based on the activation rate of the 

component. The activation process continues while the system is on. 
Thus, the activation process is akin to an operating system kernel in a 
forever loop. 

Embodiments of the system for performing pattern recognition and 
3 5 processing may comprise a general purpose computer. Such a general 
purpose computer may have any number of basic configurations. For 
example, such a general purpose computer may comprise a central 
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processing unit (CPU), one or more specialized processors, system 
memory, a mass storage device such as a magnetic disk, an optical disk, 
or other storage device, an input means such as a keyboard or mouse, a 
display device, and a printer or other output device. A system 
5 implementing the present invention can also comprise a special purpose 
computer or other hardware system and all should be included within its 
scope. 

Embodiments within the scope of the present invention also include 
computer program products comprising computer readable medium 

0 having embodied therein program code means. Such computer readable 
media can be any available media which can be accessed by a general 
purpose or special purpose computer. By way of example, and not 
limitation, such computer readable media can comprise RAM, ROM, 
EPROM, CD ROM, DVD or other optical disk storage, magnetic disk storage 

5 or other magnetic storage devices, or any other medium which can 

embody the desired program code means and which can be accessed by a 
general purpose or special purpose computer. Combinations of the above 
should also be included within the scope of computer readable media. 
Program code means comprises, for example, executable instructions and 

0 data which cause a general purpose computer or special purpose 
computer to perform a certain function of a group of functions. 

The present invention may be embodied in other specific forms 
without departing from the spirit or essential attributes thereof and, 
accordingly, reference should be made to the appended claims, rather 

5 than to the foregoing specification, as indicating the scope of the 
invention. 

Also, included as part of this application is a Support Appendix and 
associated sub-appendices. These include the following: 

SUB-APPENDIX I is the derivation of the Input and the Band-Pass 
0 Filter of the Analog Fourier Processor according to the present invention; 
SUB-APPENDIX II is the derivation of the Modulation and Sampling 
Gives the Input to the Association Mechanism and Basis of Reasoning 
according to the present invention; 

SUB-APPENDIX III is the derivation of the Association Mechanism and 
5 Basis of Reasoning according to the present invention; 

SUB-APPENDIX IV is the Ordering of Associations: Matrix Method 
according to the present invention; 
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SUB-APPENDIX V is the GENOMIC DNA SEQUENCING METHOD/MATRIX 
METHOD OF ANALYSIS according to the present invention: 

SUB-APPENDIX VI is the derivation of the Input Context according to 
the present invention, and 
5 SUB-APPENDIX VII is the derivation of the Comparison of Processing 

Activity to Statistical Thermodynamics/Predominant Configuration 
according to the present invention. 

SUPPORT APPRNDTX 

1 0 

The methods and systems of the present invention are herein 
defined as the "processor" which is capable of storing, retrieving, and 
processing data to form novel conceptual content according to the present 
invention. The "processor" comprises systems and associated processes 

1 5 which serve specific functions which are collectively called "layers". The 

"layers" are organized so as to receive the appropriate inputs and 
produce the appropriate outputs according to the present invention. In a 
preferred embodiment, the memory layer is organized in a hierarchical 
manner according to the significance of the stored information. The 

2 0 significance may be measured by how frequently the information is 

recalled during processing, or it may be significant because it represents 
reference or standard information. The most significant information may 
be stored in a layer called "High Level Memory". Unlike a conventional 
processor such as a Turing Machine, the "processor" of the present 

2 5 invention may constantly change its state such that the output to a given 

input may not be identical. The "processor" may be governed by a 
principle similar to the entropy principle of thermodynamics whereby a 
chemical system achieves a state representative of a predominant 
configuration, most probable state in time. The "predominant 

3 0 configuration" of the present "processor" is the total systems of the 

"processor" and the total state of their components in time. The following 
invention of Pattern Recognition, Learning, and Processing Methods and 
Systems comprises analog or digital embodiments of: 

1.) an Input Layer which receives data representative of physical 
3 5 characteristics or representations of physical characteristics of the 
environment and transforms it into a Fourier series in £, co -space 
wherein input context is encoded in time as delays which corresponds to 
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modulation of the Fourier series at corresponding frequencies. The 
derivation of the input comprising a Fourier series in k,co -space is given 
in SUB-APPENDIX I-The Input and the Band-Pass Filter of the Analog 
Fourier Processor. The derivation of the encoding of input context in time 
5 as delays which corresponds to modulation of the Fourier series at 

corresponding frequencies is given in SUB-APPENDIX Vl-Input Context. 
A flow diagram of an exemplary transducer data structure of a time 
delay interval subdivision hierarchy is shown in FIGURE 3. The 
corresponding derivations are also given in SUB-APPENDIX VI; 
1 0 2.) an Association Filter Layer which receives multiple Fourier 

series from the Input layer, and High Level Memory, and forms a series 
(called a "string") of multiple Fourier series each representative of 
separate information by establishing "associations" between "string" 
member Fourier series. In k y co -space, the Fourier series are sampled and 

1 5 modulated via time delayed Gaussian filters called "association filters" or 

"association ensembles" that provide input to form the "associations". The 
derivation of the time delayed Gaussian filters which provide sampling 
and modulation (frequency shifting) of the Fourier series in k,co -space is 
given in SUB-APPENDIX II-Modulation and Sampling Gives the Input to 

2 0 the Association Mechanism and Basis of Reasoning. The derivation of the 

"association" of Fourier series is given SUB-APPENDIX Ill-Association 
Mechanism and Basis of Reasoning; 

3.) a "String" Ordering Layer which receives the "string" as input 
from the Association Filter Layer and orders the information represented 

2 5 by the "string" as a nested set of subsets of information with a Matrix 

Method of Analysis Algorithm via Poissonian probability based 
associations with input from High Level Memory. The methods of 
ordering the "string" comprising associated information are given SUB- 
APPENDIX IV~Ordering of Associations: Matrix Method, and 

3 0 4.) an Output of the Ordered "String" to High Level Memory Layer 

with Formation of the "Predominant Configuration" which is analogous to 
statistical thermodynamics and arises spontaneously because the 
activation of any association filter, input to the Association Filter Layer to 
form a "string", and the input to the "String" Ordering Layer are based on 
3 5 their activation history whereby activation is effected by probability 

operators. The derivation of the predominant configuration structure is 
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given in SUB-APPENDIX VII--Comparison of Processing Activity to 
Statistical Thermodynamics/Predominant Configuration. 

A flow diagram of an exemplary hierarchical relationship between 
the characteristics and the processing and storage elements of the 
5 present "processor" is shown in FIGURE 18. FIGURE 19 is a flow diagram 
of an exemplary hierarchical relationship of the signals in Fourier space 
comprising "FCs", "SFCs", "groups of SFCs", and a "string" accordance with 
the present invention. An exemplary layer structure is shown in FIGURE 
20. A flow diagram of an exemplary layer structure and exemplary 
1 0 signal format which demonstrates the relationships of the inputs and 
outputs of the processing layers is shown in FIGURE 21. 

All layers comprise processor elements called "P elements" each 
with a system function response defined as the "impulse response" (Eqs. 
(37.22-37.24)) and an output (herein defined as the "P element 

1 5 response") shown in FIGURE 6 comprising a "pulse train of impulse 

responses"— an integer number of traveling dipole waveforms (each 
called an "impulse response"). The Fourier transform of this signal is the 
convolution of a sine function with a periodic series of delta functions 
where the amplitude and the width of the sine function is determined by 

2 0 the integer number of "impulse responses" of the signal. In a preferred 

embodiment, the amplitude of the "impulse response", the temporal and 
spatial spacing or repetition frequency of the "impulse responses", and 
the integer number of "impulse responses" of the "P element" signal is 
proportional to rate of voltage change called "depolarization" of the "P 

2 5 element". This rate is determined by the amplitude and rate of change of 

the input. Thus, in the preferred embodiment, each "P element" is a 
linear differentiator— the output (pulse train of "impulse responses") is 
the sum (superposition) of the derivative of the inputs. Additionally in 
the embodiment, the "P element" has a threshold of "depolarization" to 

3 0 generate an output. In this case, the Fourier transform of "P element 

response" comprises a repeated series of a Fourier component herein 
defined as a "FC" with quantized frequency and phase angle. In another 
embodiment, the amplitude is also quantized. In k,co -space, the Fourier 
transform of the "impulse response" function filters the "FC" of a "P 
3 5 element" and is a band-pass when the spatial frequency of the "FC" is 
equal to the temporal frequency (i.e. the "FC" is band-passed when 
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An exemplary output signal of an analog "P element" to an input of 
the form given by Eq. (37.26) is given in time by Eq. (37.27) (the 
parameters p 0 , z 0 , and N may encode quantitative information such as 
intensity and rate of change of a physical parameter such as 
5 temperature) and in k,co -space by Eq. (37.32). The latter equation is that 
of a series of a Fourier component with information encoded in the 
parameters p 0 and N of the Fourier component. "P elements" are 
directionally massively interconnected in terms of the inputs and the 
outputs of the present invention which may superimpose. Multiple "P 
1 0 elements" input into any given "P element" which then outputs to 

multiple n P elements. The Fourier transform of the superposition of the 
output of multiple M P elements" is a repeating Fourier series— a repeating 
series of trigonometric functions comprising a series of Fourier 
components "FCs" herein referred to as a " SFCs". Exemplary 

1 5 representations are given by Eq. (37.33) and Eq. (37.33a). Thus, the 

present "processor" may function as an analog Fourier processor. 

All layers also comprise memory elements called "M elements" that 
store an input such as a "P element response". The stored "P element 
response" may be recalled from the "M element". Each "M element" has a 

2 0 system function response defined as the "impulse response" (Eqs. (37.22- 

37.24)) and an output (herein defined as the "M element response") also 
shown in FIGURE 6 comprising a "pulse train of impulse responses"— an 
integer number of traveling dipole waveforms (each called an "impulse 
response"). In a preferred embodiment, the output, the "M element 

2 5 response", is the product of the "pulse train of impulse responses" and a 

time ramp. In this case, the Fourier transform of "M element response" 
comprises a repeated series of a Fourier component herein defined as a 
"FC" with quantized amplitude, frequency, and phase angle. An 
exemplary output signal of a group of analog "M elements" to an input 

3 0 time ramp is given in £,<y-space by Eq. (37.33a) (the parameters p 0w , 

N "po 9 and N *n* of the ^called function are typically the same as those 
stored). The "M elements" are directionally massively interconnected in 
terms of the inputs and the outputs of the present invention which may 
superimpose. Multiple "M elements" input into any given "M element" 
3 5 which then outputs to multiple "M elements". The collective of multiple 
"M elements" including their stored inputs is referred to as "memory" of 
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the "processor". The collective storage of a signal such as a " SFCs" having 
an exemplary representation given by Eq. (37.33) to multiple "M 
elements" is called "store to memory". The collective activation of 
multiple "M elements" to provide a signal such as a " SFCs" having an 

5 exemplary representation given by Eq. (37.33a) is referred to as "recall 
from memory". An exemplary representation of information "recalled 
from memory" with input context encoded by specific modulation is 
given by Eq. (37.110). 

The Association Layer and the "String" Ordering Layer comprise 

0 cascaded processor stages which are herein defined as "stages". The 
"stages" need not be identical. Let h^t) be the impulse response of the 
stage and assume that ^.(0^0, so that the step response of each stage (or 
indeed of any number of cascaded stages) is monotonic. Cascaded stages 
form filters. The Central Limit Theorem of probability theory states in 

5 effect that, under very general conditions, the cascade of a large number 
of linear-time-invariant (LTI) systems will tend to have a delayed 
Gaussian impulse response, almost independent of the characteristics of 
the systems cascaded. Sufficient conditions of the Central Limit Theorem 
are given by Eqs. (37.52-37.55) of SUB-APPENDIX II-Modulation and 

0 Sampling Gives the Input to the Association Mechanism and Basis of 
Reasoning. The collective of multiple cascaded "stages" comprises an 
"association ensemble" that receives input such as a " SFCs". Each 
"association ensemble" serves as a heterodyne having an exemplary 
representation given by Eq. (37.50) by modulating the Fourier series in 

5 k,co -space. It further samples the Fourier series in k,co -space. The 
modulation and sampling functions correspond to a delayed Gaussian 
filter in the time domain having an exemplary representation given by 
Eq. (37.51). 

The "stages", "P elements", and "M elements" in one embodiment of 
0 the present "processor", are directionally massively interconnected in 
terms of the inputs and the outputs of the present invention which may 
superimpose. Multiple "stages", "P elements", and "M elements" input 
into any given "stage", "P element", or "M element" which then outputs to 
multiple "stages", "P elements", and "M elements". 
5 The Input Layer comprises transducers that convert physical 

signals from the environment into measurements called "data" which in 
an analog circuit embodiment, is processed into an analog time signal 
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which corresponds to a Fourier series in k 9 a> -space. In a digital 
equivalent embodiment, the "data" is further transformed by a Fourier 
transform processor into a Fourier series in k,co -space. According to the 
present invention information is encoded in a Fourier series in 

5 k,co -space. Information is not limited to that corresponding to data, but 
is meant include all forms of information such as conceptual information, 
temporal order, cause and effect relationships, size order, intensity order, 
before-after order, top-bottom order, left-right order, and knowledge 
derived from study, experience, or instruction. Data which are 

0 transducer measurements is processed into a Fourier series in k 9 a> -space 
to form input to higher layers such as the Association Layer shown in 
FIGURE 21 whereby: 

i. ) "Data" such as the intensity and the rate of change of a physical 
signal such as the surface roughness, or the intensity of sound, light, or 

5 temperature recorded by a transducer is represented in terms of the 

frequency and amplitude parameters, p 0w and N , of each component of 

the Fourier series (e.g. Eq. (37.33a)). Information is represented in terms 
of the parameters p 0m and of each component of the Fourier series in 

the sense that if the transducer and Fourier processor were each a 
0 reciprocal device, then inputting the Fourier series into the output of the 
Fourier transform processor would yield the measured physical signals at 
the input of the transducers. 

ii. ) The input from the Input Layer to other layers can be an analog 
waveform in the analog case and a matrix in the digital case. Input 

5 context of a given transducer can be encoded in time as delays which 
correspond to modulation of the Fourier series in k,co -space at 
corresponding frequencies whereby the data corresponding to each 
transducer maps to a distinct memory location called a "register" that 
encodes the input context by recording the data to corresponding specific 

0 time intervals of a time division structured memory. 

iii. ) Input context of a complex transducer system can be encoded 
in time by the mapping of data from the components of the transducer 
system to a memory structured according to a corresponding hierarchical 
set of time intervals representative of each transducer system with 

5 respect to different transducer systems, a transducer element's rank 

relationship relative to other transducer elements, and the response of a 
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transducer element as a function of time. In terms of digital processing, 
the data from a transducer having n levels of subcomponents is assigned 
a master time interval with « + i sub time intervals in a hierarchical 
manner wherein the data stream from the final n th level transducer 
5 element is recorded as a function of time in the n + 1 th time coded sub 
memory buffer. During processing the time intervals represent time 
delays which are transformed into modulation frequencies which encode 
input context. A flow diagram of an exemplary transducer data structure 
of a time delay interval subdivision hierarchy is shown in FIGURE 3. An 
1 0 exemplary complex transducer which may represented by a data 

structure comprising a hierarchical set of time delay intervals is a video 
camera which is comprised of a multitude of charge coupled devices 
(CCDs), transducer elements each responsive to light intensity of a given 
wavelength band at a given spatial location in a grid. Another example is 

1 5 an audio recorder comprising transducer elements each responsive to 

sound intensity of a given frequency band at a given spatial location or 
orientation. A signal within the band 300-400 MHz may encode and 
identify the signal as a video signal; whereas, a signal within the band 
500-600 MHz may encode and identify the signal as an acoustic signal. 

2 0 Furthermore, a video signal within the band 315-325 MHz may encode 

and identify the signal as a video signal as a function of time of CCD 
element (100,13) of a 512 by 512 array of CCDs. An exemplary 
representation of a " SFCs" output of "P elements" or "M elements" with 
input context encoded by specific modulation is given by Eq. (37.110); 
25 iv -) Th e relationship between the "data" and the parameters p 0 and 

N m n of each component of the Fourier series, may be learned by the " 
"processor" by applying standard physical signals to each transducer 
together with other information that is "associated" with the standard. 
The information that is "associated" with the standard can be recalled and 

3 0 may comprise input to the Association Layer and the "String" Ordering 

Layer during processing according to the present invention. 

The Association Filter Layer receives multiple Fourier series from 
the Input Layer, and High Level Memory, and forms a series (called a 
"string") of multiple Fourier series each representative of separate 
3 5 information by establishing associations between "string" member 

Fourier series. The "association" between one or more Fourier series that 
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form the "string" occurs with Poissonian probability based on the spectral 
similarity of each association filtered Fourier series member with that of 
one or more others filtered by the same or different association filters as 
described further below. 
5 The process of storing output from multiple transducers to memory 

further comprises creation of "transducer strings". In one embodiment of 
this case, associations occur at the transducer level, and "SFCs" are 
mapped to distinct "registers" from the corresponding distinct 
transducers responding simultaneously, for example. Consider a 

1 0 "transducer string" made up of multiple "groups of SFCs" where each 

"SFCs" represents information of the transducer system with respect to 
different transducer systems, a transducer element's rank relationship 
relative to other transducer elements, and the response of a transducer 
element as a function of time. These aspects of each transducer are 
15 encoded via time delays corresponding to modulation in k y co -space 

within a frequency band corresponding to each aspect of the transducer. 

Two or more "transducer string" Fourier series such as two or more 
"SFCs" may become "linked" which is defined according to a 
corresponding linkage probability weighting parameter wherein 

2 0 activation of one "string" Fourier series may cause other "string" Fourier 

series to become active in according to the linkage probability weighting 
parameter. The probability that other "string" Fourier series are 
activated when any given "string" Fourier series is activated defines the 
"linkage". "Active" in this case of an analog embodiment is defined as 

2 5 providing an output signal; thus, "activate" is defined as causing an 

output signal. "Active" in a digital embodiment is defined as recalled 
from memory; thus, "activate" is defined according to causing a Fourier 
series to be recalled from memory. 

In a general sense, the "string" in k, co -space is analogous to a 

3 0 multidimensional Fourier series. The modulation within each frequency 

band may encode context in a general sense. In one embodiment, it 
encodes temporal order, cause and effect relationships, size order, 
intensity order, before-after order, top-bottom order, left-right order, etc. 
which is relative to the transducer. Further associations are established 
3 5 between "groups of SFCs" (i.e. a new "string" is created) by the 
Association Filter Layer. 
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The Association Filter Layer receives multiple Fourier series from 
the Input layer, and High Level Memory, and forms a series (called a 
"string") of multiple Fourier series each representative of separate 
information by establishing "associations" between "string" member 
5 Fourier series. FIGURE 19 is a flow diagram of an exemplary hierarchical 
relationship of the signals in Fourier space comprising "FCs", "SFCs", 
"groups of SFCs", and a "string" in accordance with the present invention. 
Each "FC" is "carried" (processed as a response to an input) by a "P 
element" or stored into and/or recalled from a "M element" as shown in 

10 FIGURE 18 which is a flow diagram of an exemplary hierarchical 

relationship between the characteristics and the processing and storage 
elements of the present "processor". Each Fourier series such as a "SFCs" 
representing information is filtered and delayed in the time domain 
(modulated and sampled in the frequency domain or k,co -space) as it is 

15 recalled from memory and "carried" (processed as a response to the 
memory input) by a series of cascaded association "stages" called an 
"association ensemble" or "association filter". Since the Fourier series is in 
k, co —space, the modulation corresponds to a frequency shift. Each 
"association ensemble" is weakly linked with multiple other "association 

2 0 ensembles" at the level of the "stages". The "association ensembles" 

produce interference or "coupling" of the "SFCs" of one set of "stages" with 
that of another by producing frequency matched and phase locked 
Fourier series -sums of trigonometric waves that are frequency matched 
and periodically in phase-that give rise to "association" of the 

2 5 corresponding recalled or prior processed information. 

"Coupling" gives rise to the formation of "associations" between one 
or more Fourier series that form the "string". "Coupling" refers to 
interference or energy exchange between "association ensembles" in an 
analog embodiment. In a digital embodiment, "coupling" refers to 

3 0 calculating an "association" probability parameter based on the spectral 

similarity of the each Fourier series such as a "SFCs" filtered by an 
"association filter" with that of one or more other Fourier series filtered 
by the same or different "association filters". The statistics may be 
Poissonian. "Association" refers to recording "coupled" Fourier series to 
3 5 memory based on the probability of the "coupling". In a digital 

embodiment, "association" refers to marking two or more Fourier series 
as associated based on a zero or one outcome of a probability operand 
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applied to the "association" probability parameter and recording the 
"associated" Fourier series to memory. The "association" probability 
parameter based on Poissonian probability is derived from a correlation 
function in SUB-APPENDIX Ill-Association Mechanism and Basis of 
5 Reasoning. The "association" probability parameter has a "coupling cross 
section" amplitude and a "frequency difference angle" as parameters. 
The former is a weighting parameter of the spectral similarity of Fourier 
series which may become "associated". The "frequency difference angle" 
is the fractional difference in the frequencies of the Fourier series which 
1 0 may become "associated" expressed as an angle. The derivation of these 
parameters as well as the derivation of the "association" of Fourier series 
that "couple" with Poissonian probability is also found in SUB-APPENDIX 
III. 

In a preferred embodiment, the "string" is formed by the 

1 5 Association Filter Layer with input context. In this case, "association" 

occurs whereby the "SFCs" or "groups of SFCs" such as those comprising 
"transducer strings" comprise a transducer specific frequency modulation 
factor. Exemplary representations of "string" outputs of "P elements" or 
"M elements" with input context encoded by specific modulation are 

2 0 given by Eq. (37.114) and Eq. (37.115). In this case, an exemplary 

representation of the "coupling cross section" amplitude and the 
"frequency difference angle" based on the spectral similarity of the each 
"SFCs" filtered by an "association filter" with that of one or more other 
"SFCs" filtered by the same or different "association filters" is given by Eq. 

2 5 (37.111) and Eq. (37.112). 

The "String" Ordering Layer receives the "string" as input from the 
Association Filter Layer and orders the information represented in the 
"string" via a method developed by Mills for sequencing DNA called the 
"Matrix Method" which is herein presented as a mechanism used by the 

3 0 "processor" to sequence information temporally, conceptually, or 

according to causality. First, the "string" (multiple Fourier series) is 
stored in memory. The "string" is recalled and processed by further sets 
of specific "association ensembles" that "couple" with other "higher level 
associations", information with conceptual significance established by a 
3 5 previous execution of the present procedure. In it, co -space, each specific 
"association ensemble" samples the "string", a Fourier series in jfc, co -space. 
It also serves as a heterodyne by modulating the Fourier series in 
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Kco -space. The sampling in the frequency domain is dependent on the 
particular half-width parameter, a s9 of each specific "association 
ensemble". The collective sampling of the specific "association 
ensembles" provides a nested set of subsets of information where each 
5 subset maps to a specific time point corresponding to the specific delay, 

— K of the specific Gaussian filter of the "association ensemble" (Eqs. 

s 

(37.50-37.51)). The nested set of subsets of information is ordered by 
the Matrix Method of Analysis Algorithm of Mills with Poissonian 
probability based associations with input from High Level Memory. Each 
10 "group of SFCs" of the input "string" has the corresponding time delay 

parameter, and the half- width parameter, a s , of the Gaussian filter 

of the "association ensemble" (Eqs. (37.50-37.51)) that resulted in the 
"coupling" and "association" to form the "string". The process of ordering 

J If 

assigns a particular time delay, and half- width parameter, a y , to 

1 5 each "group of SFCs" of the output "string". The half-width parameter, a s . 
corresponds to each specific delayed Gaussian filter that samples the 
input "string" in the frequency domain to provide each "group of SFCs" of 
the output "ordered string". Each corresponding particular time delay, 

— , encodes and corresponds to the time domain order of each "group 



a, 



2 0 of SFCs" of the output "ordered string". An order processed "string" called 
a "P string" may comprise complex information having conceptual 
content. 

The Output of the Ordered "String" to High Level Memory Layer 
with Formation of the Predominant Configuration receives ordered 

2 5 strings from the High Level Memory and forms more complex ordered 

strings as shown in FIGURE 20. This layer also activates components 
within other layers. The Output of the Ordered "String" to High Level 
Memory Layer with Formation of the Predominant Configuration is 
analogous to statistical thermodynamics and arises spontaneously 

3 0 because the activation of any "processor" component such as any "P 

element", "M element", "stage", "association ensemble", "SFCs", "string", 
"ordered string", "transducer string" having "linkages", Fourier series 
"linkage", input to the Association Filter Layer to form a "string", and the 
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input to the "String" Ordering Layer are based on their past activation 
frequency whereby activation is effected by probability operators. In 
one embodiment, an activation probability parameter is generated and 
stored in memory for each "processor" component. A probability operand 
5 is generated having a value selected from a set of zero and one, based on 
the activation probability parameter. If the value is one, the component 
is activated. Thus, any "processor" component is randomly activated 
wherein the activation is based on the activation probability parameter. 
The activation probability parameter is weighted based on an activation 
1 0 rate. "Processor" components may become "linked" which is defined 

according to the corresponding probability weighting parameter wherein 
activation of one "processor" component may cause other "processor" 
components to become "active" according to the probability weighting 
parameter. The probability that other "processor" components are 

1 5 activated when any given "processor" component is activated defines the 

"linkage". 

The processing of information depends on and dynamically alters 
(through feedback) the total state of stored information, the cascades of 
association "stages", and the hierarchical relationships of association 

2 0 "stages" and stored information (memory). "Memory linkages" occur 

whereby recalling any part of a string from a distinct memory location 
thereby causes additional Fourier series of the string to be recalled. 
"Linkages" between "stages" occur whereby activating any "stage" 
thereby causes additional "stages" to become "active". A strongly "linked 

2 5 group of cascaded association "stages" comprises an "association 

ensemble", and a strongly linked group of memory elements comprises a 
"memory ensemble". Repetitive activation of a "memory element" or 
association "stage" increases the probability of its future activation. A 
configuration of "couplings" between "memory ensembles" and 

3 0 "association ensembles" increases the probability of future activation of 

the configuration. Analogously to statistical thermodynamics, a 
predominant configuration arises from the ensemble level. Consider the 
"processor" on a higher level. The activation history of each ensemble 
relates to a hierarchical relationship of coupled "memory and association 
3 5 ensembles" which gives rise to a precedence of higher order predominant 
configurations. Pattern recognition, learning, and the ability to associate 
information and create novel information is a consequence. Machine 
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learning arises by the feedback loop of transducer input to the coupled 
predominant configurations which increases the basis for machine 
intelligence. 

Pattern recognition and learning arise from the massive directional 
5 connectivity in terms of the output to input relationships of the 

"processor" which in one embodiment functions as an analog Fourier 
processor wherein a superposition of "P element responses" becomes a 
superposition of trigonometric functions in frequency space (£, co -space). 
Information is digitized in amplitude, frequency, and phase in k,co -space 

1 0 via the "P element response". It is then modulated, sampled, associated, 

and ordered via the properties of cascaded groups of association "stages" 
with "couplings" governed by Poissonian probability. For the "processor", 
since information is encoded in Fourier series in k,co -space, specific time 
delays achieve the specific modulations equivalent to that of heterodynes 
15 of conventional signal processing circuits. In other words, a clock 

substitutes for a multitude of heterodyne circuits to encode input context 
wherein aspects of each transducer are encoded via time delays 
corresponding to modulation in k, co -space within a frequency band 
corresponding to each aspect of the transducer. The modulation and 

2 0 sampling functions correspond to a delay and filter (delayed Gaussian 

filter) in the time domain analogous to the key components of amplitude 
modulation (AM) radio except that the Fourier series of the signal and its 
modulation occurs co -space in the case of the present "processor" 
versus the time domain in the AM signal processing case. The filtering 

2 5 function occurs in the time domain in both cases. The unique processing 

features of the "processor" further permits ordering of information by a 
method developed by Mills for sequencing DNA called the "Matrix 
Method" which is herein presented as a mechanism used by the 
"processor" to sequence information temporally, conceptually, or 

3 0 according to causality. According to the Fourier theorem any waveform 

can be recreated by an infinite series of trigonometric functions, and any 
aspect of the universe can be represented by an infinite series of sine 
and cosine functions as processed by the "processor". For the "processor" 
of the present invention, the trigonometric function is the basis element 
3 5 of information. The quantity of information such as "inputs" that can be 
associated into ordered "strings" ("P strings") is essentially infinite based 
on it being encoded in Fourier series in k,co -space. And, the number of 
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terms necessary to represent most objects is not overwhelming. In fact, 
even a potentially challenging object having sharp edges such as a square 
pulse poses no difficulty in that is fairly accurately represented by only 
seven terms of a Fourier series in the time domain comprising the prior 
5 art [1], The same principle applies to information represented as a 
Fourier series in k, co -space. 

The following invention of Pattern Recognition, Learning, and 
Processing Methods and Systems comprises analog or digital 
embodiments. In one embodiment, analog circuit elements store, 
1 0 retrieve, and process input waveforms wherein the circuit elements have 
the system functions or impulse responses or comprise the operators and 
structures which transform input to output as described herein. In 
another embodiment, the mathematical functions corresponding to the 
waveforms of any stage of storage, retrieval, or processing are 

1 5 represented digitally, and the digital waveforms are digitally processed 

in a manner equivalent to the analog embodiment according to signal 
processing theory such as the Nyquist theorem. In a preferred 
embodiment, a digitally based "processor" comprises simulations methods 
and systems according to the analog systems and processes of the present 

2 0 invention. The Nyquist theorem states that all of the information in any 

waveform can be conserved and recovered by digital processing with 
frequency components equal to twice the maximum frequency of any 
waveform [2]. Thus, the analog and digital embodiments perform 
equivalently. 

25 

Exemplary Layer Structure and Exemplary Signal Format 

FIGURE 20 shows an exemplary layer structure in accordance with 
the present invention. FIGURE 21 shows a flow diagram of an exemplary 

3 0 layer structure and exemplary signal format in accordance with the 

present invention. The present invention comprises an analog Fourier 
"processor" wherein the basis element of information in k,co- space is the 
Fourier component. In a preferred embodiment, the analog systems and 
processes are implemented using the corresponding digital embodiments. 
3 5 The "processor" is applicable to standard computers comprising digital 
processors, and digital memory, storage, and retrieval systems where 
discrete values of the continuous functions evaluated at selected 



WO 99/34322 



PCT/US98/27624 



38 

frequencies and/or at the Nyquist rate [2] form matrices upon which the 
operations of the exemplary signal format are performed in place of the 
continuous functions. Exemplary embodiments of the present invention 
according to the layer structure of FIGURE 20 and the exemplary layer 
structure and exemplary signal format of FIGURE 21 comprises: 

Input Layer 



The Input Layer receives data and transforms it into a Fourier 
10 series in fc,a>-space wherein input context is encoded in time as delays 
which corresponds to modulation of the Fourier series at corresponding 
frequencies. Data is processed into a Fourier series in jfc, co -space that 
represents information as given by Eq. (37.33) and Eq. (37.33a) 



4k 



-a 0 sin 



k p ~n 



In 



sin 



2k 



0J 



1 5 



(37.33) 



Hi 



M ~ 4k 4 



k 2 



. Po m ^o M1 



a — a ± sin 



k p -n 



2k 



J o m J 



sin 



k—n 



2k )N„Zo m 



M CO 



4k 



a n N N sin 



^2 "0« i¥ m po iv m 2n 



k p -n 



2k 



Po, 



'« J 



sin 



k z —n 



2k 

Vo„ J 



(37.33a) 



2 0 whereby i.) data such as intensity and rate of change recorded by a 

transducer is represented in terms of the parameters p Qm and N of each 

component of the Fourier series; ii.) input context is encoded in time by a 
hierarchical set of time delay intervals representative of each transducer 
system with respect to different transducer systems, a transducer 
2 5 element's rank relationship relative to other transducer elements, and 
the response of a transducer element as a function of time, and iii.) the 
input from the Input Layer to other layers shown in FIGURE 21 can be an 
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analog waveform in the analog case and a matrix in the digital case 
wherein input context of a given transducer can be encoded in time as 
delays which correspond to modulation of the Fourier series in k,(0 -space 

,- jk p {Pfi>,.. + P;.,„) 
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whereby the data corresponding to each transducer maps to a distinct 
memory location called a "register" that encodes the input context by 
recording the data to corresponding specific time intervals of a time 
division structured memory, and iv.) the relationship between the "data" 
and the parameters p 0m and N mpQ of each component of the Fourier series, 

may be learned by the "processor" by applying standard physical signals 
to each transducer together with other information that is associated with 
the standard. The information that is "associated" with the standard can 
be recalled and may comprise input into the Association Layer and the 
"String" Ordering Layer during processing according to the present 
invention. In terms of digital processing, the data from a transducer 
having n levels of subcomponents is assigned a master time interval with 
n + l sub time intervals in a hierarchical manner wherein the data stream 
from the final n th level transducer element is recorded as a function of 
time in the n + l th time coded memory buffer. During processing the 
time intervals represent time delays which are transformed into 
modulation frequencies which encode the input context. FIGURE 3 is a 
flow diagram of an exemplary transducer data structure of a time delay 
interval subdivision hierarchy wherein the data from a transducer 
having n levels of subcomponents numbering integer m per level is 
assigned a master time interval with n + l sub time intervals in a 
hierarchical manner wherein the data stream from the final n th level 
transducer element is recorded as a function of time in the n + l th time 
coded sub memory buffer in accordance with the present invention. 

The process of storing output from multiple transducers to memory 
further comprises creation of "transducer strings". In one embodiment, 
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associations occur at the transducer level, and "SFCs" are mapped to 
distinct "registers" from the corresponding distinct transducers 
responding simultaneously, for example. Consider a "transducer string" 
made up of multiple "groups of SFCs" where each "SFCs" represents 
5 information of the transducer system with respect to different transducer 
systems, a transducer element's rank relationship relative to other 
transducer elements, and the response of a transducer element as a 
function of time. These aspects of each transducer are encoded via 
delays corresponding to modulation in k,co -space (Eq. (37.109)) within a 
1 0 frequency band corresponding to each aspect of the transducer. 

x(t) = j x(f¥ 2 *df x(t) = J x{ty j2Kft dt 

-<*> -oo 

(37.109) 

Delay ^(*~*o) <=> e -J 2 *fto 

Two or more "transducer string" Fourier series such as two or more 

1 5 "SFCs" may become "linked" which is defined according to a 

corresponding linkage probability weighting parameter wherein 
activation of one "string" Fourier series may cause other "string" Fourier 
series to become "active" according to the linkage probability weighting 
parameter. The probability that other "string" Fourier series are 

2 0 activated when any given "string" Fourier series is activated defines the 

"linkage". 

The "string" in k,a> -space is analogous to a multidimensional 
Fourier series. The modulation within each frequency band may further 
encode context in a general sense. In one embodiment, it encodes 

2 5 temporal order, cause and effect relationships, size order, intensity order, 

before-after order, top-bottom order, left-right order, etc. which is 
relative to the transducer. 

A "FC" of Eq. (3732) is a series of a Fourier component. A distinct 

3 0 superposition or series of "FCs" is called a "SFCs" which further 

superimpose to form "groups of SFCs". The data is digitized according 
to the parameter N of Eqs. (37.33), (37.33a), and (37.87). Input to 
higher layers is in a Fourier series format in k,co -space or data is 
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processed with a FFT (Fast Fourier Transform) routine and stored in 
memory as a series of a Fourier component in it, co -space with 
quantized amplitude, frequency, and phase angle (Eq. (37.33a)). Or, 
data is processed with a FFT (Fast Fourier Transform) routine and 
5 stored in memory as a series of a Fourier component in k, co -space 

with quantized frequency, and phase angle of the form of Eq. (37.33). 
In this case, "groups of SFCs" representing information are recalled 
from memory with a time ramp multiplication of each "FC" of a "SFCs" 
to give the form of Eq. (37.33a). In the digital case, multiplication is 
10 performed via multiplication of corresponding matrices formed from 
the continuous functions by evaluating them at discrete frequency 
values. A summary of an exemplary method of inputting data follows: 

a. ) data is recorded by one or more transducers each having one or 

1 5 more levels of component elements; 

b. ) the data recorded by each transducer is encoded as parameters 
such as p 0m and //„ % of a Fourier series in Fourier space with input 

context representing the information based on the physical 

2 0 characteristics and the physical context; 

c. ) the data from a transducer having n levels of subcomponents is 
assigned a master time interval with n + l sub time intervals in a 
hierarchical manner wherein the data stream from the final n th level 

2 5 transducer element is recorded as a function of time in the n + l th 

time coded memory buffer; 



d.) the time intervals represent time delays which are transformed 
into modulation frequencies which encode input context (e.g. the 
transducer element relationship of more than one transducer 
elements, its rank in the transducer hierarchy, and the time point of 
data recording); 



e.) the representation of the data is given by Eq. (37.110) 
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(37.110) 
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f.) in the digital case, the function of Eq. (37.110) comprising a 
"SFCs" is evaluated at discrete frequencies at twice the rate of the 



g. ) "SFCs" are mapped to distinct "registers" from corresponding 
distinct transducers responding simultaneously to form "transducer 
strings" having a representation given by Eq. (37.113) wherein input 
context is encoded by the transducer modulation factor e ~ ikp ^ +Pt -) ; 

h. ) in the digital case comprising "memory linkages" of a 
"transducer string", recalling any part of a "transducer string" from a 
distinct memory location may thereby cause additional "linked" 
Fourier series of the "transducer string" to be recalled. In one 
embodiment, a linkage probability parameter is generated and stored 
in memory for each "string" Fourier series such as a "SFCs". A 
probability operand is generated having a value selected from a set of 
zero and one, based on the linkage probability parameter. If the value 
is one, the corresponding Fourier series is recalled. Thus, when any 
part of a "transducer string" is recalled from memory, any other 
"string" Fourier series is randomly recalled wherein the recalling is 
based on the linkage probability parameter. The linkage probability 
parameter is weighted based on the linkage rate. 

Association Filter Layer to Form a "String" 

Each "SFCs" is filtered and delayed in the time domain (modulated 
and sampled in the frequency domain) as it is processed by a cascade of 
association filters (subprograms in the digital case) called an "association 
ensemble". Each "association ensemble" is weakly linked with multiple 



highest discrete frequency ( 



) to form a matrix for each "SFCs"; 



2 
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other such "association ensembles". These "association ensembles" 
produce interference or "coupling" of one "SFCs" with another by 
producing frequency matched and phase locked Fourier series —sums of 
trigonometric waves that are frequency matched and periodically in 
phase— that give rise to "association" of the recalled or prior processed 
information "carried" by the cascade. The Poissonian probability of such 
"association" (Eq. (37.106c)) is given by a correlation function given in the 
SUB-APPENDIX III— Association Mechanism and Basis of Reasoning 
wherein Eq. (37.87) and Eq. (37.89) are parameters. 



cos(5j +25111^) 



(37.106c) 



2 1 



2k y a,+^ 



m l =1 



til j =1 



2 2 



2v„ 



2v„ 



(37.87c) 
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2v„ ^ 



2v 



(37.89) 



The set of "associated" "groups of SFCs" is herein called a "string". The 
"string" comprises a Fourier series, a linear sum of "FCs". FIGURE 19 is a 
flow diagram of an exemplary hierarchical relationship of the signals in 
2 0 Fourier space comprising "FCs", "SFCs", "groups of SFCs", and a "string" in 
accordance with the present invention. Each "FC" is encoded by a "P 
element" or stored into and/or recalled from a "M element" as shown in 
FIGURE 18 which is a flow diagram of an exemplary hierarchical 
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relationship between the characteristics and the processing and storage 
elements of the present "processor". 

A summary of an exemplary method of establishing "associations" 
between "groups of SFCs" (i.e. a creating a "string" ) by "coupling" with 
Poissonian probability between "association ensembles" "carrying" the 
"groups of SFCs" comprising a transducer frequency band modulation 
factor according to Eq. (37.110) follows: 

a. ) n (n an integer) inputs each comprising a "SFCs", the function of 
Eq. (37.110) which in the digital case is evaluated at discrete 

frequencies at twice the rate of the highest discrete frequency (— y^) 
to form a "SFCs" matrix, is recalled from memory; 

b. ) in the digital case, discrete values are determined at twice the 

N p 

rate of the highest discrete frequency ( ) of the Fourier series 

inputs of up to n different Fourier transforms of delayed Gaussian 
filters functions (37.50) to form up to n different association filter 
matrices; 

c. ) in the digital case, the discrete values of each of n (n an integer) 
inputs each comprising a "SFCs", the function of Eq. (37.110) which is 
evaluated at discrete frequencies to form a "SFCs" matrix, are 
multiplied on a matrix element by matrix element basis corresponding 
to the same frequency with one or more of the n different association 
filter matrices each comprising the Fourier transform of a delayed 
Gaussian filter (Eq. (37.50)); 

H N {f)»e*">e Uj (37.50) 

d. ) the "coupling cross section" amplitude, & 2 , and frequency 
difference angle, </> s , of the harmonic "coupling", is calculated for two or 
more filtered inputs. In the case of input context, the amplitude, &\ 
which follows from Eq. (37.87c) is given by Eq. (37.111b), and the 



20 
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frequency difference angle, 0 5 , which follows from Eq. (37.89) is given 
by Eq. (37.112a); 



exp- 



2 1 



a, 2 + a 1 . 
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(37.111b) 
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(37.112a) 



1 0 



e.) the Poissonian probability of "association" is calculated (Eq. 
(37.106c)) with the "coupling cross section" amplitude, ft 2 , and 
frequency difference angle, <j> s , as parameters; 



1 5 



f.) a Poissonian probability operand with the expectation value 
given by the Poissonian probability of "association" (step e) is 
activated to return a value of zero or one; 



20 



g. ) if the output of the Poissonian probability operand is one, then 
the two or more filtered inputs are marked as "associated" and this 
status is stored in memory; 

h. ) the process of forming "associations" (Steps a-g) are repeated 
including processing the "SFCs" inputs and "associated" "SFCs" inputs 
with multiple "association ensembles" comprising Gaussian filters each 

of different delay, and half-width parameter, a s to extend the 

2 5 number of associated "SFCs" to form a string; 
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i.) in one analog embodiment, the output in Fourier space is the 

m 

"string" given by Eq. (37.113) comprising the superposition of S 
"groups of SFCs" wherein each "SFCs" corresponds to the response of M 
"M or P elements", with input context. In another embodiment, the 
output is the "string" of Eq. (37.114) 



s=l m-l n= — 



J >JE„«.,, 



(f 



k p -n 



2k 



f f 



sin 



k—n- 



2k 



(37.114) 

wherein each "SFCs" is multiplied by the Fourier transform of the 
delayed Gaussian filter (Eq. (37.50)) (i.e. the modulation factor 

4(^) 2 -i(^y -j^u^\ u . u 

e J e y a) e 2K a) e ^ a) ) which gave rise to "coupling" and 
"association" to form the "string". In the digital case, the output in 

m 

Fourier space is the "string" given by Eq. (37.113) comprising the 
superposition of 5 "groups of SFCs" wherein each "SFCs" corresponds to 
a matrix digitized according to Eq. (37.110), with input context. In 
another embodiment of the digital case, the output is the "string" 

Ifl 

of Eq. (37.114) wherein each "SFCs" corresponds to a matrix digitized 
according to Eq. (37.110) that is multiplied by a digitized matrix 
according to the Fourier transform of the delayed Gaussian filter (Eq. 
(37.50)) which gave rise to the "coupling" and "association" to form the 
"string". 



"String" Ordering Layer 



25 



The "string" representing information is temporally or conceptually 
ordered via the Matrix Method of Analysis of Mills [3, 4]. Each "group of 
SFCs" of the input "string" has the corresponding time delay parameter, 

, and the half-width parameter, a s , of the Gaussian filter of the 



a. 
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"association ensemble" (Eq. (37.51)) that resulted in the "coupling" and 
"association" to form the "string". 



■-fj 



MO' 



a 



(37.51) 



The "string" comprises a Fourier series, a linear sum of "FCs" each 
5 multiplied by its corresponding Gaussian filter modulation factor and 

modulation factor which encodes input context (Eq. (37.114)). Therefore, 
new series of "FCs", "SFCs" or "groups of SFCs", may be formed using 
additional "association filters" that sample the input "sting" in Kw- space. 
In a preferred embodiment, the string is sampled with specific 
1 0 "association ensembles" which provide a "nested set of subsets" of 
information comprised of a "SFCs" and "groups of SFCs" where each 
"subset" sampled from the input "string" maps to a specific time point 



corresponding to the specific delay, of the specific Gaussian filter of 

the "association ensemble" (Eqs. (37.50-37.51)). The process of ordering 

Jff 

15 assigns a particular time delay, — , and half-width parameter, a jM to 

each "subset" of the output "string" using the "nested set of subsets" as 
input to the Matrix Method which is herein presented as a mechanism 
used by the "processor" to sequence information temporally, 
conceptually, or according to causality. 
2 0 Consider Eqs. (37.33) and (37.33a) which represent a "SFCs" in 

k,co- space comprising a Fourier series. A "string" is a sum of Fourier 
series which follows from Eqs. (37.33) and (37.33a) and is given by Eqs. 
(37.107) and (37.108). 
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(37.108) 

The corresponding equations in the time domain are a sum of multiple 
finite series of traveling dipoles (each an "impulse response") wherein 
5 each dipole series is periodic in space and time. In frequency space, each 
time delayed Gaussian filter ("association ensemble" corresponding to a 
"SFCs") modulates and samples the Fourier series representing 
information. Thus, the time delayed Gaussian filter selects information 
from the "string" and provides input for the association mechanism as the 
1 0 "processor" implements the Matrix Method of Analysis to find the order 
of the associated pieces of information represented by each "SFCs" or 
"group of SFCs" of the "string". 

Consider the time interval * = /. to t = t f of a "string" associated by 

"association ensembles" and recorded to memory. By processing the 

1 5 "string 1 with multiple "association ensembles" comprising Gaussian filters 

each of different delay, -^-> and half-width parameter, a s , the "string" 

can be "broken" into "groups of SFCs" each having a center of mass at a 
time point corresponding to the delay and frequency composition 

corresponding to a 5 which form a nested set of "sequential subsets" of 

2 0 "groups of SFCs" of the "string" in £,o>-space which map to time points 

which are randomly positioned along the time interval from the /^/,-side 

and the r = r 2 -side as shown in FIGURES 8, 10, 12, and 14. This nested set 

of "sequential subsets" of random "groups of SFCs" mapping to random 
time points from the r = r r side and the / = reside is analogous to the 

2 5 nested set of "sequential subsets" of random DNA fragments from the 5' 
end and the 3' end. The order in both cases can be solved by the 
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Genomic DNA Sequencing Method/Matrix Method of Analysis of Mills [3, 
4] described in SUB-APPENDIX V. 

The output of an association filter is the convolution of the input 
"groups of SFCs n ( each "SFCs" is given by Eqs. (3733) and (37.33a)) of a 
5 "string" (Eq. 37.108) or the string itself with a delayed Gaussian. In 
terms of the matrix method of analysis (hereafter "MMA"), the filter 
parameter a of the time delayed Gaussian filter corresponds to the 
acquisition of the composition of a polynucleotide member of a nested set 
of subsets. The time delay (time domain) and modulation (frequency 



1 0 domain) parameter J— determines the center of mass of the output, and 



it corresponds to the terminal nucleotide data. By forming "associations" 
with input from "High Level Memory", the "processor" determines the 
relative position of the center of mass of each Fourier series such as a 
"group of SFCs" as either "before" or "after" the center of mass of the 

1 5 preceding and succeeding Fourier series "associated" with Fourier series 

input from "High Level Memory". The complete set of Fourier series 
"associated" with Fourier series input from "High Level Memory" covers 
all of the frequencies of the "string". By Parseval's theorem, by 
processing the entire interval in k, co -space, the information is entirely 

2 0 processed in the time domain. The order such as temporal order of the 

Fourier series representing information is determined using the MMA. 

"Groups of SFCs" such as the "groups of SFCs" represented by Eq. 
(37.110) comprising a transducer frequency band modulation factor 
"carried" by "association ensembles" "couple" with Poissonian probability. 

2 5 "Associations" are established between "groups of SFCs" that result in the 

output of a second ordered "string" created from the input "string". In 
this case of input context, the "coupling cross section" amplitude, & 2 , 
which follows from Eq. (37.87) is given by Eq. (37.111). And, the 
frequency difference angle, 0 S , of the "coupling" which follows from Eq. 

3 0 (37.89) is given by Eq. (37.112a). 

Input to form "associations" is provided by changing the decay 
constant a and the number of "stages" in the cascade N y or by processing 
"a SFCs" of a "string" using an "association ensemble" with different 
parameters a and N over all "groups of SFCs" that make up the entire 
3 5 "string". Each "group of SFCs" is determined to be on the / = f r side or the 
r = r r -side of the "axis" of the "string" corresponding to the S'-side or 3'- 




WO 99/34322 



PCT/US98/27624 



50 

side of the "axis" of a polynucleotide to be sequenced via the Matrix 
Method of Analysis. A feedback loop comprises sequentially switching to 
different "known", "set", or "hardwired" delayed Gaussian filters which 
corresponds to changing the decay constant, a s , with a concomitant 
5 change in the half-width parameter, a s , and the number of elements, N s , 

with a concomitant change in the delay, where each a and ^ is 

"known" from past experiences and associations. The feedback loop 
whereby information from memory encoded in the "string" is filtered and 
delayed (modulated and sampled in frequency space) to provide "FCs", 
1 0 "SFCs" or "groups of SFCs" which are "associated" with input from "High 
Level Memory" provides the data acquisition and processing equivalent 
to the formation, acquisition, and analysis of the composition and 
terminal nucleotide data of a set of "sequential subsets" of the Matrix 
Method of Analysis. Changing the filters which process the "string" 

1 5 corresponds to changing the "guess" of the "known" nucleotides, 

K X K 2 K 3 K A -K n ., as well as the "unknown" nucleotides, A:,,X 2 ,X 3 ,X 4 .., of the 
Matrix Method of Analysis as applied to DNA sequencing. The order of 
the "groups of SFCs" of the "string" is established when "associations" with 
the "High Level Memory" are achieved for a given set of delayed Gaussian 

2 0 filters. Then each Fourier series of the ordered "string" is recorded to the 

"High Level Memory" wherein each Fourier series of the ordered "string" 
may be multiplied by the Fourier transform of the delayed Gaussian filter 

represented by ;^ k \<^ that established the 

correct order to form the ordered "string". The total output response V y 

M 

2 5 in Fourier space comprising the superposition of S "groups of SFCs" 

wherein each "SFCs" corresponds to the response of M "M or P elements", 
with input context, is the "string" given by Eq. (37.113). 

A summary of a method of ordering the nested set of subsets of 

3 0 Fourier series (e.g. each a "group of SFCs") follows: 

a.) the "string" of the Association Filter Layer to Form a "String 
Section is recalled from memory; 
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b. ) the recalled "string" is filtered and delayed (modulated and 
sampled in frequency space) to provide input to form "associations" 
with "High Level Memory" as given in the Association Filter Layer to 
Form a "String Section; 

c. ) a feedback loop is used to sequentially switch as described 
below to different "known", "standardized", "set", or "hardwired" 
delayed Gaussian filters which corresponds to changing the decay 
constant, a s , with a concomitant change in the half-width parameter, 
a s , and the number of elements, N s , with a concomitant change in the 



delay, 3—*-, where each a s and 3L_l is "known" or "standardized" from 
past "associations" ; 

d. ) "associations" are established between Fourier series by their 
"coupling" with Poissonian probability with Fourier series input from 
"High Level Memory" as given in the Association Filter Layer to Form a 
"String Section. The "associations" establish the relative position of the 
center of mass of each Fourier series such as a "group of SFCs" as 
either "before" or "after" the center of mass of the preceding and 
succeeding Fourier series "associated" with Fourier series input from 
"High Level Memory". The complete set of Fourier series "associated" 
with Fourier series input from "High Level Memory" covers all of the 
frequencies of the "string"; 

e. ) "groups of SFCs" of a sequential "set of subsets" are sequentially 
mapped to a time line by being added to the "before" or "after" end of 
the emerging temporally, conceptually, or causally ordered "string" 
wherein assignment of each "group of SFCs" is consistent with the 
frequency compositional and center of mass data to arrive at the order 
of the entire "string"; 

f. ) steps c)- e) are performed reiteratively until an order can be 
assigned without contradiction; 

g. ) the order of the associated "groups of SFCs" is established when 
"associations" with the "High Level Memory" are achieved for a given 
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set of delayed Gaussian filters (i.e. the order is established when 
internal consistence is achieved with input from ordered "strings" of 
High Level Memory); 

5 h.) the "groups of SFCs" of the " P string" of the form of Eqs. 

(37.113-37.115) that are parameterized according to their relative 
order are recorded to the "High Level Memory". For example, each 
Fourier series of the ordered string is recorded to the "High Level 
Memory" wherein each Fourier series of the ordered "string" is 
1 0 multiplied by the Fourier transform of the delayed Gaussian filter 

represented by e^^K^'^^K^'^ that established 
the correct order to form the ordered "string" represented by 



1 5 



f 2k 
k p ~n 



(37.115) 



Output of the Ordered " String" to High Level Memory Layer with 
Formation of the Predominant Config uration 

The activation of a "P element" increases its excitability or 
2 0 probability of future activation with input. Each "P element" has an 

activation memory with a finite half-life. Repetitive activation of a "P 
element" results in a longer half-life of the increased excitability; thus, 
the activation memory becomes long term. The same principle applies to 
cascade of association "stages" ("association ensembles") and M elements 

2 5 ("memory ensembles") and "configurations" of "couplings" of ensembles. 

For example, each "association ensemble" is comprised of "stages" in 
different states of "activity" where each state is equivalent to a 
microstate of statistical thermodynamics. A predominant configuration 
arises for any "association ensemble". Of the immense total number of 

3 0 microstates that can be assumed by an "association ensemble", an 

overwhelming proportion arises from one comparatively, small set of 
configurations centered on, and only minutely different from, the 
predominant configuration— with which they share an empirically 
identical set of macroscopic properties. On a higher level, a configuration 
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of "couplings" between "association ensembles" increases the activation of 
the "stages" comprising the "association ensembles". Analogously to 
statistical thermodynamics, a predominant configuration arises from the 
"association ensemble" level. Consider the "processor" on a higher level. 
5 The activation history of each "association ensemble" relates to a 

hierarchical activation relationship of coupled "association ensembles" 
which gives rise to a precedence of higher order predominant 
configurations. The ability to associate information and create novel 
information, is a consequence. Machine learning arises by the feedback 
1 0 loop of transducer input to the coupled predominant configurations which 
increases the basis for creating information with novel conceptual 
content. 



A summary of the method of Output of the Ordered "String" to High 
Level Memory Layer with Formation of the Predominant Configuration 
follows: 



a. ) the "groups of SFCs" of the " P string" of the form of Eqs. 
(37.113-37.115) that are parameterized according to their relative 

2 0 order are recorded to the "High Level Memory"; 

b. ) a counter corresponding to each "P string" and each "association 
ensemble" increases its stored count each time the "P string" or 
"association ensemble" is activated. In one embodiment, the count is 
also proportional to the length of time the "P string" or "association 

2 5 ensemble" is "active", and the count decays over time; 

c. ) the count is transformed into an expectation value and stored in 
a probability register which corresponds to each "string" and each 
"association ensemble"; 

d. ) during the process of establishing "associations" a probability 

3 0 operand causes a given "P string" or "association ensemble" to become 

"active" with an expectation value according to the value stored in its 
corresponding probability register; 

f.) on a lower level, the mechanism whereby past activation 
increases the probability of future activation applies to "P and M 
3 5 elements" as well; 

e. ) as more "P strings" are created, more "P elements", "M elements", 
and "stages" are activated, and more "association ensembles" are 
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created and activated, the relationship of the probability of future 
activation based on past activation gives rise to a processing 
predominant configuration of the "processor" analogous to that of 
statistical thermodynamics. 
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SUB-APPENDIX I 

The Input and the Band-Pass Filter of the Analog Fourier Processor 

The "P element" "impulse response" is a traveling wave in one 

5 spatial dimension (p) plus time (/ = -) where the wave function is a 

v 

dipole traveling at a constant velocity v. The magnitude of the potential, 
V, in cylindrical spacetime coordinates at the point (p,z) due to an 
"impulse response" centered at the position (p 0 ,z 0 ) is 

( 2 (z-Zo) 2 -(p-Po) 2 ) 

10 V = ^2 + ~/ ] 1 ®Hp-Po,z-z 0 ) (37.23) 

where 

z ° = v 'o (37.24) 
The potential is the convolution of the system function, h(p, z ), (the left- 
handed part of Eq. (37.23)) with the delta function (the right-hand part 

1 5 of Eq. (37.23)) at the position (p 0 ,Zo). A very important theorem of 

Fourier analysis states that the Fourier transform of a convolution is the 
product of the individual Fourier transforms, and the Fourier transform 
of a product is the convolution of the individual Fourier transforms [5]. 
The Fourier transform of the system function, h(p,z), is given in Box 16.1 

2 0 of the Superconductivity Section of Mills [6]. Also, see Mills [7]. 

An "impulse response" has the system function, h(p,z), which has 
the Fourier transform, H[k p ,k z ], which is shown in FIGURE 7. 

The output of a "P element", V, r , to an input of a pulse train of one 

2 5 or more "impulse responses" is another pulse train of "impulse 

responses". The spacetime "P element response", a pulse train function, is 
the convolution of the array pattern with the elemental pattern. The 
elemental pattern is the system function, h(p,z),-the spacetime potential 
function of an "impulse response". And, the array pattern is a finite 

3 0 periodic array of delta functions each at the center position of an 

"impulse response". 
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V, r (p,z(t)) = 



lp 2 +z 2 ] 



2 n 5/2 



' X a n 5 {P ~ n P 0 > Z - "V/ 0 )X 



P 2 ' Z 2 

(37.26) 



where a„ is a constant and C/fp + ^o,7 + il^ 

V 2 2 



Nvt n 



J is the unitary step function 
at P = zfl andz ^ = Z^ and u( p _^ z _N^ is tne unitary step 

function at position p = ^ a ndz = ^ = ^. Multiple "P elements" input 

into any given "P element" which then outputs to multiple "P elements". 
And, the amplitude, frequency, and length of the "P element response" 
(pulse train) is proportional to the length and rate of voltage change-the 
amplitude and rate of change of the input. Thus, each "P element" is an 
linear differentiator-the output, V 0ut , is the sum (superposition) of the 
derivative of the inputs. An exemplary output signal of an analog "P 
element" to an input of the form given by Eq. (37.26) is 
V 0u ,{p,z{t)) 

(2* 2 -p 2 ) 
[p 2 + z 2 ] 5 ' 2 



8p8z 



) X a » 5 (P~"Po>z- m *o) 



z + - 



-U\p 



Np 0 



,z-- 



(37.27) 



The Fourier Transform of the periodic array of delta functions of Eq. 
(37.27) is also a periodic array of delta functions in k,co -space 
1 



PoZont 



v of, 2k . 2k^ 
Z, a -° k o- n —>K ~n 



v 



(37.28) 



where z 0 = w 0 . The Fourier Transform of the window function given by 
the difference of the unitary step functions of Eq. (37.27) is the product 
of two sine functions in k,a> -space 



P 2 1 2 



(37.29) 



By the Fourier Theorem, the Fourier Transform of Eq. (37.26) is the 
product of the Fourier Transform of the elemental function, system 
function given by Eq. (37.25), and the Fourier Transform of the array 
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function given by Eq. (37.28) convolved with the Fourier transform of the 
window function given by Eq. (37.29). 



An 1 



1 +^z_ Po^O n 

kl 



7 2K 1 

k -n ,k -n 



2n 



. . Np 0 . . Nz 0 
sin/c„— ^ sin/< — - 
"2 z 2 



Po * vt 0 ) k p k z 



(37.30) 



Each "P element" is an linear differentiator— the output is the sum 
(superposition) of the derivative of the inputs. The differentiation 
property of Fourier transforms [8] is 



x(t)=]x(fy^df 



Differentiation 



(37.31) 

-j^ j2xfX(f) 

From Eqs. (37.30) and (37.31), the Fourier transform of a "P element 
response", v{k p ,k z ), called a "FC" is 

V{k p ,k z ) = k p k z p-— X a n$\ k p -n~,k z -n— <g>4 — 2 —Z- 

l + ^.PoZo.=. V Po vt 0 ) k k z 



k 2 2, k p -n—,k z -n~\®A & mk p ^mk z -f- 



(37.32) 



4tv 4 



sin 



k p -n 



2k 

ft. 



Np 0 



\ ( 
sin 



U vt °) 2 J 



Information "carried" by "P elements" may be represented by a Fourier 
series called a "SFCs" (series of Fourier components) comprising the 
superposition of the "P element responses" of multiple "P elements". Each 
"P element" contributes a Fourier component comprising an amplitude, 

No No 

15 a 0jn , at a specific frequency, m ™ m , "'^ , which is repeated 



as a series 



with a specific phase, 



nN 



A "SFCs" comprising the Fourier transform of 



the superposition of the "P element responses" of M "P elements", V y , is 



represented by 
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a n sin 
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2k 



N n p oPQ. 



sin 



k—n 



2k 



N „ lz0 Zo m 



o J 



(37.33) 

Each "FC M of Eqs. (37.33) is a series of a Fourier component with quantized 
frequency and phase angle. 
5 Consider the case that the amplitude of all "P element responses" , 

are equal where each amplitude is represented by a 0m . The "P element 

response" function given by Eq. (37.33) corresponds to recording to 
memory ("writing"). Consider the case that memory elements are 
activated to read the stored information. In one embodiment, this "read" 
1 0 operation is effected by a voltage ramp that is linear with time. The 
Fourier transform of the response is given by the differentiation and 
duality properties of Fourier transforms [8]. The "read" total response 
in Fourier space comprising a "SFCs", the superposition of M "FCs" 

HI 

wherein each "FC" corresponds to the response of a "M or P element" is 
V (k p ,k z ) 



-XX 



1 5 



47C 4 N m n Po m N mZ0 m . 

m — rr a 0m — ffl — - — -sin 

"i n^-<x> | _|_ zz- Po m z o m m 2 2 



k p -n 



2k 



Po, 



N ^P 



'« J 



(f 



sin 



k z ~n 



m 0 W J 



4k 



TT^n N m AT, sin 



k p ~n 



2k 
PoJ 



N m p 0 



(f 



sin 



k—n 



2k 



20 



(37.33a) 



Each "FC" of Eqs. (37.33a) is a series of a Fourier component with 
quantized amplitude, frequency, and phase angle. 

The relationship between k,co -space and real space is 



^ _ 2^ 2k _ 2k 
P ^ P P Wo 
k _2n __2k _ 2k 



(37.34) 



nvt n 



In k,co -space, the Fourier transform of the "impulse response" function 
(the left-hand side of Eq. (37.33)) filters each "FC" of a "P element". In 
the special case that 

(37.35) 
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the Fourier Transform of the system function (the left-hand side of Eq. 
(37.33)) is given by 

H = Ak (37.36) 
Thus, the Fourier Transform of the system function band-passes the 
5 Fourier Transform of the time dependent "P element response" function 
when the spatial frequency of the "FC M is equal to the temporal 
frequency. In one embodiment, "FC" filtering may be provided by 
adjusting the "P element" response corresponding to k p versus k z such 

that the band-pass condition of Eq. (37.35) is not met. In an analog 
1 0 embodiment, the "FC" may be filtered by adjusting the "impulse 
response" frequency as a function of time and therefore space 
corresponding to k p since the "impulse response" is a traveling wave. In 

another analog embodiment, the "FC" may be filtered by adjusting the 
conduction velocity which alters the output corresponding to k z , 

1 5 When the band-pass condition is met (Eq. (3735)), the Fourier 

transform of the superposition of a series of pulse trains of "impulse 
responses" of multiple "P elements" representing information is a series 
of trigonometric functions. Thus, in one embodiment of the present 
invention, the "processor" is an analog Fourier processor. According to 

2 0 the Fourier theorem any waveform can be recreated by an infinite series 

of trigonometric functions. 

*(*) = a 0 + X a n COS MS + Y, b n SiUC0 J (37.37) 
n-l n=l 

where a 0 , a n , and b n are constants. And, any aspect of the universe can be 
represented by an infinite series of sine and cosine functions as 

2 5 processed by the "processor". For the present "processor", the 

trigonometric function is the basis element of information. And, the 
complexity or information content of any analog waveform or digital 
equivalent is reducible to the number of Fourier components required for 
its assimilation. 

3 0 A unique feature of the present invention is that information is 

encoded in a Fourier series in k, co -space versus a conventional Fourier 
series in time and space. 
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SUB-APPENDIX II 

Modulation and Sampling Gives the Input to the Association Mechanism 

and Basis of Reasoning 

5 Each "P element" connects to multiple other n P elements" which 

further connect to association "stages" that propagate the "P element 
responses" as input along these "stages" in a linear cascade. Consider an 
amplifier made up of cascaded stages. The stages need not be identical. 
Let h^t) be the impulse response of the i th stage and assume that /*.(/)> 0, 
1 0 so that the step response of each stage (or indeed of any number of 

cascaded stages) is monotonic. Assuming that both integrals exist, T,., the 
normalized first moment of h t {t) is defined as 

J'MO* 

T g ==S (37.38) 

which can be interpreted as the center of gravity of a mass distributed 

1 5 along the t-axis with density h^t). If h^t) is positive, it is analogous to a 

probability density function, and T. corresponds to the statistical analog— 
the mean of A.(0* Thus, T i is considered as the measure of the delay in 
the impulse or step response of the /'* stage. The delay resulting from a 
cascade of n stages is the sum of the delays of each stage [9]; that is if 

2 0 MO = MO® MO®"-® MO (37.39) 

where <g) is the convolution operator, then 

T=T X +T 2 + — +T n (37.40) 

Similarly, assuming that both integrals exist, > the normalized 

moment of inertia about a center of gravity of a mass distribution h.(t) is 
2 5 defined as 



(AT,.) 2 =4 



(37.41) 

\{t-T,)\{t)dt 
= 4^_ 

\h.xt)dt 
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If h,{t) is positive, it is analogous to a probability density function, and 
can be interpreted as the statistical analog-the variance or 

dispersion of h,(t). AT; is twice the radius of gyration of the mass 
distribution. Thus, AT) is a measure of the duration of A.(f) or of the rise 
time of the step response of the /'* stage. The rise time resulting from a 
cascade of n stages is the sum of the rise times of each stage [10]; that is 
if h(t) is given by Eq. (37.39), then 

(AT) 2 =(AT 1 ) 2 +(A7 2 ) 2 +- + (ATJ 2 (37.42) 
Thus, in particular, for identical stages, the rise time is proportional to the 
square root of the number of stages. If h;(t) is not positive, rather than 
the definition of Eq. (37.41), the measure of duration is better defined as 



(AT) 2 



]t 2 h\t)dt 
]h\t)dt 



jth 2 {t)dt 
]h\t)dt 



= 4- 



(37.43) 



)hMdt 



In many ways AT of Eq. (37.43) is the most analytically satisfactory 
simple general measure of duration; for virtually any h t (t) for which the 
integrals exist, Eq. (37.43) will give a reasonable estimate of duration. 
Equivalently, possibly the best simple measure of bandwidth for real 
lowpass waveforms is 

]f 2 Hffd f 



(AW) 



jHffdf 



(37.44) 



From the definitions of AT and AW given by Eq. (37.43) and Eq. (37.44), 
respectively, it is possible to prove the following Uncertainty Principle 
[9]: 



25 



For any real waveform for which A7 and AW of Eq. (37.41) and Eq. 
(37.43) exist, 



A7AW> 



1 



K 



(37.45) 
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In other words, AT and AW' cannot simultaneously be arbitrarily small: A 
short duration implies a large bandwidth, and a small-bandwidth 
waveform must last for a long time. 

Consider a cascade of association "stages". The Uncertainty 
Principle given by Eq. (37.45) applies to the "P element response" as it is 
transmitted from one "stage" to another in the cascade. In one 
embodiment, the "voltage" decays exponentially at the junction or linkage 
of any two "stages". The cascade forms a filter, and an ideal filter 
response is that which has the smallest duration-bandwidth product in 
the sense of Eqs. (37.43) and (37.44). Such a response is a Gaussian pulse 
which also has the same form in the time and space domain [9]. 
However, a Gaussian pulse cannot be the impulse response of any casual 
system, even with substantial delay. Consider, for example, an //-stage 
amplifier with the impulse response of each stage equal to 

h(t) = ajNe- am u{t) (37.46) 
The frequency response of the cascade of N such stages is 

/ "\ w 
1 



//„(/) = ["(/)]" = 



1 + 



J2ttf 



(37.47) 



The shape of H N (f) for large N can be determined by taking logarithms 
and using the power series expansion 



ln(l + x) = x- — + — ... 

2 3 

The power series expansion of the In H N is 

In H N — —N 



(37.48) 



(37.49) 



a 2\ a ) 

where the remaining terms vanish as fast as for lar S e N - Thus > th e 

frequency response tends to 

H N (f)~e*«)e t«J (37.50) 
for large N, and the impulse response of the cascade tends to 
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h »W~J^ e ai (37.51) 



that is, a Gaussian pulse delayed by ^E-. This result is a very special case 

of a remarkable theorem [10]-the Central Limit Theorem of probability 
theory-which states in effect that, under very general conditions, the 
5 cascade of a large number of linear-time-invariant (LTI) systems will 
tend to have a Gaussian impulse response, almost independent of the 
characteristics of the systems cascaded. Sufficient conditions of the 
Central Limit Theorem are that 

10 1- The absolute third moments, 

„ 3 

Jl'KM* (37.52) 
exist for all components of the systems and are uniformly bounded; 

2. The durations, AT,, of the component systems in the sense of Eq. 

1 5 (37.43) satisfy the relation 

JS^X^) 8 * 0 (37.53) 

For large N, the first condition allows the higher order terms in the 
expansion such as Eq. (37.49) to be ignored, and the second condition 

2 0 guarantees that no finite subset of the component systems will dominate 

the result because the remainder have relatively wide bandwidths. 
Given this theorem, it follows from Eqs. (37.38-37.45) [9] that the overall 
impulse response of N cascaded stages is approximately 

k 

h{t)~^L—e 2(Ar > (37 54) 

2 5 where T and AT are given by Eq. (37.40) and Eq. (37.42), respectively, 
and 

N ~ 

k = H\ h Xt)dt (37.55) 

Eq. (37.54) is a filter function. Consider Eq. (37.33) where the 
Fourier transform of the superposition of "P element responses" (the sum 
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of multiple pulse trains of "impulse responses" representing information) 
is given by a series of trigonometric functions wherein the "processor" 
can be an analog Fourier processor. The input of information to the 
association mechanism arises as the Fourier series is modulated and 
5 sampled. Consider the output of a cascade of association "stages"--each 
with an "impulse response". The "stages" are cascaded as an //-stage 
amplifier with the transmission impulse response of each stage in one 
embodiment equal to that of a decaying exponential given by Eq. (37.46). 
The filtered signal is the sum of the convolution of the response of each 

1 0 transmission stage of the cascade with each "P element response". Using 

the distributive, commutative, and associative laws of the convolution 
operation and using the Central Limit Theorem, the filtered signal is the 
convolution of the superposition of the "P element responses" over the 
cascade of "stages" also given by Eq. (37.27) with the Gaussian response 
15 given by Eq. (37.51). A very important theorem of Fourier analysis 

states that the Fourier transform of a convolution is the product of the 
individual Fourier transforms [5]. Thus, the output of a cascade of N 
stages each with a transmission decay constant a (corresponding to the 
transmission impulse response) is the product of Eq. (37.33) and Eq. 

2 0 (37.50). By changing the decay constant a and the number of "stages" N 

in the cascade, Fourier series representing information including that 
from memory can be filtered and delayed (modulated and sampled in 
frequency space) to provide input to form associations of the Association 
Filter Layer. For example, consider the result on exemplary filter 

2 5 functions and the corresponding Fourier transforms shown in FIGURES 8 

to 15 as the decay constant a and the number of "stages" N of each 
corresponding cascade are altered. In frequency space, the time delayed 
Gaussian filter corresponds to modulation and sampling of the Fourier 
series representation of the memory output comprising the superposition 

3 0 of multiple "M element responses". Thus, the time delayed Gaussian filter 

selects memory output and provides input for the association mechanism 
and basis of reasoning. 

In another embodiment, the time delayed Gaussian filter may be 
modulated in the time domain to effect a frequency shift in k, co -space. 
3 5 The shift follows from Eq. (37.109) and the duality property of Fourier 
transforms [8]. 
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SUB-APPENDIX III 
The Association Mechanism and Basis of Reasoning 

"Coupling" 

5 A cascade of association "stages'* called an "association ensemble" is 

"activated" with input from "M elements", "P elements", or "stages" of a 
different "association ensemble". The "association ensemble" is "active" if 
it is "carrying" a Fourier series such as a "SFCs" wherein "active" in the 
digital case may refer to a recall of an "SFCs" from memory followed by 
1 0 steps a-i of the Association Filter Layer to Form a "String" Section. The 
"association ensemble" is "inactive" if it has no output and is not 
"carrying" a Fourier series such as a "SFCs" wherein "inactive" in the 
digital case may refer to no recall of an "SFCs" from memory. 

In an analog embodiment, the "stages" of an "association ensemble" 

1 5 are intraconnected and interconnected. A first "active" cascade of 

association "stages" can interfere with and "couple" with a second set, 
third set, etc. The probability distribution function of "coupling" between 
a first "active" "association ensemble" and at least one other "active" 
"association ensemble" is Poissonian. Each "association ensemble" is 

2 0 comprised of a large number of cascaded association "stages" each weakly 

linked to one or more "stages" of the one or more different "association 
ensembles". (The "coupling" is analogous to interference between 

coherent or harmonic states.) The probability P T 



first "active" cascade of association "stages" with modulation e' '^'^ 

2 5 given by Eq. (37.50) will interfere with and "couple" with s separate 

"active" cascades of association "stages" ("association ensembles") each 

with modulation e v<x > } given by Eq. (37.50) can be derived from the 
correlation function (Eq. (37.78) for the statistical average of the large 
number of possible "couplings" between the individual weakly linked 

3 0 "stages". 

The physical behavior of a large number of "active" cascaded 
association "stages" (an "association ensemble") each weakly linked to 
provide a Poissonian probability of "coupling" to one or more "stages" of 
one or more different "association ensembles" is equivalent to that of the 
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interaction of ultrasound with Mossbauer gamma rays. Each "association 
ensemble" "carries" a Fourier series in k,ca -space such as a "M or P 
element response" which comprises a sum of harmonic functions. Thus, 
physically, the former case corresponds to interference of a first Fourier 
5 series input filtered by an "association ensemble" with a second, third, or 
s th Fourier series input filtered by s th "association ensemble". The 
latter case corresponds to interference of an electronic transition and an 
oscillator transition. In both cases, a harmonic energized state interferes 
with another. 

1 0 Consider the Lamb-Mossbauer formula for the absorption of a 7 

ray of energy E by a nucleus in a crystal given by Maradudin [11]. 

{E Q -E + E n -E m ) 2 +^ 

In this equation, E 0 is the energy difference between the final and initial 

nuclear states of the absorbing nucleus, E m and E n are the energies of the 

1 5 eigenstates \m) and \n) of the crystal, respectively, T is the natural width 

of the excited state of the nucleus, p is the momentum of the 7 ray, R(/) 
is the instantaneous position vector of the absorbing nucleus, Z is the 
crystal's partition function, T = (kpy\ and cr 0 is the resonance absorption 
cross section for the absorbing nucleus. By expressing the denominator 

2 0 of Eq. (37.56) as an integral, Eq. (37.56) is equivalent to 

- \ cr 0 r]dte iox ~^x(cxp[~ik * u(/;0]exp[/k • u(/;0)]) (37.57) 

wherein the position vector R(Z) is 

R(/) = x(/) + u(Z) (37.58) 
For, Eq. (37.58), x(/) is the position vector of the mean position of the 

2 5 absorbing nucleus, and u(/) is its displacement from the mean position. 

Eq. (37.57) follows from Eq. (37.56) with the following substitutions: 

Q~) p = k (37.59) 

hco=E-E 0 (37.60) 

r~ (37.61) 

3 0 and u(/;/) denotes the Heisenberg operator, 

u(l;t) = e {h) u(/;0)e UJ (37.62) 
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where H is the Hamiltonian. The angular brackets in Eq. (37.57) denote 
an average over the canonical ensemble of the crystal. 

The correlation function for the statistical average of a large 
number of "active" cascaded association "stages" (an "association 
5 ensemble") each weakly coupled to one or more "stages" of one or more 
different "active" "association ensembles" is equivalent to that of the 
interaction of ultrasound with Mossbauer gamma rays. From Eq. (37.57), 
the correlation function Q{t) of acoustically modulated gamma ray 
absorption by Mossbauer nuclei is 
1 0 G(0 = (exp[-/k*u(/;0]exp[/k.u(/;O)]) (37.63) 

In the present case, u(/) corresponds to the delay of an "association 
ensemble" s comprising a time delayed Gaussian filter. In Jk, co -space, the 
time delay corresponds to a modulation of the s th Fourier series (e.g. "P 
or M element response" given by Eq. (37.33)) that is "carried" by the 

1 5 "association ensemble" s). Since the Fourier series is a sum of 

trigonometric functions in k,co -space, the modulation corresponds to a 
frequency shift of the Fourier series "carried" by the "association 
ensemble" 5. k of Eq. (37.59) corresponds to the wavenumber of the 

TP 77 

frequency shifted s th Fourier series. ~ 0 of Eq. (37.60) is the shifted 

2 0 frequency of a first Fourier series that is "carried" by a first "association 

ensemble". 

In the case of acoustically modulated gamma ray absorption by 
Mossbauer nuclei, u(/;r) of Eq. (37.62) is 



2 5 The matrix elements of Eq. (37.63) are calculated by using the theorem 




u(/;0)e 




(37.64) 




(37.65) 



For a harmonic oscillator, the commutator of k-u(/;/) and k-u(/;0) is a c 
number; thus, 



30 



Q{t) = (exp[-/k ■ u(/;f)]exp[zk • u(/;0)]) 

= (exp[-/k-[u(/;0-u(/;0)]])Xexp |{[k-u(/;/), k.u(/;0)]) 



(37.66) 



Since the correlation function applies to an ensemble of harmonic 
oscillator states, the first thermodynamic average can be simplified as 
follows: 
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(exp[-ik • [u(/;0 - u(/;0)]]) = exp[-±({k • [u(/;/) - u(/;0)]} 2 )] (37.67) 

This theorem is known in lattice dynamics as Ott's theorem [13] or 
sometimes as B loch's theorem [14]. Using the time independence of the 
harmonic potential, Eq. (37.67) is 

exp[-I({k-[u(/;0-u(/;O)]} 2 }] = exp-[i([k-u(/;0] 2 ) + |([k-u(/;O)f)l (37.68) 

= exp-([k-u(0] 2 } (37.69) 
Substitution of Eqs. (37.67-37.69) into Eq. (37.66) gives 



e(0 = exp(-[k-u(/;/)] 2 )xexp i([k-u(/;/), k-u(/;0)]) 



(37.70) 



Expanding u a (/;f) in terms of the normal coordinates of the harmonic 
1 0 potential and the phonon operators of that harmonic potential gives 



(37.71) 



where a labels the Cartesian components, M l is the mass of the ion in the 
/th experiment, co s is the frequency of the sth normal mode, B {s) (l) is the 
associated unit eigenvector, and foj and b s are the phonon creation and 
1 5 destruction operators for the sth normal mode. By use of the coordinate 
expansion, the exponential of the correlation function appearing in Eq. 
(37.70) can be written as 



(k-u(/;/)k-u(/;0)) _ 



n 



(37.72) 



J 0 {2c]) + ^J n {2c]) 



T+(r.)' 



where the following substitutions were made: 



20 



. 1 hc °s 

7s=- = ** r 



n r = 



tt<a s 
,~kT~ 



(37.73) 
(37.74) 



WO 99/34322 



PCT/US98/27624 



69 



1 0 



1 5 
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CO, 



2M, '« ^ 



(37.75) 



kT 



e" -1 

and where the Bessel function relationship [15] 



— xi v+v ) , 



(37.76) 

was used. n s is the mean number of phonons in the slh mode at 
temperature T. 

In the case of "coupling" between a first "active" "association 
ensemble" and at least one other "active" "association ensemble", the 
correlation function is independent of time-not a function of e /GV and 
e -iav Thus, the time dependent factors are dropped in Eq. (37.72), and 
combining Eqs. (37.70-37.72) and Eq. (37.75) gives the correlation 
function as 

Q{c?) = txp-clYlJ 0 (2c 2 M ) (37.77) 

s 

For the "coupling" of "active" "association ensembles", the partition 
function of Eq. (37.56) is equal to one. By the Central Limit Theorem, s=l 
in Eq. (37.72) corresponds to each cascade of association "stages" giving 
rise to a specific frequency shift. The correlation function for each 
"association ensemble" is 

Q(c 2 ) = exp- [c 2 }j 0 (2c 2 ) (37.78) 



The probability P T 



that a first "active" "association 



ensemble" will "couple" with s "active" "association ensembles" can be 
derived from the correlation function, Eq. (37.78). The expansion of the 
Bessel function is 



f -x 2 ^ 



7vW "V2j5[m!r(m + v + l)] 



^x 2 ^ 



(37.79) 



y ° W = 5[m!r(m + l)]^5 



[m!m!] 

where r(m + l)-m! was used. The probability distribution function of 
"coupling" between "association ensembles" (coherent states) is 
Poissonian. From SUB-APPENDIX II-Modulation and Sampling Gives the 
Input to the Association Mechanism and Basis of Reasoning, the output of 
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a cascade of TV "stages" with a transmission decay constant a is the 
product of Eq. (37.33) and Eq. (37.50). From Eqs. (37.40), (37.42), 
(37.46), (37.50), (37.51), and (37.54), the frequency shift of "association 



ensemble" one is 2LL-L 9 and the frequency shift of "association ensemble" s 



5 is - L ^- where the impulse response of each "stage" in both "association 
ensembles" is 

h(t) = a4Ne- a ^u{t) (37.80) 
"Coupling" of filtered Fourier series is based on their spectral 
similarity. In one embodiment, the "coupling cross section" amplitude, /J s 2 , 
1 0 is given by the integral of the product of the spectrum of the first Fourier 
series sampled and modulated by the first "association ensemble" and the 
complex conjugate of the spectrum of the s th Fourier series sampled and 
modulated by the s th "association ensemble". The spectrum of a Fourier 
series ("SFCs") sampled and modulated by an "association ensemble" is 
1 5 given by the product of Eq. (37.33) and Eq. (37.50). Thus, Eq. (37.75) is 

# = JJ (37.81) 
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Let k p = k z , then p 0 = z 0 = vt 0 . Thus, Eq. (37.82) is 



(37.82) 
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(37.83) 
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Substitution of k z =—+- and sinG = e~ J0 into Eq. (37.83) gives 



4 



+7 



(37.84) 



The integral o f Eq. (37 .84) is given by Hogg and Tanis [16] 
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where a 2 = — ^-r- and / = - j 
cc 2 + cc 2 s 



2v„ 



2v„ 



(37.85) 



in corresponding 



integrals. It was given previously (Eq. (37.83)) that p 0 =^ 0 = v / 0 ; thus, Eq. 
1 0 (37.85) simplifies to 
1 



a 2 a 2 A 4 



^ia 2 + a 2 ^ t p 2 ^ 



Z X \ cos7V ( nN m, ~ «' JV m , )|exp- 



2 "0. 

Po„. " 



2 2 

a, a; 



a? + a? 



(37.86) 



where cr 2 = 



2 2 

a? + a 2 ; 



and r = 



■ + 



in corresponding 



1 5 integrals. Consider the case that the amplitude of all "P element 

responses" are equal, thus = a K for all m x and m s in Eq. (37.86). In the 

case that each "SFCs" is represented by Eq. (37.33a), Eq. (37.86) is 
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(37.87a) 
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where a 2 



2„2 



and r = -y 



a:, a, 2 2 



in corresponding 



integrals. 

In one embodiment, the present "processor" is an analog Fourier 
processor wherein the data is digitized according to the parameter N of 
Eqs. (3733), (37.33a), and (37.87). Each "FC" of Eqs. (37.33) is a series of 
a Fourier component with quantized frequency and phase angle. Each 
"FC" of Eqs. (37.33a) is a series of a Fourier component with quantized 
amplitude, frequency, and phase angle. Each "SFCs" represented by Eq. 
(37.33) and Eq. (37.33a) is filtered and delayed in the time domain 
(modulated and sampled in the frequency domain) as it is recalled from 
memory and processed by an "association ensemble". "Association 
ensembles" produce interference or "coupling" of the "SFCs" of one set of 
"M or P elements" with that of another by producing frequency matched 
and phase locked Fourier series —sums of trigonometric waves that are 
frequency matched and periodically in phase— that give rise to 
"association" of the corresponding recalled or processed information. The 
Poissonian probability of such "coupling" (Eq. (37.106)) can be derived 
from the correlation function (Eq. (37.78) wherein Eq. (37.87) is a 
parameter. The magnitude of the "coupling cross section" of Eq. (37.87a) 
and Eq. (37.86) is independent of any phase matching condition because 
the phase angle is quantized. Thus, the argument of the cosine function 
of Eq. (37.87a) and Eq. (37.86) is zero or an integer multiple of n. 
Consequentially, in each case, the corresponding time convolution of Eq. 
(37.84) is a cyclic convolution, and the sum over n is eliminated. 
Whereas, the frequency matching condition provided by the frequency 
shifts of the cascades of association "stages" comprises the zero argument 
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of the exponential function of Eq. (37.87a). Thus, the magnitude of the 
"coupling cros s section " follows from Eq. (37.87a) 

#=(87T) 2 1 
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(37.87b) 



In terms of the relationship p Q = Zo = vto , Eq. (37.87b) is 
#=(8*) 2 



2 2 



ni, =1 



2 2 



a, 2 + a 2 



a, a 2v . 2v 



(37.87c) 

"Coupling" between "active" "association ensembles" further 
1 0 depends on the frequency difference angle, ^, between the two or more 
Fourier series "carried" by the corresponding "association ensembles". In 
k,co- space, the information is represented as Fourier series which 
comprise sums of harmonic functions. Thus, the "coupling cross section" 
is a complex number with a projection in k,co- space that is a function of 



1 5 the frequency shift 



a, 



of the first "association ensemble" and the 



frequency shift ^-jl G f t h e s th "association ensemble". The frequency 

shift of each "association ensemble" corresponds to the respective 
modulation function given by the Fourier transform of the delayed 
Gaussian filter (Eq. (37.50)). The resultant "coupling cross section", 
20 (& 2 (0j)> as a function of frequency difference angle, 0 S , is given by 

(m<Ps)) = fte i2 *> ? (37.88) 

where the frequency difference angle, 0 s , is 
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(37.89) 



Thus, the "coupling cross section" given by Eq. (37.88) is a dimensionless 
complex number that comprises a "coupling cross section" amplitude, & 2 , 
and frequency difference angle, <p s , of the harmonic "coupling". In other 
embodiments of the present invention, further operations may be 
performed on such as phase shifting, normalizing to a given 

parameter, scaling, multiplication by a factor or parameter such as a gain 
factor, addition or subtraction of a given parameter or number such as an 
offset, etc. In a further embodiment, may be represented by 

different equations than those such as Eq.(37.87c) and Eq. (37.81) that 
also represent the spectral similarity and difference of the frequencies of 
filtered or unfiltered Fourier series that may "couple". 

In the case that p 0 =z 0 = vt 0 , the frequency difference, is 



K 



a. 



+ 



M - N L 



•X 

m,=l 



a, ~, 2 



(37.90) 



The probability distribution of "coupling" between two "association 
ensembles" each "carrying" a Fourier series such as a "SFCs" is Poissonian 
with mean number of "stage" "couplings" 

<"> = £ 2 (37.91) 
The probability [17] of "coupling" with a second "association ensemble" 
with m "couplings" between "stages" is 

p __(nTe-<«^{P 2 ) m e-e 2 fj^ 



ml 



ml 



ml 



with mean number of "stage" "coupling" events {n) = p 2 



(37.93) 
The probability 



can be derived by factoring Eq. (37.93) from the 
Bessel function of the correlation function (Eq. (37.78)) and its expansion 



2 5 which follows from Eq. (37.79). 
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(37.94) 



f3 2m e 



Combining Eq. (37.93) and Eq. (37.94) demonstrates that the probability 



= P r (J3x) is proportional to 



5 Let i = /, then the change of variable in Eq. (37.95) is 



/ 5 \ w/2 



m=o m - nil 



Let rri = m/2 9 then the change of variable in Eq. (37.96) is 

2 \m/2 ✓ 2 xm' 



"4 ) 



"4 y 



m =o ™! ^ (2m 1 )! 
The series expansion of cos(jc) is 

< ^ v M" 



(37.95) 



(37.96) 



(37.97) 



(37.98) 



Combining Eq. (37.78) and Eqs. (37.95-37.98) gives the probability 



a, ' a 



* J 



proportional to 



°= COS 



where ), = V* = V^F- From Eqs. (37.81-37.90), 
Combining Eq. (37.99) and Eq. (37.100) gives the probability 



(37.99) 



(37.100) 



proportional to 
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a. 



a 



COS 



(2sin0,) 



(37.101) 



where (p s is the frequency difference angle. Combining Eq. (37.78), Eq. 
(37.100), and Eq. (37.101) gives the probability P x (<p) proportional to 



: ex P[~A" 2 sin 2 &]cos(2sin (p 5 ) = exp -/?; 2 ^- 



-cos2^ 



cosi 



(2 sin 0 5 ) 



(37.102) 

where 0 S is the frequency difference angle and £ 2 is the "coupling cross 
section" amplitude. 

According to the time delay property of Fourier transforms [8], a 
time delay, S(t-t 0 ) 9 during independent activation of a given "association 

ensemble" with recall from memory is equivalent to a phase shift of the 
correlation function given by Eq. (37.63) 

2(0 = (exp/<S exp[-zk • u(/;/)]exp[/k • u(/;0)]) (37.103) 
Thus, Eq. (37.102) is phase shifted. 



Mr 



■ exp 



2 f l-cos2^ 



-jJcos(5 s 



+ 2sinfc) (37.104) 



where <p s is the frequency difference angle, & 2 is the "coupling cross 
section" amplitude, and S s is the phase shift. 

In an analog embodiment, each of the s separate "association 
ensembles" that may "couple " with the first "active" "association 
ensemble" may be "inactive" before "coupling". The "coupling" causes the 
corresponding "association ensemble" to become "active". Eq. (37.104) 
represents the probability that a first "active" "association ensemble" will 
"couple" with s "active" "association ensembles" as a function of the 
frequency difference angle, <j> s , the "coupling cross section" amplitude, ft 2 , 
and the phase shift, S s , Eq. (37.104) also represents the probability that a 
first "active" "association ensemble" will "couple" with and "activate" s 
"inactive" "association ensembles" as a function of the frequency 
difference angle, <p s , the "coupling cross section" amplitude, /J 2 , and the 
phase shift, S s . In the latter case, the Fourier series such as a "SFCs" 
"carried" by the "activated" s th "association ensemble" may be "linked" 
with the "association ensemble". The "linkage" is as described for 
"transducer strings" in SUB-APPENDIX Vl-Input Context. 
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"Association" 

Given that a first "association ensemble" is "active", the probability 
of the occurrence of either the "active" state or the "inactive" state of the 
s th "association ensemble" is one. In one embodiment, in the absence of 
interference (i.e. "coupling") between the "association ensembles", the 
probability of the occurrence of the "active" state of the s th "association 
ensemble" is the same as the probability of the occurrence of the 
"inactive" state- 1/2. However, in the event that "coupling" between the 
first and s th "association ensemble" may occur, the s th "association 
ensemble" may be "activated". The probability of the occurrence of the 
"active" state of the s th "association ensemble" with the possibility of 
"coupling" with the first "active" "association ensemble" is 1/2 plus the 
probability function, Eq. (37.104), normalized to 1/2. Therefore, given 
that the first "active " "association ensemble" may "couple" with the s th 
"association ensemble", the probability function for the occurrence of the 
"active" state of the s th "association ensemble" is 

l-cos2<ft/ 
2 



A 



1 + exp 



cos(<5 5 -i^sin^) 



(37.105) 



where <p s is the frequency difference angle, p 2 s is the "coupling cross 
section" amplitude, and 8 S is the phase shift. 

In an embodiment, the two Fourier series (e.g. each a "SFCs") are 
"associated" if they are "active" simultaneously. Thus, given that the first 
"active" "association ensemble" may "couple" with the s th "association 
ensemble", Eq. (37.105) is the probability function for the occurrence of 
the "association" of the Fourier series of the first "active" "association 
ensemble" with that which may be "carried" by the s th "association 
ensemble" as a function of the frequency difference angle, <j> s , the 
"coupling cross section" amplitude, /J 5 2 , and the phase shift, d s . 

The probability of the occurrence of "association" between a first 



Fourier series and s other Fourier series P A 



a 0 



A 



wherein 



the first "active" "association ensemble" may "couple" with each of s 
"association ensembles" is the product of the probabilities 

? _ 2 f l-cos20 5 



a. 



-A 



1 + exp 



r2 p-cos2^ - 



cos(5 J + 2sin^) 



(37.106a) 
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wherein the first "association ensemble" provides modulation e K<Xl } 
given by Eq. (37.50), the s th "association ensembles" provides 

modulation e K 1 } given by Eq. (37.50, (f> s is the frequency difference 
angle, (5 2 s is the "coupling cross section" amplitude, and 5 S is the phase 



5 shift. The plot of the probability P / 



— ,5_ of the occurrence 



of "association" of the first Fourier series with the s th Fourier series 
according to Eq. (37.106a) is given in FIGURES 16 A-C and FIGURES 17 A- 
D. 

In another embodiment, in the absence of "coupling" between the 
"association ensembles", the probability of the occurrence of "association" 
is p T . With the replacement of 1/2 of Eq.(37.106a) with /? T , the 

probability of the occurrence of "association" of the corresponding Fourier 
series based on a first "active" "association ensemble" "coupling" with s 
separate "association ensembles" is 

(37.106b) 

where p T is the probability of the occurrence of "association" in the 
absence of "coupling", <j) s is the frequency difference angle, $] is the 
"coupling cross section" amplitude, and 5 S is the phase shift. 

Eq. (37. 106b) gives one as the maximum probability of the 
occurrence of "association". In other embodiments, the probability 
maximum may be less than one. In this case, Eq. (37.106b) is 



n 



cos[S s +2sin0 5 ) 



(37.106c) 



where P is the maximum probability of the occurrence of "association". 
Eq. (37.105) and Eq. (37.106) represent the "association" probability 
parameter. 

The probability of "association" of Fourier series was herein derived 
for Poissonian statistics using delayed Gaussian filters; however, the 
invention is not limited to Poissonian statistics and the use of Gaussian 
filters. In other embodiments, the "association" can be based on 
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alternative statistics corresponding to their respective distributions. 
Examples are Gaussian or normal statistics, binomial statistics, Chi-square 
statistics, F statistics, and t statistics. Other statistical distributions are 
given in Hogg and Tanis [18] which are herein incorporated by reference. 
5 Furthermore, in other embodiments, the "association" can be base on 

alternative filters such as Butterworth, band pass, low pass, and high pass 
filters. Other filters are given in Siebert [19] which are herein 
incorporated by reference. 

In an analog embodiment, "coupling" may potentiate the two or 

1 0 more Fourier series. For example, each signal in the time domain 

corresponding to the Fourier series k,co -space may repeat in time and 
therefore increase in duration. In an embodiment, potentiated Fourier 
series are recorded to memory as "associated" with a probability that 
depends of the potentiation. In an embodiment, the probability is given 
15 by Eq. (37.106c) wherein the potentiation is via "coupling". 

Eq. (37.106c) also applies to the probability of "association" 
between a first "active" "association ensemble" and s "active" "association 
ensembles". In this case, an equivalent digital embodiment comprises the 
recall of Fourier series such as two or more "SFCs" from memory followed 

2 0 by steps a-i of the Association Filter Layer to Form a "String" Section. 



Eq. (37.106) gives the probability P A 



of the 



occurrence of "association" of the corresponding Fourier series based on a 

first "active" "association ensemble" with modulation e V(Xl) given by 
Eq. (37.50) "coupling" with s separate "association ensembles" each with 

2 5 modulation e K<X ' J given by Eq. (37.50) and independent phase shift, S s . 

The process of first establishing "associations" between different Fourier 
series representative of different pieces of information is the basis of 
producing information with novel conceptual content. The formation of 
"associations" is also the basis of reasoning. The generation of 

3 0 "associations" depends on the statistics of "coupling" of multiple 

"association ensembles" each comprised of cascaded association "stages". 
Then the "associated" information is ordered or further processed to 
provide general context such as cause and effect relationships by a 
mechanism involving the half-width parameters, a s , the time delay 
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parameters, and potentially the independent phase shifts, 8 S , of Eq. 

(37.106). The ordering of "associated" information is described in SUB- 
APPENDIX IV-Ordering of Associations: Matrix Method. 



5 
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SUB-APPENDIX IV 
Ordering of Associations: Matrix Method 

The set of "associated" Fourier series such as at least two "groups of 
5 SFCs" and/or at least two "SFCs" is herein called a "string". The "string" is 
a superposition of Fourier series; thus, it comprises a Fourier series, a 
linear sum of "FCs". FIGURE 19 is a flow diagram of an exemplary 
hierarchical relationship of the signals in Fourier space comprising "FCs", 
"SFCs", "groups of SFCs", and a "string" in accordance with the present 
1 0 invention. Each "FC" is the output of a "P element" or is stored into 
and/or recalled from a "M element" as shown in FIGURE 18. The 
information of "string" may be ordered to provide cause and effect, 
chronology, and hierarchical relationships. The ordered "string" is 
retained in memory to provide successive associative reference or 

1 5 further ordering of information. The information of the "string" is 

ordered or sequenced temporally, conceptually, or according to causality 
via the Matrix Method of Analysis of Mills [3, 4]. 

Consider Eqs. (37.33) and (37.33a) where each represents a "SFCs" 
in k, cq — space comprising a Fourier series. A "string" is a sum of Fourier 

2 0 series which follows from Eqs. (37.33) and (37.33a) as follows: 
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(37.108) 
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The corresponding equations in the time domain are a sum of multiple 
finite series of traveling dipoles ("impulse responses") wherein each 
dipole series is periodic in space and time. In frequency space, each time 
delayed Gaussian filter ("association ensemble" corresponding to a "SFCs") 
5 modulates and samples the Fourier series which encodes information. 
Thus, the time delayed Gaussian filter selects information from the 
"string" and provides input for the association mechanism as the 
"processor" implements the Matrix Method of Analysis to find the order 
of the associated pieces of information represented by each "SFCs" of the 
1 0 "string". 

Consider the time interval * = to t = t f of a "string" associated by 
"association ensembles" and recorded to memory. By processing the 
"string' with multiple "association ensembles" comprising Gaussian filters 



JN 

each of different delay, -*-^ L , and half- width parameter, a s , the "string" 
1 5 can be "broken" into "groups of SFCs" each having a center of mass at a 



A N 

time point corresponding to the delay - and frequency composition 

corresponding to a s which form a nested set of "sequential subsets" of 
"groups of SFCs" of the "string" in k,& -space. The set members map to 
time points which are randomly positioned along the time interval from 
2 0 the f = r r side and the r = r 2 -side as shown in FIGURES 8, 10, 12, and 14. 

This nested set of "sequential subsets" of random "groups of SFCs" 
mapping to random time points from the r = r r side and the r = r / -side is 

analogous to the nested set of "sequential subsets" of random DNA 
fragments from the 5' end and the 3' end. The order in both cases can be 

2 5 solved by the Genomic DNA Sequencing Method/Matrix Method of 

Analysis of Mills [3, 4] described in SUB-APPENDIX V. 

The output of an association filter is the convolution of the input 
"groups of SFCs" (each "SFCs" given by Eqs. (37.33) and (37.33a)) of a 
"string" (Eq. 37.108) or the "string" itself with a delayed Gaussian. In 

3 0 terms of the matrix method of analysis (hereafter "MMA"), the filter 

parameter a of the time delayed Gaussian filter corresponds to the 
acquisition of the composition of a polynucleotide member of a nested set 
of subsets. The time delay (time domain) and modulation (frequency 

[N 

domain) parameter J— determines the center of mass of the output, and 
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it corresponds to the terminal nucleotide data. By forming "associations" 
with input from "High Level Memory" as given in SUB-APPENDIX Ill- 
Association Mechanism and Basis of Reasoning, the "processor" 
determines the relative position of the center of mass of each Fourier 
5 series such as a "group of SFCs" as either "before" or "after" the center of 
mass of the preceding and succeeding Fourier series "associated" with 
Fourier series input from "High Level Memory". The complete set of 
Fourier series "associated" with Fourier series input from "High Level 
Memory" covers all of the frequencies of the "string". By Parseval's 
10 theorem, by processing the entire interval in k 9 co -space, the information 
is entirely processed in the time domain. The order such as temporal 
order of the Fourier series representing information is determined using 
the MMA. 

Input to form associations is provided by changing the decay 

1 5 constant a and the number of "stages" in the cascade N 9 or by processing 

each "group of SFCs" of a "string" using an "association ensemble" with 
different parameters a and N over all "groups of SFCs" that make up the 
entire "string". Each "group of SFCs" is determined to be on the / = /,-side 
or the r = reside of the "axis" of the "string" corresponding to the S'-side 

2 0 or 3*-side of the "axis" of a polynucleotide to be sequenced via the Matrix 

Method of Analysis. A feedback loop comprises sequentially switching to 
different "known", "set", or "hardwired" delayed Gaussian filters which 
corresponds to changing the decay constant, a s , with a concomitant 
change in the half-width parameter, a s9 and the number of elements, N s9 

2 5 with a concomitant change in the delay, where each a and is 

"known" from past experiences and associations. The feedback loop 
whereby information from memory encoded in the "string" is filtered and 
delayed (modulated and sampled in frequency space) to provide "FCs", 
"SFCs" or "groups of SFCs" which are associated with input from "High 

3 0 Level Memory" provides the data acquisition and processing equivalent 

to the formation, acquisition, and analysis of the composition and 
terminal nucleotide data of a set of "sequential subsets" of the Matrix 
Method of Analysis. Changing the filters which process the "string" 
corresponds to changing the "guess" of the "known" nucleotides, 
35 K l K 2 K 3 K 4 ■ > K n . 9 as well as the "unknown" nucleotides, X X ,X 2 ,X 3 ,X 4 of the 
Matrix Method of Analysis as applied to DNA sequencing. The order of 



WO 99/34322 



PCT/US98/27624 



84 



1 0 



1 5 



20 



the "groups of SFCs" of the "string" is established when "associations" with 
the "High Level Memory" are achieved for a given set of delayed Gaussian 
filters (i.e. the order of Fourier series representing information is solved 
when internal consistence is achieved according to the MM A). Then each 
Fourier series of the ordered "string" is recorded to the "High Level 
Memory" wherein each Fourier series of the ordered "string" may be 
multiplied by the Fourier transform of the delayed Gaussian filter 



correct order to form the ordered "string". 

Also, multiple other cascades of association "stages" ("association 
ensembles") may act as delay-line memory actuators that produce a time 
delay, <5(/w 0 ), during independent "activation" of a given "association 

ensemble" with recall from memory. In k, co -space, the time delay is 
equivalent to a modulation of the correlation function given by Eq. 
(37.63) corresponding to the independent phase shifts, S s , of the 

correlation function (Eq. (37.106)) of the separate "associated" "groups of 
SFCs". During "string" ordering by the Matrix Method of Analysis, the 
independent phase shifts, 8 S , may modify the order of the Fourier series 

of the "string" representing information. In addition, the independent 
phase shifts, S s , may initially modify the content of the "string" by 

altering the correlation function (Eq. (37.106)) to cause information to be 
"associated" which otherwise would not likely be or inhibit the 
"association" of information which otherwise would be. These 
mechanisms further provide for information with novel conceptual 
content. 



represented by 




that established the 
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SUB-APPENDIX V 

GENOMIC DNA SEQUENCING METHOD/MATRIX METHOD OF ANALYSIS 

5 Abstract 

As an overview, the Mills method of sequencing DNA comprises the steps of: 

a) preparing from segments of a DNA strand to be sequenced, families of 
polynucleotides, each family including all polynucleotides, complementary to at least a portion of 
the DNA segment and at least a portion of the 3' flanking DNA segment of the DNA strand to be 

1 0 sequenced, of the formula: 

K n , ■ • ■ K^K^X^X^ ■ X n 
ranging in length from K } X Y to K n . - X n wherein K X K 2 K Z K 4 - -K n . represents the nucleotides 5' 
to an internal reference point, the reference point defined as the dividing line between K x and X x ; 
wherein X l X 2 X 3 X 4 - •■X n represents the nucleotides 3' to the internal reference point; wherein n 

1 5 and ri are integers and n + ri , the number of nucleotides in a polynucleotide, is less than or 

equal to the number of nucleotides in a polynucleotide of length within the analyzable limit of the 
method for determining base composition and identity of the 3' terminal nucleotide of a 
polynucleotide; and wherein each polynucleotide in the family conforms to the criterion that if the 
polynucleotide contains X n it also contains X n _ x , X n __ 2 * -X x \ and the criteria that if the 

2 0 polynucleotide contains K n it also contains K n ^ x ,K n ._ 2 -K x \ 

b) determining the base composition (the number of As, T's, C's, and G's) and the 
identity of the 3 1 terminal base of each polynucleotide of each family; 

c) determining the base sequence of the longest polynucleotide in each family from 
the determined base composition and identity of the 3' terminal base of each polynucleotide in the 

2 5 family and derived change in base composition and terminal base between polynucleotides in 

each family; and 

d) determining the base sequence of the entire DNA strand to be sequenced based 
upon the overlapping sequences of the longest polynucleotides in each family. 

3 0 The base sequence of the longest polynucleotide of each family is determined by the Matrix 

Method of Analysis of the base composition of each polynucleotide in the family and the identity 
of the 3 f terminal base of each polynucleotide. 

The base sequence of the longest polynucleotide in each set is determined by: 
3 5 a) setting up a matrix consisting of -^-M + 1 columns and ~M rows where M is the 

number of nucleotides in the longest polynucleotide of the set; 
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b) assigning the longest polynucleotide a coordinate position in the matrix of column 
1, row 1; 

c) assigning polynucleotides which are successively one nucleotide shorter on the 5' 
end to each column position and polynucleotides which are successively one nucleotide shorter 

5 on the 3' end to each row position; 

d) determining all paths through the matrix from position 1,1 to position + 1, 

~ M which are consistent with the base composition and the 3' terminal base of the 

polynucleotide assigned to each position in the matrix and with the change in base composition 
and 3 f terminal base between polynucleotides; and 
10 e) from position + 1 , — M determining the path back to position 1,1 which 

permits the assignment of specific bases at each step either the 5' or 3' end of a polynucleotide, 
consistent with the compositional and terminal base data, to arrive at the sequence of the longest 
polynucleotide wherein the K^K^^K^ - -K n . is guessed and steps d) and e) are performed 

reiteratively until a sequence can be assigned without contradiction. 

1 5 

Mills [3, 4] has developed a method of determining the nucleotide 
sequence of a DNA molecule of arbitrary length as a single procedure by 
sequencing portions of the molecule in a fashion such that the sequence 
of the 5 1 end of the succeeding contiguous portion is sequenced as the 3' 

2 0 end of its preceding portion is sequenced, for all portions, where the 

order of contiguous portions is determined by the sequence of the DNA 
molecule. Sequencing of the individual portions is accomplished by 
generating a family of polynucleotides under conditions which determine 
that the elements are partial copies of the portion and are of random 

2 5 nucleotide length on the 3' and 5' ends about a dinucleotide which is an 

internal reference point; determining the base composition and terminal 
base identity of each element of the family and solving for the sequence 
by a method of analysis wherein the base composition and terminal base 
data of each element is used to solve for a single base of the sequence by 

3 0 assigning the base to either the 5' or 3' end of the partial sequence about 

the internal reference point as the entire sequence of the portion is built 
up from a dinucleotide. 

The molecules generated from the DNA to be sequenced comprise 
families of polynucleotides. Each family corresponds to a segment of the 
3 5 DNA to be sequenced and is made up of a longest polynucleotide (the 
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length of which is selected to be within the analyzable limit of the 
procedure used to determine base composition and identity of the 
terminal base) and shorter polynucleotides which form a "sequential 
subset" of the longest polynucleotides. Grouped hierarchically from the 
5 longest to the shortest polynucleotide of the family is progressively one 
nucleotide shorter than the preceding polynucleotide and has the same 
sequence except that it lacks the one nucleotide. A further restraint on 
the elements of the family is that there is a specific dinucleotide of the 
sequence contained in each element. The molecules can be envisioned as 
1 0 being built around an "axis" which is at the mid position of the common 
dinucleotide. The "axis" constitutes an internal reference point. The 
polynucleotides vary around the "axis," each containing one less 
nucleotide on the 3' or 5' end than its longer predecessor in the group. 
All such molecules are included in the family, from the longest to the 

1 5 shortest, a dinucleotide. 

The sequence of the DNA portion from which each family of 
polynucleotides has been made can be solved by determining the base 
composition (the number of A ! s, T's, C's, and G's) and the identity of the 3 ! 
terminal base of each polynucleotide of the family. The composition and 

2 0 terminal nucleotide data of the elements of each family of 

polynucleotides are used to solve the sequence of the corresponding DNA 
portion template by a method of first generating all polynucleotides 
which can be obtained from a guessed solution of the sequence by 
successive removal of a 3' or 5' nucleotide consistent with the data of the 

2 5 change in composition between set elements and with the further 

constraint that a specific dinucleotide of the sequence must be present in 
all polynucleotides. The terminal nucleotide data is used to determine if 
a subset of the hypothetical family of polynucleotides exists such that the 
elements have a one to one correspondence with the data of terminal 

3 0 nucleotide as well as composition. If no such subset exists, the process is 

repeated for improved guesses until convergence to the correct solution 
for the sequence occurs. 

A method which performs this analysis by testing for the validity 
of a guess for part of the sequence while solving for the remaining part 
3 5 using the composition and terminal base data independently to execute 

binary hypothesis testing decisions compatible with computer logic is the 
matrix method of analysis algorithm. 
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The matrix method of analysis is analogous to solving a system of n 
equations in n unknowns where the knowns are: 1) the structural 
properties of the polynucleotides, 2) the base composition and the 
identity of the terminal base, 3) the change in composition and change in 
5 terminal base between a polynucleotide and the next in the family. The 
method exploits the given information by implementing a reiterative 
procedure to find a path through a matrix of the possible polynucleotides 
having sequences consistent with the data. Final assignment of the 
sequence is made when the entire path finding procedure can be 
1 0 accomplished without contradictions between sequence assignment and 
actual data. 

Strategy of the Sequencing Method 

The strategy is to create a group of molecules which contain a 
reference point which is internal. Initially, location of the reference point 
is unknown, but it exists in all of the molecules. The molecules are a 
family of polynucleotides comprising complementary copies of a portion 
of the parent molecule from which they are generated and superimpose 
on the parent by alignment of this internal point of reference. The 
location of the point of reference or "axis," and the sequence of the parent 
molecule is solved for simultaneously by an algorithm called the matrix 
method of analysis. 

The family of polynucleotides can be thought of as being all 
molecules which result from the sequential loss of nucleotides from the 5' 
and 3' end of the longest polynucleotide of the group. An ordered 
pattern of terminal nucleotide change and nucleotide compositional 
change occurs between members of sequential subsets. This algorithm 
exploits the pattern of ordered systematic nucleotide compositional 
change and terminal nucleotide change that a designated longest 
polynucleotide with a given internal reference point and given nucleotide 
loss constraints can produce. 

Criteria of Polynucleotides 

The nucleotide sequence of a DNA strand can be solved by 
3 5 generating a family of polynucleotides overlapping portions of the DNA to 
be sequenced. Each family of polynucleotides forms a "sequential subset" 
of the longest polynucleotide of the group. The molecules are identical 
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less one nucleotide from either the 5' or 3' end of a given molecule, and 
the former are defined as sequential subsets of the latter. 

The molecules can be depicted as follows: 
K n . . . . K 4 K 3 K 2 K x X x X 2 X 3 X 4 . . . X n 
Where the series K x , K 2 , K 3 ,K 4 ... K n . represent the nucleotides of the 
polynucleotide 5* to the internal reference point, or axis, and the series 
X l9 X 29 X 3 ,X 4 --.X n represents the nucleotides of the polynucleotide on the 
3' side of the axis. The 5' end with respect to the axis is designated as 
the "known" portion of the molecules (this does not necessarily imply 
that this sequence is initially known), and the 3* end of the 
polynucleotide is designated as the "unknown" portion . Thus, 
K X ,K 2 ,K 3 ,K 4 -.. represent the "known" sequence and X X ,X 2 ,X 3J X 4 ••• 
represent the "unknown" sequence. The distinction is that in the 
matrix, as described below K X ,K 2 ,K 3 ,K 4 -.- appear as nucleotides, where 
as the X's represent variables. The nucleotides of the "known" portion 
can be known extrinsically or they can be guessed. 

The polynucleotides are governed by the following constraints. No 
polynucleotide contains X 2 without containing X x . In general terms, no 
polynucleotide contains X n without containing X n _ ly X tt _ 2 ,~ In 
addition, no polynucleotide contains K 2 without containing K x . That is, 
polynucleotide contains an unknown with out containing all preceding 
unknowns and, every polynucleotide contains all succeeding knowns if 
it contains any given known. As a set, all the polynucleotides satisfy 
these criteria and vary randomly at the 3' and 5' ends. 
The criteria can be represented symbolically as follows: 

X n ^X x (X n implies X,) 
30 

K H -*K { {K n < implies K x ) 

•'* K n'- X n'~ ( th e polynucleotides are random at the 5' and 3*; ends; 
the knowns and unknowns are variables where K = Known, X = 
3 5 Unknown, «* = 1 to 4**- and n = lro4---) 
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Principles of Matrix Method of Analysis 

The matrix method of analysis entails setting up a rectangular 
matrix where the designated longest polynucleotide appears at position 
(1,1). The sequence of one half of this molecule is "known". The 
nucleotide sequence at the other one half of the molecule is designated 
"unknown" and is represented by variables. The term "known" does not 
necessarily imply that the nucleotide sequence of the parent molecule is 
known initially. The division between the "knowns" and "unknowns" is 
the internal reference point. The location of the internal reference point 
is not necessarily known initially and can be changed by changing the 
knowns so that this sequence superimposes a different region of the 
parent molecule. That is, when the sequence is solved, it will 
superimpose a region of the parent and the location of the internal 
reference point will be fixed. The location on the parent is at the line 
dividing the "knowns" and the "unknowns". If the 5 1 end of the sequence 
(and consequently the entire sequence) superimposes on a different 
region of the parent, the location of the internal reference point would be 
different. Thus, the location of the internal reference point relative to the 
parent molecule is determined by the "knowns". 

An exemplary matrix is shown below for polynucleotides which 
conform to the criteria set forth. For a designated longest polynucleotide 
which contains a total of eight (8) nucleotides the matrix consists of 5 
rows and 4 columns. 



The matrix columns contain polynucleotides which have lost 
nucleotides at the 5' end; the rows are formed of polynucleotides which 
have lost nucleotides from the 3' end. Nucleotides are lost from the 5' 
end down any column and lost from the 3' end across any row. The 
matrix is constructed such that all the constraints governing the 




X2 X^ 



K 4 K j K-2 K\ X\ X 2 ^3 ^4 

K 3 K.2 K\ X± %2 X^ X 4 
K-2 X^ X-2 X-^ X 4 
K l X l X 2 X 3 X 4 
X x X 2 XyX At 




K 4 K^K 2 K l X l 
K 2 K 2 K l X l 

K.X. 
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polynucleotides are satisfied, and all possible polynucleotides are 
recorded in the matrix according to the describe format. 

The determination of the sequence of the polynucleotides proceeds 
as follows: starting at position (1,1) in the matrix, the base which has 
5 been lost is determined by the difference in base composition between 
the longest polynucleotide and the next longest of the set. The change is 
consistent with a move to position (1,2) and/or (2,1) of the matrix. The 
step is repeated for each polynucleotide of the family. These moves are 
down a column and/or across the row from left to right. Moves down a 

1 0 column or across a row from left to right are designated from/to moves. 
The result can be recorded, e.g. in a "lattice" which contains all coordinate 
positions arranged in levels such that each successive level from top to 
bottom corresponds to all possible from/to moves, and each successive 
level from bottom to top corresponds to all possible to/from moves. A 

1 5 to/from move is a movement up a column and/or across a row from right 
to left. 

General Lattice 
II 

/ \ 
12 21 

/ \ / \ 
13 22 31 

/ \ / \ / \ 
14 23 32 41 

/\/\/\/\ 
K4K3K2K1 24 33 42 51 

/ \ / \ / \ / 
K3K2K1 34 43 52 

/ \ / \ / 
K 2 Kj 44 53 
/ \ / 
Ki 54 

Lattice 

Polynucleotide Coordinate position 



K4K3K2K1 15 

K3K2K1 25 

K2K1 35 

Kj 45 



For each step, the base which could have been lost from the 3 T or 5 f 
2 0 end is determined, and the appropriate move to a position in the matrix 
is made. This establishes the appropriate path in the matrix which can 
be designated by connecting the corresponding coordinates in the lattice. 
This procedure is repeated until all consistent from/to moves are 
recorded in the lattice. At least one path is formed from coordinate 
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position (1,1) to a point of convergence, i. e., a coordinate position from 
which no further from/to moves can be made. 

The next step is to determine which path is the correct path. This is 
accomplished by starting at a point of convergence and determining 
5 which to/from steps for all single or binary decisions are consistent with 
the terminal base data as moves are made back to position (1,1) from the 
point of convergence. Assignment of a base to the 3* or 5' end is made by 
a to/from move which does not contradict the change in base. For all 
to/from moves, if the path that is chosen from one coordinate to another 
1 0 corresponds to a move across a row from right to left, then the base is 

assigned to the 3 1 end which is consistent with the move. That is the base 
change determined from the data occurred from the 3' end. A 
contradiction arises if this assignment is inconsistent with terminal base 
data for the polynucleotide represented at the coordinate position or if 

1 5 the change in terminal base for this step is inconsistent with the data. 

For all to/from moves, if the path that is chosen from one coordinate to 
another corresponds to a move up a column then the base change for that 
step indicates which base to assign to the 5' end. A contradiction would 
arise if the next "known" up the column in the matrix is different from 

2 0 that indicated by the base change. 

The sequence is solved when at least one path is found from (1,1) 
to a point of convergence by from/to moves and to the (1,1) position 
from the point of convergence by to/from moves at each data step 
without contradictions. The matrix method of analysis yields a unique 

2 5 solution for a matrix of all possible polynucleotides of size (^M + l, ^A/) 

that conform to the constraints for polynucleotides, for any set of data of 
M-l polynucleotides that are successively one nucleotide less and are 
sequential subsets from M-l nucleotides to a dinucleotide. (The longest 
polynucleotide is M nucleotides in length.) 

3 0 The key to the matrix method of analyze is that there is 

convergence to at least one of the terminal possibilities (point in the 
matrix at which no further from/to moves can be made). It may 
converge to more than one (e.g., if the sequence contained only A, or T, or 
C, or G bases, then it would converge to all possible termini of the matrix 
3 5 that yields the solution of the sequence). Once any terminus is 

determined to be correct, it can serve as an initiation point, that is, a 
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point, or coordinate position from which the initial to/from move is made. 
A terminus representing a single nucleotide or single variable in the 
matrix is correct if it is consistent with the data. The sequence can be 
deciphered by making decisions at branch points and by taking the 
5 return path that is determined to be correct by the data, i. e. the terminal 
base and the change in the terminal base at each step. If more than one 
path is correct, anyone of the correct paths will yield the sequence. 



1 0 



1 5 



Examples of Solving Sequences bv the Matrix Method of Analysis 

To further illustrate the matrix method of determining sequence, 

examples of its application are given below. In each example a matrix for 

a polynucleotide family of eight nucleotides in length is shown. The 

lattice diagram shows all possible matrix from/to moves consistent with 

the change in composition data. The column labeled "path" represents 

the possible to/from moves in the matrix which are consistent with the 

terminal base data and the change in terminal base. The path which 

determines the solution to the sequence is read from bottom to top. 
EXAMPLE 1 

1 4 6 7 5 3 2 
ATTCGCTA 
X 1X3X3X4 



1. 

2. 
3. 
4. 
5. 



I 

ATTCX 1X2X3X4 
TTCX ^2X3X4 
TCX 1X2X3X4 
CXiX 2 X 3 X4 
X 1X2X3X4 



ATTCXjX2X3 
TTCX I X 2 X 3 
TCX 1X2X3 
CX 1X2X3 
X1X2X3 



ATTCX1X2 
TTCX1X2 
TCX iX 2 
CXiX 2 

X1X2 



ATTCXi 
TTCXi 
TCX 1 
CXi 
Xi 



Lattice 



Composition 




Terminal 






Data 


A 


Nucleotide 


Path 


Sequence 


3T.2C.IG.2A 


A 


A 


1,1 


ATTCGCTA 


3T.2C.1G.1A 


A 


A 


2,1 


TTCGCTA 


3T.2C.1G 


T 


T 


2,2 


TTCGCT 


2T.2C1G 




C 


2,3 


TTCGC 




T 






IT,2C,1G 




C 


3,3 


TCGC 




C 






ITJCtG 


T 


G 


3,4 


TCG 


1C.1G 


C 


G 


4,4 


CG 


1G 




G 


5.4 


G 



11 

/ \ 

12 21 

/ \ / 
13 22 
/ / \ 

14 23 32 

/ / \ / \ 

ATTC 24 33 42 

/ / / \ 

TTC 34 43 52 

/ \ / / 
TC 44 53 

/ \ / 

C 54 
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EXAMPLE 2 



1 2 5 7 6 4 3 
X^XjX* 





1 


2 


3 


4 


1. 


AGTCX1X2X3X4 


AGTCX1X2X3 


AGTCX1X2 


AGTCX i 


2: 


GTCXiX 2 X3X4 


GTCX t X2X3 


GTCX1X2 


GTCXj 


3. 


TCX 1X2X3X4 


TCX1X2X3 


TCX1X2 


TCXi 


4. 


CX1X2X3X4 


CX1X2X3 


CX1X2 


CXi 


5. 


X 1X2X3X4 


X1X2X3 


X1X2 


Xi 




Composition Terminal 








Lattice Data 


A Nucleotide 


Path 


Sequence 



It 

/ \ 

12 21 

/ / \ 

13 22 31 

/ / \ / 

14 23 32 

/ / / \ 

AGTC 24 33 42 

/ / \ / 

GTC 34 43 
/ \ / 
TC 44 

/ \ 
C 54 



3T,2G,1C,2A 




G 


1,1 1,1 


AGTCATTG 


AGTCATTG 




A 










3T,2G.IC!A 




G 


2,1 2,1 


GTCATTG 


GTCATTG 




G 










3T,1G,1C,IA 




G 


3.1 3,1 


TCATTG 


TCATTG 




G 










3T,1C,1A 




T 


3,2 3,2 


TCATT 


TCATT 


2T.1C1A 


T 
T 


T 


3,3 4,2 


TCAT 


CATT 


1T,1C,1A 




T 


4,3 4,3 


CAT 


CAT 


T 










1C,1A 




A 


4,4 4,4 


CA 


CA 




C 










1A 




A 


5,4 5,4 


A 


A 
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SUB-APPENDIX VI 
Input Context 

An Input Layer receives data and transforms it into a Fourier series 
5 in k y co -space wherein input context is encoded in time as delays which 
corresponds to modulation of the Fourier series at corresponding 
frequencies. The Fourier series in Fourier space represents information 
parameterized according to the data and the input context. The 
information is the data and the input context. The information is based 
1 0 on physical characteristics or representations of physical characteristics 
and physical context. Data from transducers responding to an input 
signal representative of the physical characteristics and the physical 
context is used to parameterize the Fourier series in k, co -space whereby 

i.) "Data" such as intensity and rate of change recorded by a 

1 5 transducer is represented in terms of the parameters p Qm and N of each 

component of the Fourier series wherein the input context corresponds to 
the physical context based upon the identity of a specific transducer and 
its particular transducer elements. The physical context maps on a one to 
one basis to the input context. The processed signals from each 

2 0 transducer which can be input from the Input Layer to other layers such 

as the Association Layer and the "String" Ordering Layer, and the 
Predominant Configuration Layer comprises a Fourier series as given by 
Eq. (37.33) and Eq. (37.33a) wherein: 

each of the factors and N of the Fourier series component is 

2 5 proportional to the rate of change of the signal response of each 

transducer which is proportional to the rate of change of the physical 
signal such as the surface roughness, or the intensity of sound, light, or 
temperature; and 

each of the factors p 0 ^ and z 0m of each Fourier component is 

3 0 inversely proportional to the amplitude of the signal response of each 

transducer which is proportional to the physical signal such as the 
surface roughness, or the intensity of sound, light, or temperature; or 

each of the factors N m ^ and N of the Fourier series component is 

proportional to the amplitude of the signal response of each transducer 
3 5 which is proportional to the physical signal such as the surface roughness, 
or the intensity of sound, light, or temperature; and 
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each of the factors p 0m and z 0m of each Fourier component is 
inversely proportional to the rate of change of the signal response of each 
transducer which is proportional to the rate of change of the physical 
signal such as the surface roughness, or the intensity of sound, light, or 
5 temperature; or 

each of the factors N m ^ and N of the Fourier series component is 

proportional to the duration of the signal response of each transducer; 
and 

each of the factors p 0 ^ and z 0m of each Fourier component is 
1 0 inversely proportional to the amplitude of the signal response of each 
transducer which is proportional to the physical signal such as the 
surface roughness, or the intensity of sound, light, or temperature. 

ii.) The input from the Input Layer to other layers shown in FIGURE 
21 can be an analog waveform in the analog case and a matrix in the 

1 5 digital case. Input context of a given transducer can be encoded in time 

as delays which correspond to modulation of the Fourier series in 
K co -space at corresponding frequencies whereby the data corresponding 
to each transducer maps to a distinct memory location called a "register" 
that encodes the input context by recording the data to corresponding 

2 0 specific time intervals of a time division structured memory. The input 

context maps on a one to one basis to an Input Layer section of a 
memory. Thus, there is a one to one map of physical context to input 
context to Input Layer section of a memory. The representation of 
information as a Fourier series in Fourier space allows for the mapping. 

2 5 iii.) Input context of a complex transducer system can be encoded 

in time by the mapping of data from the components of the transducer 
system to a memory structured according to a corresponding hierarchical 
set of time intervals representative of each transducer system with 
respect to different transducer systems, a transducer element's rank 

3 0 relationship relative to other transducer elements, and the response of a 

transducer element as a function of time. In terms of digital processing, 
the data from a transducer having n levels of subcomponents is assigned 
a master time interval with n + 1 sub time intervals in a hierarchical 
manner wherein the data stream from the final n th level transducer 
3 5 element is recorded as a function of time in the n + \ th time coded 

memory buffer. During processing the time intervals represent time 



WO 99/34322 



PCT/US98/27624 



97 

delays which are transformed into modulation frequencies which encode 
input context. FIGURE 3 is a flow diagram of an exemplary transducer 
data structure of a time delay interval subdivision hierarchy wherein the 
data from a transducer having n levels of subcomponents numbering 
5 integer m per level is assigned a master time interval with n + l sub time 
intervals in a hierarchical manner wherein the data stream from the final 
n th level transducer element is recorded as a function of time in the n + l 
th time coded sub memory buffer in accordance with the present 
invention. 

1 0 The "processor" may be taught the relationship between the "data" 

such as intensity and rate of change recorded by a transducer and the 
parameters such as p 0 and N of each component of the Fourier series 

by inputting standard physical signals to each transducer together with 
other information that is "associated" with the standards. A data base 

1 5 may be established. The information that is "associated" with the 

standard may be recalled and can comprise input into the Association 
Layer and the "String" Ordering Layer shown in FIGURE 20 and FIGURE 
21 during "processing" according to the present invention. 

The process of storing output from multiple transducers to memory 

2 0 further comprises creation of "transducer strings". In one embodiment, 

associations occur at the transducer level, and "SFCs" are mapped to 
distinct memory "registers" from the corresponding distinct transducers 
responding simultaneously, for example. In one embodiment, two or 
more Fourier series such as two or more "SFCs" of the "string" are "linked" 

2 5 whereby activation of any Fourier series such as a "SFCs" of the "string" 

may thereby activate other or all Fourier series of the "string" stored in 
the corresponding "registers". The activation may be based on 
probability. The activation probability may depend on the "strength of 
the linkage" which is defined in terms of a linkage probability parameter 

3 0 which increases with the linkage rate, the rate at which the activation of 

a Fourier series of a "string" thereby causes the activation of another 
Fourier series of the "string". Probability operators may activate other or 
all Fourier series of the "string" when any Fourier series of the "string" is 
"active" based on the linkage probability parameter. 
3 5 In a digital embodiment comprising "memory linkages" of the 

"transducer string", recalling any part of a "transducer string" from a 
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distinct memory location may thereby cause additional "linked" Fourier 
series of the "transducer string" to be recalled. In one embodiment, a 
linkage probability parameter is generated and stored in memory for 
each "string" Fourier series such as a "SFCs". A probability operand is 
5 generated having a value selected from a set of zero and one, based on 
the linkage probability parameter. If the value is one, the corresponding 
Fourier series is recalled. Thus, when any part of a "transducer string" is 
recalled from memory, any other "string" Fourier series is randomly 
recalled wherein the recalling is based on the linkage probability 
1 0 parameter. The linkage probability parameter is weighted based on the 
linkage rate. 



*(*)= ]x(f)e j2j *df 



Delay 



x(0= ]x{ty n ^dt 



(37.109) 



1 5 Consider a "transducer string" made up of multiple "groups of SFCs" 

where each "SFCs" represents information of the transducer system with 
respect to different transducer systems, a transducer element's rank 
relationship relative to other transducer elements, and the response of a 
transducer element as a function of time, space, or space and time. (The 

2 0 latter case applies to a transducer which is responsive to changes in the 

intensity of a parameter over time and spatial position). These aspects of 
each transducer are encoded via delays corresponding to modulation in 
k,co -space within a frequency band corresponding to each aspect of the 
transducer. 

2 5 The "string" in &,a)-space is analogous to a multidimensional 

Fourier series. The modulation within each frequency band may further 
encode context in a general sense. In one embodiment, it encodes 
temporal order, cause and effect relationships, size order, intensity order, 
before-after order, top-bottom order, left-right order, etc. which is 

3 0 relative to the transducer. 

Eq. (37.33a), the "read" total response in Fourier space 

comprising the superposition of M "FCs" wherein each "FC" corresponds to 
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the response of a "M or P element" with input context encoded by the 
modulation factor e Jkp ( P/b *->" +p '*-»>} becomes 



m=ln=-«1 j z_ V 



sin k "» ° m -n ^ 



20 



5 (37.110) 
where p, = v, r, is the modulation factor which corresponds to the 

physical time delay t, m and p^ = v^j^ is the modulation factor which 

corresponds to the specific transducer time delay t^, and are 

constants such as the signal propagation velocities.. 
1 0 "Associations" are established between Fourier series such as "SFCs" 

and "groups of SFCs" (i.e. a "string" is created) by "coupling" with 
Poissonian probability between "association ensembles" "carrying" the 
"SFCs" and "groups of SFCs". Input context is encoded by the transducer 
frequency band modulation factor e ~ Jk ^ p ^ +Pt ^) according to Eq. (37.110). 
15 In this case, Eq. (37.87b) is 
1 



exp— 



2 2 f 



+ t fh + r. 

\ z J 



(N„t n 



2"! 



(37.111a) 



And, Eq. (37.87c) is 



exp- 



2 2 
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The corresponding frequency difference angle, <p s , which follows from Eq. 



(37.89) is 



n 



<Ps = 



or, 



a 



_ v 



s m] =1 



N *Po. . P» m . P, m 



2v_ 



- + - 



V v <„„ J 



m, =1 



N m,Po m Pjb m Pi m 

2v„, 



V 



'mj y 



V 



V, 



2v„. 



mi y 



(37.112a) 

The corresponding frequency difference angle, <p s , which follows from Eq. 
(37.90) is 



71 



a. 



m. =1 



M < f N t 

2 An 



^ M < f N t 



A/. 



1 m^I 



(37.1 12b) 



Eq. (37.108), the "read" total response in Fourier space comprising 

the superposition of 5 "SFCs" wherein each "SFCs" corresponds to the 
1 0 response of M s "M or P elements", with input context encoded by the 

modulation factor e ~M^ + ^) 9 becomes the following "string". 



5 M s 



s—\ m=\ n=— <x> 



1 + ^ 



ff 



sin 



k n - n 



2k 



(f 



(37.113) 



sin 



K - n- 



2k 



Wo 



y 



where p, - v, /, is the modulation factor which corresponds to the 

's y m l s,m 's,m A 

physical time delay t t ^ m and P/b sm ~ v jb s J/b sm 1S modulation factor which 
1 5 corresponds to the specific transducer time delay . v, and v A are 

constants such as the signal propagation velocities. In another 
embodiment, the output is the Gaussian sampled and modulated 

"string" of Eq. (37.113) wherein each "SFCs" is multiplied by the Fourier 
transform of the delayed Gaussian filter (Eq. (37.50)) (i.e. the modulation 



2 0 factor e 



■^^A^L- j &^ 



) which gave rise to "coupling" 



and "association" to form the "string". V v is given by 
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1 5 



e -M'^ + O sin 



2k 



sin 



2;r 



(37.114) 



wherein input context is encoded by the modulation factor e 
Eq. (37.114) is also an exemplary "string" with each Fourier series 
5 multiplied by the Fourier transform of the delayed Gaussian filter 

t a k -tf'-r-T ->J?<^> 

represented by * ^ **J e v *° e 21 e ^° that established the 

correct order to form the ordered "string" given in SUB -APPENDIX IV- 
Ordering of Associations: Matrix Method. The index over s is 
independent of m since each "FC" of a given "SFCs" is filtered by the same 
1 0 Gaussian filter. In embodiments, the index for the Gaussian filter is not 
independent of m. In one case, some "FCs" may be filtered by the same 
Gaussian filters; whereas, other "FCs" may be filtered by different 
Gaussian filters. In another case, each "FC" may be filtered by a different 
Gaussian filter. 

For the case where v Smm t 0tM =p 0 ^ and k p = k t . t the "string" in Fourier 
space is one dimensional in terms of k p and is given by 



S M s 



SMl s-l m-\ n=-» 



(f 



sin 



vv 



2lZ 
Po 



'w. J 



(37.115) 

The "string" comprises a Fourier series, a linear sum of "FCs" each 
2 0 multiplied by its corresponding Gaussian filter modulation factor and 
modulation factor which encodes input context (Eqs. (37.114-37.115)). 
FIGURE 19 is a flow diagram of an exemplary hierarchical relationship of 
the signals in Fourier space comprising "FCs", "SFCs", "groups of SFCs", and 
a "string" in accordance with the present invention. Each "FC" is encoded 
by a "P element" or stored into and/or recalled from a "M element" as 
shown in FIGURE 18. 



25 
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SUB-APPENDIX VII 
Comparison of Processing Activity to Statistical 
Thermodynamics/Predominant Configuration 

5 The quantity of information that can be "associated" into ordered 

"strings" called M P strings" is essentially infinite based on the input to the 
layers of the "processor" comprising Fourier series in k,a> -space. 
Consider Eq. (37.33a). In the case that the parameter N sm spans 1 to 100, 
p 0jw spans 1 to 1000, and there are 1000 modulation bands, the number 
10 of distinct inputs W is 

W = 1000!1000!100! (37.116) 
Using Sterling's approximation 

In AM = N\nN-N (37.117) 
W is approximately 

15 ^ I2 ' 360 (37.1 18) 

In essence an infinite amount of information can be represented as 
distinct Fourier series in Kco -space according to this method of encoding 
it. 

According to statistical thermodynamics [20], a macroscopic 
2 0 thermodynamic system is viewed as an assembly of myriad 

submicroscopic entities in ever changing quantum states. Consider the 
number of distinct ways each called a microstate that a set number 
quanta of energy can be distributed between a set number of energy 
levels. The total number of microstates W associated with any 

2 5 configuration involving N distinguishable units is 

W =K§JF- <37119) 

where r) a represents the number of units assigned the same number of 
energy quanta (and, hence, occupying the same quantum number), and 
t] b represents the number of units occupying some other quantum level. 

3 0 As the number of units increases, the total number of microstates 

skyrockets to unimaginable magnitudes. Thus, one can calculate that an 
assembly of 1000 localized harmonic oscillators sharing 1000 energy 
quanta possesses more than 10 600 different microstates. This explosive 
expansion of the total number of microstates with increasing N is a direct 
3 5 consequence of the mathematics of permutations, from which arises also 
a second consequence of no less importance. Statistical analysis shows 
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that the emergence of a predominant configuration is characteristic of 
any assembly with a large number (N) of units. Of the immense total 
number of microstates that can be assumed by a large assembly, an 
overwhelming proportion arises from one comparatively, small set of 
5 configurations centered on, and only minutely different from, the 
predominant configuration— with which they share an empirically 
identical set of macroscopic properties. 

Eq. (37.119) is equally valid for the number of distinct ways that a 
set of "active" states at any given time can be distributed over N "P 
1 0 elements" and "M elements" where rj a represents one set of 

indistinguishable "P elements", or "M elements", and r} b represents 
another set of indistinguishable "P elements", or "M elements". Eq. 
(37.119) is equally valid for the number of distinct ways that 
interference or "coupling" can occur between cascades of association 

1 5 "stages" at any given time distributed over N "links" where rj a represents 

one set of indistinguishable "links", and 7] b represents another set of 
indistinguishable "links". Of the immense total number of microstates 
that can be assumed by a large assembly of "active" states distributed 
over a large set of "P or M elements" or by a large assembly of 

2 0 "couplings" distributed over many cascades of association "stages", an 

overwhelming proportion arises from one comparatively, small set of 
configurations centered on, and only minutely different from, the 
predominant configuration— with which they share an empirically 
identical set of macroscopic properties. Due to the large numbers of "P 

2 5 and M elements" and cascades of association "stages" involved in 

information processing, the present "processor's" performance is stable. 

Consider the "processor" on a component level such as that of a "P 
element". In an embodiment, the activation of a "P element" increases its 
excitability or probability of future activation with input. Each "P 

3 0 element" has an "activation" memory with a finite half-life. Repetitive 

"activation" of a "P element" results in a longer half-life of the increased 
excitability; thus, the "activation" memory becomes long term. The same 
principle applies to ensembles of association "stages", "processor 
elements" ("P elements"), and "memory elements" ("M elements") and 
3 5 "configurations" of "couplings" of ensembles. Each ensemble is comprised 
of "stages", "P elements", or "M elements" in different states of "activity" 
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where each state is equivalent to a microstate of statistical 
thermodynamics. A predominant configuration arises for any ensemble. 
Of the immense total number of microstates that can be assumed by an 
ensemble, an overwhelming proportion arises from one comparatively, 
small set of configurations centered on, and only minutely different from, 
the predominant configuration-with which they share an empirically 
identical set of macroscopic properties. On a higher level, a configuration 
of "couplings" between ensembles increases the activation of the "stages" 
"P elements", or "M elements" comprising the ensembles. Analogously to 
statistical thermodynamics, a predominant configuration arises from the 
ensemble level. Consider the "processor" on a higher level. The 
activation history of each ensemble relates to a hierarchical activation 
relationship of "coupled" ensembles which gives rise to a precedence of 
higher order predominant configurations. The ability to associate 
information and create novel information, is a consequence. Machine 
learning arises by the feedback loop of transducer input to the coupled 
predominant configurations which increases the basis for creating 
information with novel conceptual content. 
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Modifications and substitutions of the system elements and process 
1 0 steps by one of skill in the art are considered within the scope of the 

present invention which is not to be limited except by the claims. What 
is claimed is: 
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CLAIMS 

What is claimed is: 

1. A method for pattern recognition and processing of information, the 
5 information comprising data representative of physical characteristics or 
representations of physical characteristics within an input context of the 
characteristics, the method comprising the steps of: 

a.) encoding the data as parameters of at least two Fourier 
components in Fourier space; 
1 0 b -) adding the Fourier components to form at least two Fourier 

series in Fourier space, the Fourier series representing the information; 

c. ) sampling at least one of the Fourier series in Fourier space with 
a filter to form a sampled Fourier series; 

d. ) modulating the sampled Fourier series in Fourier space with the 

1 5 filter to form a modulated Fourier series; 

e. ) determining a spectral similarity between the modulated Fourier 
series and another Fourier series; 

f. ) determining a probability expectation value based on the 
spectral similarity; 

20 S-) generating a probability operand having a value selected from a 

set of zero and one, based on the probability expectation value; 

h. ) repeating steps a-g until the probability operand has a value of 

one; 

i. ) adding the modulated Fourier series and the another Fourier 

2 5 series to form a string of Fourier series in Fourier space, and 

j.) recording the string of Fourier series to memory. 

2. The method of claim 1, wherein each Fourier component of the 
Fourier series comprises a quantized amplitude, frequency, and phase 

3 0 angle. 

3. The method of claim 2, wherein the Fourier series in Fourier space 
is selected from one of: 



^ 4k 
2,2, JT'h^N 

m=I /!=-<» i ^ z 

k 2 
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and 



^ ^ 4tt 4 

Zu — ~ a * sin 



(f 



1 + 4Po„Zo„ 



K 



2k 
k -n 

Po m J 



sin 



2k 
k z — n 



wherein a 0m is a constant, k p and ^ are the frequency variables, re, in, and 
M are integers, and 7V„ % , , p 0 ^ , and z Qm are data parameters selected 

5 from one of: 

i. ) each of the data parameters N m and TV of the Fourier series 

component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p 0w and z 0m of each Fourier component 

1 0 is inversely proportional to the amplitude of the physical characteristic; 
or 

ii. ) each of the data parameters N m and N m of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 

1 5 each of the data parameters p Qm and z 0m of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic ; or 

iii. ) each of the data parameters N and N m of the Fourier series 

component is proportional to the duration of the signal response of each 

2 0 transducer; and 

each of the data parameters p 0w and z Qm of each Fourier component 
is inversely proportional to the physical characteristic. 

4. A method of claim 1 wherein the step b further comprises encoding 

2 5 the input context in time as a characteristic modulation frequency band 

in Fourier space of the Fourier series. 

5. The method of claim 4 wherein the characteristic modulation 
frequency band in Fourier space represents the input context according 

3 0 to at least one of a specific transducer or transducer element, and 

fundamental relationships including a temporal order, a cause and effect 
relationship, a size order, an intensity order, a before-after order, a top- 
bottom order, and a left-right order. 
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6. The method of claim 5 wherein the transducer has n levels of 
subcomponents and is assigned a master time interval with « + l sub time 
intervals in a hierarchical manner that parallels and corresponds to the n 
levels of the transducer subcomponents wherein the data stream from 

5 the n th level transducer sub component is recorded as a function of time 
in the az + 1 sub time interval and the time intervals represent time delays 
which correspond to the characteristic modulation frequency band in 
Fourier space which represents the input context according to the specific 
transducer or transducer subcomponent. 

1 0 

7. The method of claim 4 wherein the characteristic modulation 
having a frequency within the band in Fourier space is represented by 
e-^o w hich corresponds to the time delay 8{t-t 0 ) wherein / is the 
frequency variable, t is the time variable, and / 0 is the time delay. 

1 5 

8. A method of claim 7 wherein the step b further comprises storing 
the characteristic modulation frequency in a distinct memory location 
within the band encoded as a delay in time. 



2 0 9. The method of claim 7, wherein the Fourier series in Fourier space, 
has a characteristic modulation having a frequency within the band 
represented by e ~ Jk *( p *» +Pt ») and is selected from one of: 



XX 



4k 



a n N m AT e 



,2 u 0. 



sin 
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2kN„ 



■ — n- 



^ sin k. 



2kN 



25 



30 



and 

M 



E-cp 4k 
Zj 



N m p n 



■ — n- 



N m z Q 2kM \ 
sin k — n 



k 2 °" n 7 
k 2 

wherein p, =v, t t is the modulation factor which corresponds to the 

'm *m m 

physical time delay t , = t fi is the modulation factor which 
corresponds to the specific transducer time delay , v and are 
constants such as the signal propagation velocities, a 0 is a constant, k p 
and k z are the frequency variables, n y m, and M are integers, and N 9 
N m , p 0w , and z Qm are data parameters selected from one of: 
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1 0 



L) each of the data parameters N and N m of the Fourier series 

component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p 0n and z 0m of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

ii.) each of the data parameters N mpQ and N of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 



each of the data parameters p Qm and z K of each Fourier component 
is inversely proportional to the rate of change of the physical 
characteristic; or 

iii.) each of the data parameters N and TV of the Fourier series 

component is proportional to the duration of the signal response of each 
1 5 transducer; and 

each of the data parameters p 0 ^ and z 0m of each Fourier component 
is inversely proportional to the physical characteristic. 



10. The method of claim 7, wherein the string has a characteristic 
2 0 modulation having a frequency within the band represented by 
e~ Jkp ^ +p '"^ is selected from one of: 
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2 5 wherein Pt Stm =\Jt sm is tlie modulation factor which corresponds to the 
physical time delay /, , p = v t is the modulation factor which 
corresponds to the specific transducer time delay t- , v and v, are 
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constants such as the signal propagation velocities, a 0 is a constant, k p 
and k z are the frequency variables, n, m, s, M 5> and S are integers, and 
^s w m m y N s > p 0 w , and z 0sm are data parameters selected from one of: 

i.) each of the data parameters N sm ^ and N of the Fourier series 

5 component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p 0 ^ and z 0sjn of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

10 ii.) each of the data parameters N Km and of the Fourier series 

S ' m P0 s * m ^ 

component is proportional to the amplitude of the physical characteristic; 
and 

each of the data parameters p 0 ^ and Zo tm of each Fourier component 

is inversely proportional to the rate of change of the physical 

1 5 characteristic; or 

iii.) each of the data parameters N smpQ and N of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0am and z 0 of each Fourier component 

2 0 is inversely proportional to the physical characteristic. 



11. The method of claim 10, wherein each Fourier series of the string is 
multiplied by the Fourier transform of the delayed Gaussian filter 

, , 4 V ^1 4f*-r-T -;J?(^> 

represented by e ^ aspo) e i 1 ** e 2{ e Va *° wherein the filter 
2 5 established the association to form the string, wherein the string is 
represented by: 



S M t 



s—\ m=\n=-°° 



An 



•cu N„„ N. 
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2 0j.,« s >™po 



-jJS^) -if^T -^('-aj 



k p -n 



vv 



2k 



sin 



2k 



wherein v spo and v sz0 are constants such as the signal propagation 
velocities in the p and z directions, respectively, and are delay 
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parameters and a spo and a sz0 are half-width parameters of a 
corresponding Gaussian filter in the p and z directions, respectively, 
Pt sm ~ v t s Jt sm i s the modulation factor which corresponds to the physical 
time delay r^, Pjb sm =v jb s J/b sni is the modulation factor which corresponds 
5 to the specific transducer time delay t, , v, and v A are constants such 
as the signal propagation velocities, a 0 ^ is a constant, k p and k z are the 
frequency variables, «, m, j, M 5 , and 5 are integers, and AT , N rm , p ft , 
and are data parameters selected from one of: 

i.) each of the data parameters N m and TV of the Fourier series 

1 0 component is proportional to the rate of change of the physical 

characteristic; and 

each of the data parameters p 0 ^ and z 0jm of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

15 ii.) each of the data parameters TV and N m of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 

each of the data parameters p 0 ^ and Zq of each Fourier component 

is inversely proportional to the rate of change of the physical 

2 0 characteristic; or 

iii.) each of the data parameters and N m of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 ^ and Zq of each Fourier component 

2 5 is inversely proportional to the physical characteristic. 

12. A method of claim 1 wherein step b) further comprises creating 
transducer strings by obtaining a Fourier series from at least two selected 
transducers and adding the Fourier series. 

30 

13. A method of claim 12 further comprises selecting transducers that 
are active simultaneously. 

14. The method of claim 13 wherein the transducer string is stored in a 

3 5 distinct memory location wherein a characteristic modulation having a 
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frequency within the band in Fourier space is represented by e ~ i270 ° 
which corresponds to the time delay 5(r-r 0 ) wherein / is the frequency 
variable, t is the time variable, and t 0 is the time delay. 

5 15. The method of claim 14 wherein the step b) further comprises 
recalling any part of the transducer string from the distinct memory 
location which thereby causes additional Fourier series of the transducer 
string to be recalled. 

10 16. The method of claim 1, wherein the filter is a time delayed 
Gaussian filter in the time domain. 



17. A method of claim 16 wherein the Gaussian filter comprises a 
plurality of cascaded stages each stage having a decaying exponential 
1 5 system function between stages. 



18. A method of claim 17 wherein the Gaussian filter is modulated in 
the time domain to produce a frequency shift of the sampling and 
modulation in the frequency domain. 

20 

19. The method of claim 16 wherein the filter is characterized in time 
by: 

(-f)' 

—e ^ 
[N 

wherein J— is a delay parameter, a is a half-width parameter, and t is 
2 5 the time parameter. 



20. The method of claim 19 wherein the filter, in frequency space, is 
characterized by: 

3 0 wherein J— and a are a corresponding delay parameter and a half- 
width parameter in time, respectively, and / is the frequency parameter. 
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2 1 . The method of claim 1 wherein the probability expectation value is 
based upon Poissonian probability. 



1 0 



1 5 



20 



25 



22. The method of claim 21 wherein the probability expectation value 
is characterized by 

■cos20 5 



n 



cos(5 ir + 2sin0 s ) 



wherein P is the maximum probability of at least one other Fourier 
series being associated with a first Fourier series, p t is a probability of at 

least one other Fourier series being associated with a first Fourier series 
in the absence of coupling of the first Fourier series with the at least one 
other Fourier series, ji] is a number that represents the amplitude of 

spectral similarity between at least two filtered or unfiltered Fourier 
series, <j> 5 represents the frequency difference angle between at least two 
filtered or unfiltered Fourier series, and 5 S , is a phase factor. 

23. The method of claim 22 wherein fi 2 s is characterized by 
1 



A a =(8*r 



1 1 

at at 



[2k y a\ + a] 



X a o„,AU X X^ ex P- 



2 2 

at at 



at + a] 



a, 



a r 



2v„, 



2v„, 



and 



correspond to delay parameters of a first and s-th time 



delayed Gaussian filter, respectively, a x and a s corresponding half-width 

parameters of a first and s-th time delayed Gaussian filter, respectively, 
M| and M s are integers, <z 0 and a Qm are constants, v mj and v mj are 

constants such as the signal propagation velocities, and N„ h , N mg9 p 0 , and 

p 0 are data parameters selected from one of: 

i.) each of the data parameters N and N ms of the Fourier series 

component is proportional to the rate of change of the physical 

characteristic; and 
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each of the data parameters p 0 ^ and p 0 ^ of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

iL) each of the data parameters N nh and N ms of the Fourier series 

5 component is proportional to the amplitude of the physical characteristic; 
and 

each of the data parameters p^ and p 0 ^ of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

1 0 Hi.) each of the data parameters N mt and N mj of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 ^ and p 0 of each Fourier component 
is inversely proportional to the physical characteristic. 



1 5 



24. The method of claim 22 wherein <p s is characterized by 

2v„ 



r 

71 



a, a r ~, 2v_. At, 



m, 

and - — correspond to delay parameters of a first and s-th time 



delayed Gaussian filter, respectively, a x and a s corresponding half-width 
2 0 parameters of a first and s-th time delayed Gaussian filter, respectively, 
M l and M s are integers, a 0 and a Q are constants, v and v are 

constants such as the signal propagation velocities, and N , N m , p 0 , and 
p Qms are data parameters selected from one of: 

i.) each of the data parameters N nh and N ms of the Fourier series 
2 5 component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p 0 ^ and p 0 of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 



or 
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1 0 



ii. ) each of the data parameters N and N m of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 

each of the data parameters p 0 ^ and p 0m of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

iii. ) each of the data parameters N mi and N of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 ^ and p 0w of each Fourier component 
is inversely proportional to the physical characteristic. 



25. The method of claim 22 wherein f5 2 s is characterized by 
1 



2 2 

at cc 



P^W'l Y a 0 N Y a 0 N 



exp— 



a, 



a. 



2v„ 



K„Po m! } Pfi m> ^ P,., 

2v„ 



m s J 



1 5 wherein p, = v. t and p, = v r t t are the modulation factors which 

corresponds to the physical time delays t tm and t tm , respectively, 
Ptb ~ v fb t n> an d Pfb ~ v fb t ib are the modulation factors which 
corresponds to the specific transducer time delay and , 
respectively, v tm , v, w , , and are constants such as the signal 

2 0 propagation velocities, and correspond to delay parameters of a 

first and s-th time delayed Gaussian filter, respectively, a x and a s 

corresponding half-width parameters of a first and s-th time delayed 
Gaussian filter, respectively, M x and M s are integers, a Qm , a Qm are 

constants, v nh and v mj are constants such as the signal propagation 
2 5 velocities, and A^, W mj , p 0 , and p 0w are data parameters selected from 

one of: 
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1 0 



1 5 



20 



25 



30 



i. ) each of the data parameters N mi and N nh of the Fourier series 
component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p 0 ^ and p 0 ^ of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

ii. ) each of the data parameters N mi and N ms of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 

each of the data parameters p^ and p 0 ^ of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

iii. ) each of the data parameters N nh and N„ h of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 ^ and of each Fourier component 
is inversely proportional to the physical characteristic. 



26. The method of claim 22 wherein ^ is characterized by 



K 



a, 



a 



Wjii P ft 

2v., 



Pf» m . A„ 



2v_ 



"fin 



' m > )) 



a 



1 m,=l 



2v„ 



- + ■ 



mi ) 



and Pt m =v t l t are the modulation factors which 



wherein 

corresponds to the physical time delays t Sm and / , respectively, 

p K^ v P> mx t fl>m x and ^ = \^ ffl| are the modulation factors which 
corresponds to the specific transducer time delav / A and 

J An, ft>m s 9 

respectively, , v,^ , , and Vy& are constants such as the signal 



propagation velocities, 



a, 



and 



correspond to delay parameters of a 



first and s-th time delayed Gaussian filter, respectively, a x and a s 
corresponding half-width parameters of a first and s-th time delayed 
Gaussian filter, respectively, Af l9 and M s are integers, a n and a n are 
constants, v mi and v mj are constants such as the signal propagation 



WO 99/34322 



PCT/US98/27624 



118 

velocities, and N mt , 7V m ,p 0 ^, and p 0n are data parameters selected from 
one of: 

i. ) each of the data parameters N mi and N of the Fourier series 

component is proportional to the rate of change of the physical 
5 characteristic; and 

each of the data parameters p 0 and p 0 of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

ii. ) each of the data parameters N mi and N of the Fourier series 

1 0 component is proportional to the amplitude of the physical characteristic; 

and 

each of the data parameters p 0 ^ and p 0n of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

15 iii.) each of the data parameters N mi and N nh of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p Q ^ and p 0 of each Fourier component 
is inversely proportional to the physical characteristic. 

20 

27. A method of linking at least two Fourier series stored in a memory 
comprising the steps of 

a. ) generating a probability expectation value that recalling any 
part of one of the Fourier series from the memory causes at least another 

2 5 Fourier series to be recalled from the memory; 

b. ) storing the probability expectation value to memory; 

c. ) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value, and 

d. ) recalling the at least another Fourier series from the memory if 

3 0 the operand is one. 

28. The method of claim 27 whereby the probability expectation value 
increases with a rate of recalling any part of any of the Fourier series. 

3 5 29. A method for pattern recognition and processing information 
comprising the steps of: 
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representing the information as a plurality of Fourier series in 
Fourier space; 

forming associations between at least two of the Fourier series by 
modulating and sampling the Fourier series with filters and by coupling 
5 the filtered Fourier series based on a probability distribution; 

adding the associated Fourier series to form a string, and 
ordering the string. 

30. The method claim 29 wherein the filter is a time delayed Gaussian 
1 0 filter in the time domain. 

31. The method claim 29 wherein the probability distribution is 
Poissonian. 



1 5 32. The method of claim 1, wherein the string is selected from one of: 



S M s oo 

EEE 

5=1 m-\ «=-oo 



4k 



1 + 



k p -n 



2k 



sin 



k—n 



25 



and 

S M * 00 A TV 

V* 47ZT 

2, 2j — ~ a " 



-sin 



k p -n 



2k 



Po 



(f 



sin 



k—n 



2k 



'*«! J 



wherein a Qjm is a constant, k p and k z are the frequency variables, n, m, s, 
2 0 Af„ and S are integers, and A^, A^,p 0 ^, and ^ are data parameters 

selected from one of: 

i. ) each of the data parameters N sm and N m of the Fourier series 

component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p Q m and Zo tm of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

ii. ) each of the data parameters N sm and N tm of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
30 and 
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each of the data parameters p 0 ^ and z^ sm of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

iii.) each of the data parameters N s and N of the Fourier series 

5 component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 w and z 0sm of each Fourier component 
is inversely proportional to the physical characteristic. 

1 0 3 3. A method of ordering a string representing information using a 
high level memory comprising a set of initial ordered Fourier series, 
comprising the steps of: 

a. ) obtaining a string from a memory; 

b. ) selecting at least two filters from a selected set of filters; 

1 5 c.) sampling the string with the filters such that each of the filters 

produces a sampled Fourier series, each Fourier series comprising a 
subset of the string; 

d. ) modulating each of the sampled Fourier series in Fourier space 
with the corresponding selected filter such that each filter produces an 

2 0 order formatted Fourier series; 

e. ) adding the order formatted Fourier series produced by each 
filter to form a summed Fourier series in Fourier space; 

f. ) obtaining an ordered Fourier series from the memory; 

g. ) determining a spectral similarity between the summed Fourier 

2 5 series and the ordered Fourier series; 

h. ) determining a probability expectation value based on the 
spectral similarity; 

i. ) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value; 

3 0 j.) repeating steps b-i until the probability operand has a value of 

one; 

k.) storing the summed Fourier series to an intermediate memory; 

1.) removing the selected filters from the selected set of filters to 
form an updated set of filters; 
3 5 m.) removing the subsets from the string to obtain an updated 

string; 
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n.) selecting an updated filter from the updated set of filters; 

o.) sampling the updated string with the updated filter to form a 
sampled Fourier series comprising a subset of the string; 

p.) modulating the sampled Fourier series in Fourier space with the 
5 corresponding selected updated filter to form an updated order 
formatted Fourier series; 

q.) recalling the summed Fourier series from the intermediate 
memory; 

r.) adding the updated order formatted Fourier series to the 
1 0 summed Fourier series from the intermediate memory to form an 
updated summed Fourier series in Fourier space; 

s.) obtaining an updated ordered Fourier series from the high level 
memory; 

t.) determining a spectral similarity between the updated summed 

1 5 Fourier series and the updated ordered Fourier series; 

u.) determining a probability expectation value based on the 
spectral similarity; 

v.) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value; 
20 w 0 repeating steps n-v until the probability operand has a value of 

one or all of the updated filters have been selected from the updated set 
of filters; 

x.) if all of the updated filters have been selected before the 
probability operand has a value of one, then clearing the intermediate 

2 5 memory and returning to step b; 

y.) if the probability operand has a value of one, then storing the 
updated summed Fourier series to the intermediate memory; 

z.) repeating steps 1-y until one of the following set of conditions is 
satisfied: the updated set of filters is empty, or the remaining subsets of 

3 0 the string of step m.) is nil; and 

aa.) storing the Fourier series in the intermediate memory in the 
high level memory. 

34. A method of claim 33 wherein the high level memory is initialized 
3 5 with standard inputs. 
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35. The method of claim 33 wherein the ordering is according to one of 
the list of: temporal order, cause and effect relationships, size order, 
intensity order, before-after order, top-bottom order, or left-right order. 

5 36. The method of claim 33 wherein each filter of the set of filters is a 
time delayed Gaussian filter having a half-width parameter a which 
determines the amount of the string that is sampled. 

37. The method of claim 33 wherein each filter of the set of filters is a 

1 0 time delayed Gaussian filter having a delay parameter J— which 

V a 

corresponds to a time point. 



1 5 



3 8. The method of claim 37, wherein each Fourier series of the string is 
multiplied by the Fourier transform of the delayed Gaussian filter 



4 

represented by e ^ e e ^ ' a -* J e wherein the filter 

established the correct order to form the string, wherein the ordered 
string is represented by: 



5=1 m~\ n=~<» J _|_ 



k 2 a Q Jjm N s.m pa N s,m to e 



20 



25 



-j k p(pfi> sm +p lsm ) . 
e v ' sn} sm 



k p -n 



2n 



sin 



2n 



wherein v spo and v^ 0 are constants such as the signal propagation 
velocities in the p and z directions, respectively, J^ei and are delay 

parameters and a spo and a SZQ are half-width parameters of a 

corresponding Gaussian filter in the p and z directions, respectively, 
p t =v, t t is the modulation factor which corresponds to the physical 

time delay * w Pjb sm =v jb 5 Jjb sm 1S the modulation factor which corresponds 

to the specific transducer time delay r~ , v, and v, are constants such 

as the signal propagation velocities, a Qgm is a constant, k p and k z are the 

frequency variables, «, m, s, M s , and S are integers, and N „ , N , p n , 

and z Qjm are data parameters selected from one of: 
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i. ) each of the data parameters A^„ % and N of the Fourier series 

component is proportional to the rate of change of the physical 
characteristic; and 

each of the data parameters p 0 ^ and of each Fourier component 

5 is inversely proportional to the amplitude of the physical characteristic; 
or 

ii. ) each of the data parameters N sm ^ and N of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 

1 0 each of the data parameters p 0jM and Zq m of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

iii. ) each of the data parameters N sm and N m of the Fourier series 

component is proportional to the duration of the signal response of each 

1 5 transducer; and 

each of the data parameters p 0 ^ and Zo am of each Fourier component 

is inversely proportional to the physical characteristic. 

39. The method of claim 38, wherein v sm t 0 sm =p 0sjn and k p = k z such that 

2 0 the string in Fourier space is one dimensional in terms of k p and is 

represented by 



(f 



sin 



2k 



„ ^ 2 

wherein v spQ is a constant such as the signal propagation velocity in the p 
direction, J^~^ is a delay parameter and cc spo is a half- width parameter of 

2 5 a corresponding Gaussian filter in the -space, p, = v A /- is the 

modulation factor which corresponds to the specific transducer time 
delay t fb sm ^ v jb sm is a constant such as the signal propagation velocity, a 0 
is a constant, k p is the frequency variable, n, m 9 s, M sJ and S are integers, 
aR d N s. mpo anc * Po SJII are data parameters selected from one of: 

3 0 i.) the data parameter N s of the Fourier series component is 

proportional to the rate of change of the physical characteristic; and 
the data parameter p 0 ^ of each Fourier component is inversely 

proportional to the amplitude of the physical characteristic; or 
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ii. ) the data parameter of the Fourier series component is 

proportional to the amplitude of the physical characteristic; and 

the data parameter p 0 ^ of each Fourier component is inversely 

proportional to the rate of change of the physical characteristic; or 

iii. ) the data parameter N s ^ of the Fourier series component is 

proportional to the duration of the signal response of each transducer; 
and 

the data parameter p 0 of each Fourier component is inversely 
proportional to the physical characteristic. 

40. The method of claim 33 wherein the probability expectation value 
is based upon Poissonian probability. 



1 5 



20 



25 



41. The method of claim 40 wherein the probability expectation value 
is characterized by 

I1[pi m + (P - J-p[-^(l^^)]cos(^ + 2sin fc) 

wherein P is the maximum probability that at least one other Fourier 
series is active given that a first Fourier series is active, p r is a 

probability of a Fourier series becoming active in the absence of coupling 
from at least one other active Fourier series, p 2 s is a number that 
represents the amplitude of spectral similarity between at least two 
filtered or unfiltered Fourier series, </> s represents the frequency 

difference angle between at least two filtered or unfiltered Fourier series, 
and S s , is a phase factor. 

42. The method of claim 41 wherein p 2 s is characterized by 



#=(8*) 3 -L. J-^E a 0 N m Y a 0 N m 

V27T V ocf + a 2 *" 1 ^ mj 



exp— < 



2 2 

atat 



a, 2 + a? 



nil / 



2v„ 



A. 



y J» m 
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wherein p = v t and p, = v t, are the modulation factors which 
corresponds to the physical time delays t tm and / , respectively, 
P/b mi = v fb m *fb mx and Pjb ms ^ v /b m /jb mj are the modulation factors which 
corresponds to the specific transducer time delay /*. and t A 
5 respectively, v > v v , and v. are constants such as the signal 

propagation velocities, and correspond to delay parameters of a 

oc x cc s 

first and s-th time delayed Gaussian filter, respectively, a x and a s 
corresponding half-width parameters of a first and s-th time delayed 
Gaussian filter, respectively, M x and M s are integers, a 0 , a Q are 

1 0 constants, v niy and v fflj are constants such as the signal propagation 

velocities, and A^, 7V mj , p 0-i , and p 0 ^ are data parameters selected from 
one of: 

i. ) each of the data parameters N mi and N ntj of the Fourier series 
component is proportional to the rate of change of the physical 

1 5 characteristic; and 

each of the data parameters p 0 ^ and p 0m of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

ii. ) each of the data parameters N mi and N nh of the Fourier series 

2 0 component is proportional to the amplitude of the physical characteristic; 

and 

each of the data parameters p 0 ^ and p 0 of each Fourier component 

is inversely proportional to the rate of change of the physical 
characteristic; or 

25 iii-) each of the data parameters N mi and N ms of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 ^ and of each Fourier component 
is inversely proportional to the physical characteristic. 

30 

43. The method of claim 41 wherein <p s is characterized by 
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K 



<t>s=- 



a. 



— + z 

m, =1 



a 



2v_ 



Pfl> m . P, 



v. 



-X 

m 5 =1 



2v_ 



a 



1 m,=i 



— + + 

2V V rr V 



wherein p^ = v,^ r,^ and P,^ = v,^f are the modulation factors which 
corresponds to the physical time delays t tm and t, , respectively, 
Pjb m{ z=iv jb mi t jb mx and Pfi> mM =Vjb n tjb ms are the modulation factors which 
corresponds to the specific transducer time delay r A and / A , 
respectively, v , v , v , and v. are constants such as the signal 



propagation velocities, 



a, 



and 



or. 



correspond to delay parameters of a 



first and s-th time delayed Gaussian filter, respectively, a } and cc s 
corresponding half-width parameters of a first and s-th time delayed 
10 Gaussian filter, respectively, Af„ and M s are integers, a 0 ^ and a 0m are 

constants, v mi and v mj are constants such as the signal propagation 
velocities, and N mi , A^p^, and p 0 ^ are data parameters selected from 
one of: 

i.) each of the data parameters N mi and N mt of the Fourier series 

1 5 component is proportional to the rate of change of the physical 

characteristic; and 

each of the data parameters p 0 ^ and p Qm of each Fourier component 

is inversely proportional to the amplitude of the physical characteristic; 
or 

2 0 ii.) each of the data parameters N mi and iV m of the Fourier series 

component is proportional to the amplitude of the physical characteristic; 
and 

each of the data parameters and p 0 ^ of each Fourier component 

is inversely proportional to the rate of change of the physical 
2 5 characteristic; or 

iii.) each of the data parameters N nh and N nu of the Fourier series 

component is proportional to the duration of the signal response of each 
transducer; and 

each of the data parameters p 0 ^ and p 0m of each Fourier component 
is inversely proportional to the physical characteristic. 



30 
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44. A system for pattern recognition and processing comprising: 
an input layer that receives data representative of physical 

characteristics or representations of physical characteristics within an 
5 input context of the physical characteristics and transforms the data into 
a Fourier series in Fourier space wherein the input context is encoded in 
time as delays corresponding to modulation of the Fourier series at 
corresponding frequencies; 

a memory comprising a set of initial ordered Fourier series; 
1 0 an association layer that receives a plurality of the Fourier series in 

Fourier space from the memory, forms a string comprising a sum of 
Fourier series, and stores the string in memory; 

a string ordering layer that receives the string and at least one 
ordered Fourier series from the memory, orders the Fourier series 

1 5 contained in the string to form an ordered string, and stores the ordered 

string in memory; and 

a predominant configuration layer that receives multiple ordered 
strings from the memory, forms complex ordered strings from the 
ordered strings, stores the complex ordered strings to the memory, and 

2 0 activates the components of any of the layers of the system. 

45. A method for pattern recognition and processing of information 
comprising the steps of: 

a. ) generating an activation probability parameter based a prior 

2 5 activation probability parameter and a weighting based on an activation 

rate of the corresponding component; 

b. ) storing the activation probability parameter in memory; 

c. ) generating a probability operand having a value selected from a 
set of zero and one, based on the activation probability parameter; 

3 0 d.) activating any component of one or more of the group consisting 

of an input layer, an association layer, a string ordering layer, and a 
predominant configuration layer, the activation being based on the 
activation probability parameter, and 
e.) repeating steps a-d. 

35 

46. A computer-readable medium on which is stored a computer 
program for providing a method for pattern recognition and processing of 
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information, the information comprising data representative of physical 
characteristics or representations of physical characteristics within an 
input context of the characteristics, the computer program comprising 
instructions which, when executed by a computer, perform the steps of: 
5 a.) encoding the data as parameters of at least two Fourier 

components in Fourier space; 

b. ) adding the Fourier components to form at least two Fourier 
series in Fourier space, the Fourier series representing the information; 

c. ) sampling at least one of the Fourier series in Fourier space with 
10 a filter to form a sampled Fourier series; 

d. ) modulating the sampled Fourier series in Fourier space with the 
filter to form a modulated Fourier series; 

e. ) determining a spectral similarity between the modulated Fourier 
series and another Fourier series; 

1 5 f-) determining a probability expectation value based on the 

spectral similarity; 

g. ) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value; 

h. ) repeating steps a-g until the probability operand has a value of 

2 0 one; 

i.) adding the modulated Fourier series and the another Fourier 
series to form a string of Fourier series in Fourier space, and 
j.) recording the string of Fourier series to memory. 

2 5 47. A computer-readable medium on which is stored a computer 

program for providing a method of ordering a string representing 
information using a high level memory, the high level memory 
maintaining a set of initial ordered Fourier series, the computer program 
comprising instructions which, when executed by a computer, perform 

3 0 the steps of: 

a. ) obtaining a string from a memory; 

b. ) selecting at least two filters from a selected set of filters; 

c. ) sampling the string with the filters such that each of the filters 
produce a sampled Fourier series, each Fourier series comprising a subset 

3 5 of the string; 
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d. ) modulating each of the sampled Fourier series in Fourier space 
with the corresponding selected filter such that each filter produces an 
order formatted Fourier series; 

e. ) adding the order formatted Fourier series produced by each 
5 filter to form a summed Fourier series in Fourier space; 

f. ) obtaining an ordered Fourier series from the memory; 

g. ) determining a spectral similarity between the summed Fourier 
series and the ordered Fourier series; 

h. ) determining a probability expectation value based on the 
1 0 spectral similarity; 

i. ) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value; 

j.) repeating steps b-i until the probability operand has a value of 

one; 

1 5 k.) storing the summed Fourier series to an intermediate memory; 

1.) removing the selected filters from the selected set of filters to 
form an updated set of filters; 

m.) removing the subsets from the string to obtain an updated 
string; 

2 0 n.) selecting an updated filter from the updated set of filters; 

o.) sampling the updated string with the updated filter to form a 
sampled Fourier series comprising a subset of the string; 

p.) modulating the sampled Fourier series in Fourier space with the 
corresponding selected updated filter to form an updated order 

2 5 formatted Fourier series; 

q.) recalling the summed Fourier series from the intermediate 
memory; 

r.) adding the updated order formatted Fourier series to the 
summed Fourier series from the intermediate memory to form an 

3 0 updated summed Fourier series in Fourier space; 

s.) obtaining an updated ordered Fourier series from the high level 
memory; 

t.) determining a spectral similarity between the updated summed 
Fourier series and the updated ordered Fourier series; 
3 5 u.) determining a probability expectation value based on the 

spectral similarity; 
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v.) generating a probability operand having a value selected from a 
set of zero and one, based on the probability expectation value; 

w.) repeating steps n-v until the probability operand has a value of 
one or all of the updated filters have been selected from the updated set 
5 of filters; 

x.) if all of the updated filters have been selected before the 
probability operand has a value of one then clearing the intermediate 
memory and returning to step b; 

y.) if the probability operand has a value of one, then storing the 
1 0 updated summed Fourier series to the intermediate memory; 

z.) repeating steps 1-y until one of the following set of conditions is 
satisfied: the updated set of filters is empty, or the remaining subsets of 
the string of step m.) is nil; and 

aa.) storing the Fourier series in the intermediate memory in the 
1 5 high level memory. 

48. A computer-readable medium on which is stored a computer 
program for providing a method for forming complex ordered strings, the 
computer program comprising instructions which, when executed by a 
computer, perform the steps of: 

a. ) recording ordered strings to a high level memory; 

b. ) forming association between Fourier series of the ordered 
strings to form complex strings; 

c. ) ordering the Fourier series of the complex strings to form 
complex ordered strings, and 

d. ) storing the complex ordered strings to the high level memory. 

49. A computer-readable medium on which is stored a computer 
program for providing a method for forming a predominant 

3 0 configuration, the computer program comprising instructions which, 
when executed by a computer, perform the steps of: 

a.) generating an activation probability parameter based a prior 
activation probability parameter and a weighting based on an activation 
rate of the corresponding component; 
3 5 b.) storing the activation probability parameter in memory; 

c.) generating a probability operand having a value selected from a 
set of zero and one, based on the activation probability parameter; 
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d.) activating any component of one or more of the group consisting 
of an input layer, an association layer, a string ordering layer, and a 
predominant configuration layer, the activation being based on the 
activation probability parameter, and 
5 e.) repeating steps a-d. 

50. A computer program product for pattern recognition and processing 
of information for use in a general purpose computer including a central 
processing unit and a memory, the memory maintaining a set of initial 

1 0 ordered Fourier series, the computer program product comprising: 

a computer readable medium; 

program code means embodied in the computer readable medium, 
the program code means comprising: 

means for receiving data representative of physical characteristics 
15 or representations of physical characteristics within an input context of 
the physical characteristics and transforming the data into a Fourier 
series in Fourier space wherein the input context is encoded in time as 
delays corresponding to modulation of the Fourier series at corresponding 
frequencies; 

2 0 means for receiving a plurality of the Fourier series in Fourier 

space including at least one ordered Fourier series from the memory, 
forming a string comprising a sum of the Fourier series and storing the 
string in memory; 

means for retrieving the string from memory, ordering the Fourier 
2 5 series contained in the string to form an ordered string and storing the 
ordered string in memory; and 

means for retrieving multiple ordered strings from the memory, 
forming complex ordered strings from the ordered strings and storing the 
complex ordered strings to the memory. 
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